%bvp4c_example.m
%Mathieu方程在邊界條件下的解
lambda = 15;                                                 %未知參數(shù)猜測解
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)]);