function [orange] = orange(r,sigma,t,S,T,E,eta,ti,red)
%returns the ns fourier transform of the boundary, on [-m,m] with 
%2m[exp(m^4)]+1 steps to [-m^2,m^2] with 2m^2[exp(m^4)]+1 steps
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;
x=linspace(-m,m,nsteps1);
k=(2*r)/sigma^2;
xsq=linspace(-m^2,m^2,nsteps2);
for (j=1:nsteps1)
    g(j)=red(2,j);
end
for (j=1:nsteps2)
    s(j)=stepsize*(-square*eta+(j-1));
    theta(j,1)=stepsize*g(1)*exp(pi*i*m*s(j));
    for (k=2:nsteps1)
        l(k)=stepsize*(-m*eta+(k-1));
        theta(j,k)=theta(j,k-1)+stepsize*g(k)*exp(-pi*i*l(k)*s(j));
    end
end
for (j=1:nsteps2)
    hatg(j)=theta(j,nsteps1);
end
for (a=1:2)
        if (a==1)
            for (j=1:nsteps2)
        orange(a,j)=xsq(j);
            end
        else
            for (j=1:nsteps2)
           orange(a,j)=hatg(j);
            end
        end
end
end