exm041453_1.m

来自「matlab mfiles实例分析」· M 代码 · 共 18 行

clear
x00=linspace(0,pi,4);							%均布初始网点
y00=inline('[cos(4*x);-4*sin(4*x)]');		%解的初始猜测形式	<3>
lmb=15;											% 的初始猜测
s0=bvpinit(x00,y00,lmb);						%得到初始猜测网  

opts=bvpset('AbsTol',0.5,'RelTol',0.38,'Stats','on');%选项设置<6>
s1= bvp4c(@mat4ode,@mat4bc,s0,opts);			%获得第一近似解
Lambda1=s1.parameters  

opts.AbsTol=1e-6;opts.RelTol=1e-3;					%重置计算精度
s2= bvp4c(@mat4ode,@mat4bc,s1,opts);			%获得第二近似解
Lambda2=s2.parameters  

plot(s0.x,s0.y(1,:),'ks--',s1.x,s1.y(1,:),'bo:',s2.x,s2.y(1,:),'r*-')
legend('\fontname{隶书}\fontsize{16}猜测解','第一近似解','第二近似解',0)  
axis([0,pi,-1,1.2]),xlabel('x'),ylabel('solution y')