foz.m

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

function  b=foz(u,t,d1,d2,n1,n2,n3);
dx=0.001;
x0=u-1.e-10;
x1=x0-dx;
f0=slab(x0,t,d1,d2,n1,n2,n3);
f1=slab(x1,t,d1,d2,n1,n2,n3);
if f0*f1<0
N=1;
eps = 1.e-10;
fnew=1;
while (abs(fnew)>=eps) & (N < 200)
  N=N+1;
  %xnew=x1-slab(x1,t,d1,d2,n1,n2,n3)*((x1-x0)/(f1-f0));
  xnew=(x0+x1)/2;
  fnew=slab(xnew,t,d1,d2,n1,n2,n3);  
  if fnew*f1 < 0 
     x0=xnew;
     f0=fnew;
  elseif fnew*f1 > 0
     x1=xnew;
     f1=fnew;
  end
[x0 x1 f0 f1];
end;
b=xnew;
%fnew
%N
if N>= 200
   'not convergent'
end
elseif f0*f1>0, b=0;
end