bvp4c_example.m

matlab7.0书中的实例程序的源代码

%bvp4c_example.m
%Mathieu方程在边界条件下的解
lambda = 15;                                                 %未知参数猜测解
solinit = bvpinit(linspace(0,pi,10),@Mathieu_initfun,lambda);%求初始值
sol = bvp4c(@bvp_Mathieufun,@Mathieu_bcfun,solinit);         %求解方程
xint = linspace(0,pi);                                       %画图
Sxint = deval(sol,xint);
plot(xint,Sxint(1,:))
axis([0 pi -1 1.1])
title('Mathieu方程在边界条件下的解'); 
xlabel('x')
ylabel('y')
text(0.1,1,['4阶Mathieu方程的特征值 = ' num2str(sol.parameters)]);