How to get the exact plot of optical ring resonator drop port output curve? which is similar to gaussian beam i.e. intensity is highest at centre and it gradually decreases at the corners on both ends.

e1=input('enter the value of ei1')

e2=input('enter the value of ei2')

w=1.5:0.01:1.6;

k1=0.25;

k2=0.25;

a=0.0005;

ne=3.19;

r=1.5*10^-6;

l=2*pi*r;

kn=(2*pi*ne*l)./w;

b=sqrt(1-k1);

c=sqrt(1-k2);

f=sqrt(k1-k2);

g=sqrt(-1)*kn*l;

ed=e2*((c-b*exp(-a*l)*exp(-g))./(1-b*c*exp(-a*l)*exp(-g)))+e1*((-f*exp(a*l/2)*exp(-g/2))./(1-b*c*exp(-a*l)*exp(-g)));

y=abs(ed);

z=y.^2;

plot(w,z)

[ABHISHEK KUMAR]

after running this matlab code i  am getting a straight line. but it not correct i need the plot of gaussian beam....