MATLAB Coding for BVP

1 次查看（过去 30 天）

显示更早的评论

Joy Salomi 2019-6-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/469118-matlab-coding-for-bvp

评论： Joy Salomi 2020-3-4

function res = mat4bc(ya,yb)

res = [ ya(1)-1

yb(2)

ya(3)

yb(4) ];

The above is for 2 boundary conditions. I need the coding for 3 conditions.

function sol=ex4

ex4init=bvpinit(linspace(0,100),[0 1 0 -1]);

sol=bvp4c(@ex4ode, @ex4bc, ex4init)

end

function f=ex4ode(x,y)

f=[y(2)

-x*x*y(2)

y(4)

-x*x*y(4)

];

end

function res=ex4bc(ya,yb,yc)

res=[ya(1)-3.462*10^(-22)

ya(3)-2.017

yb(1)-1

yb(3)

yc(2)+yc(4)

];

end

This is my program. when I run this program, I'm getting an error. Can someone tell me what is the error here?

6 个评论
显示 4更早的评论隐藏 4更早的评论

Torsten 2019-6-27

编辑：Torsten 2019-6-27

Maybe of interest for you:

https://www.wolframalpha.com/input/?i=y''(t)%2Bt%5E2*y'(t)%3D0

This code gives you the solution for your system:

function main
  fun = @(t,c1,c2)c1*gamma(1/3)*gammainc(t.^3/3,1/3,'upper')/3^(2/3)+c2;
  c2u = 1;
  c1u = (3.462e-22-1)*3^(2/3)/gamma(1/3);
  c2v = 0;
  c1v = 2.017*3^(2/3)/gamma(1/3);
  t = linspace(0,20,200);
  u = fun(t,c1u,c2u);
  v = fun(t,c1v,c2v)
  plot(t,u,t,v)
end