findroot.m

solvin Parbolic forms of PDE explicitly. The error is O( h2+k2) or O(2,2) , it doesn鈥檛 develop the

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