bvp4c_fun.m

matlab从入门到精通第6章书后源码对初学者很有帮助

% bvp4c_fun.m
% 求解微分方程边值问题
S=[ 0 0
    0 -2];
options=bvpset('SingularTerm',S);
% 得到微分方程的初场
guess=[sqrt(3)/2;0];
solinit=bvpinit(linspace(0,1,5),guess);
sol=bvp4c(@emdenode,@emdenbc,solinit,options);
% 根据分析解计算
x=linspace(0,1);
truy=1./sqrt(1+(x.^2)/3);
% 绘制图形结果
plot(x,truy,'b-');
hold on;
plot(sol.x,sol.y(1,:),'ro');
legend('analytical','computed');
xlabel('x'),ylabel('y');
hold off;