exm041452_1.m

matlab mfiles实例分析

sinit=bvpinit(linspace(0,4,5),[1;0]);    

c=1;
sol=bvp4c(@twoode,@twobc,sinit,[ ],c);						%<3>  

x=linspace(0,4,100);
y = bvpval(sol,x);
plot(x,y(1,:),'b-',sol.x,sol.y(1,:),'ro',sinit.x,sinit.y(1,:),'ks')
legend('\fontname{隶书}\fontsize{20}插值后的解曲线','解点','猜测解点',0)