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')