normal.m

ok &ouml nemli bir kaybnkao lasdhlas sdohfs naosfnsdn sdondsfgpdsmr ljrg fdbh

function  NNJh=normal(x,t,d1,d2,n1,n2,n3)
B=x;

k1h=t*2*pi*sqrt(n1*n1-B*B);
k2h=t*2*pi*sqrt(B*B-n2*n2);
k3h=t*2*pi*sqrt(B*B-n3*n3);
k2D2=d2*2*pi*sqrt(B*B-n2*n2);
k3D1=d1*2*pi*sqrt(B*B-n3*n3);
k1D2=d2*2*pi*sqrt(n1*n1-B*B);

%if k2D2 == 0
%    'Coth(k2D2)= Inf (Sonsuz) Hatasi ';
%end

N1Jh=k1h*k1h+k2h*k2h*(coth(k2D2))*(coth(k2D2))+(k1h*k1h-k2h*k2h*(coth(k2D2))*(coth(k2D2)))*sin(4*k1h)/(4*k1h)+k2h*coth(k2D2)*(1-cos(4*k1h))/2;
N2Jh=k1h*k1h*coth(k2D2)/(2*k2h)-k1h*k1D2/(2*(sinh(k2D2))^2);
N3Jh=(k1h*cos(2*k1h)+k2h*sin(2*k1h)*coth(k2D2))^2*(sinh(2*k3D1)-2*k3D1)/(4*k3h*(sinh(k3D1))^2);
NNJh1=N1Jh+N2Jh+N3Jh;
NNJh=real(NNJh1);
   
%  t=h/lambda
%  d1=D1/lambda
%  d2=D2/lambda