function [indigo] = indigo(r,sigma,t,S,T,eta,E,ti)
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
red1=red(r,sigma,t,S,T,E,eta,ti);
orange1=orange(r,sigma,t,S,T,E,eta,ti,red1);
yellow1=yellow(r,sigma,t,S,T,E,eta,ti,orange1);
green1=green(r,sigma,t,S,T,E,eta,ti,yellow1);
z=log(eta);
m=z^0.25;
square=m^2;
quart=m^4;
a=2*m*eta;
b=rem(a,1);
steps1=a-b;
nsteps1=steps1+1;
c=2*square*eta;
d=rem(c,1);
steps2=c-d;
nsteps2=steps2+1;
stepsize=1/eta;
steps3=ti^2;
nsteps3=steps3+1;
stepsizeti=1/ti;
w1=rem(sigma,stepsizeti);
sigma1=sigma-w1;
v=2/(sigma^2*ti);
w=rem(T,v);
T1=T-w;
k1=(0.5*ti*T1);
steps4=(T1*ti*sigma^2)/2;
nsteps4=steps4+1;
blue1=blue(r,sigma,t,S,T,E,eta,ti,green1);
for (j=1:nsteps1)
    f(j)=E*exp(-m^2+((j-1)/eta));
end
for (j=1:nsteps1)
    if f(j)<= S
        c(j)=j;
    else
        c(j)=c(j-1);
    end
end
w1=rem(t,v);
t1=t-w1;
t2=((sigma^2*ti*t1)/2);
t3=t2+1;
qu=c(nsteps1);
q=round(qu);
r=rem(t,v);
t2=t-r;
pu=(t2/v)+1;
p=round(pu);
indigo=blue1(3,p,q);
end

