findroot.m

来自「solvin Parbolic forms of PDE explicitly.」· M 代码 · 共 17 行

function falpha=findroot(a)

beta=2;
%for i=1:99
%for k=1:20   

[T,x]=findtemp(a);

falpha=beta+((T(2)-T(1))/0.01);

%y=fzero(beta+((T(1,1)-T(1))./(x(i+1)-x(i))),a)
%if falpha<=0.005
  %  disp(a), break, end
%a=a+0.01;
%end
%end