参考代码： % 解的初始估计solinit = bvpinit(linspace(0,1,100),[1 -50]); % BVP问题求解ode = @(t,x) [x(2); -pi^2/4*(1-x(1))];bc = @(ya,yb)[ya(2)-pi; yb(2)+pi/2];sol = bvp4c(ode,bc,solinit); % 结果绘图t=sol.x;vars={'x', 'x'''};for ...
