bvp4c--TWOBVP

1. y''+ abs(y)=0; y(0)=0 y(4)=-2

令 y(1)= y; y(2)=y'则微分方程可以写成标准形式

function dydx=ode(x,y)

dydx=[y(2)

          -abs( y(1) )]; % 列向量

边界条件bc( y(a),y(b))=0

function res=bcfun(ya,yb)
res=[ya(1)
     yb(1)+2];% ya(1)表示自变量为a时候的第一个因变量y(1),ya(2)表示自变量为a时候的第二个因变量y(2)

初始假设solinit

solinit=bvpinit(linspace（0,4,10）,[0,0]);

解方程：

sol=bvp4c(@ode,@bcfun,solinit)

源程序

function example
solinit=bvpinit(linspace(0,4,10),[0,0])
sol=bvp4c(@ode,@bcfun,solinit)
plot(sol.x(:),sol.y(1,:))
%%%微分方程
function dydx=ode(x,y)

dydx=[y(2)
      -abs( y(1) )]; % 列向量
 
%%%边界条件
function res=bcfun(ya,yb)
res=[ya(1)
     yb(1)+2];% ya(1)表示自变量为a时候的第一个因变量y(1),ya(2)表示自变量为a时候的第二个因变量y(2)

matlab 中源文件见：

MATLAB6p5\toolbox\matlab\demos\twobvp.m