21.510.50-0.5-10102030405060708090100
>> opt=odeset; opt.RelTol=1e-6;
>> [t,y]=ode45('lorenzeq',[0,t_final],y0,opt); >> plot(t,y) %刚性解法 >> opt=odeset; opt.RelTol=1e-6;
>> [t,y]=ode45('lorenzeq',[0,t_final],y0,opt); >> plot(t,y)
21.510.50-0.5-10102030405060708090100
8、 试求出习题3中给出的微分方程边值问题数值解， 绘制出y(t)曲线，并和该习题得出
的解析解比较精度。
??u?u??0?22?x?y??9、 试用数值方法求解偏微分方程?ux?0,y?0?1,u??x?0,y?0??22y?0,x?0?0，并绘制出u函数曲面。
>> pdetool