Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: need matlab code to solve transcendental equation
Replies: 0

 Nidhi Posts: 19 Registered: 9/4/10
need matlab code to solve transcendental equation
Posted: Oct 1, 2013 11:54 AM

Hi, I am new to MATLAB and need help in solving an equation.
I m trying to solve a transcendental equation to find out values of propagation constant (pconst) for TE modes of the fiber, which is of the form

tan(kappa.*d1)=(gamma1./kappa)+(gamma3./kappa)./(1-((gamma1./kappa).*(gamma3./kappa)));

Please help me out to find out the pconst by appropriate MATLAB code for it. The constants and variables used in the equation are given below

B1=0.6961663; B2=0.4079426; B3=0.8974794;
c= 3e14;
l1=0.0684043; l2=0.1162414; l3=9.896161;
w1=2*pi*c/l1; w2=2*pi*c/l2; w3=2*pi*c/l3;
n3=1.321;
L=40e+3;
D=200;
NA=0.72;
wp=1.533e15;
t= 6.58e-15;
einf=3.5;
lambda=.4:0.001:1.5 ;
d1=1.6;
k0=2*pi./lambda;
w=2*pi*c./lambda;
em=einf-(wp^2./(w.^2+sqrt(-1).*w./t));
n1=sqrt(1+(B1*w1^2)./(w1^2-w.^2)+(B2*w2^2)./(w2^2-w.^2)+(B3*w3^2)./(w3^2-w.^2));
n2=sqrt(n1.^2-NA.^2);
kappa=sqrt(((w.^2./c^2).*em)-pconst.^2);
gamma1=sqrt(pconst.^2-((w.^2./c^2).*n1.^2));
gamma3=sqrt(pconst.^2-((w.^2./c^2).*n3^2));