?? exm041453_1.m
字號:
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')
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -