Shooting Method on ODE

3 次查看（过去 30 天）

Alexander Kimbley 2019-3-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/448856-shooting-method-on-ode

评论： darova 2019-3-7

MagneticShootingMethod.m

Hi, I'm pretty new to Matlab so bare with me please. I need use the ODE below in the ode45 command for the shooting method but I'm unsure how to write it in the correct form i.e. Y_1=... and Y_2=... if anyone could be so kind to re-write this for me I should be fine.

Thanks.

2 个评论
显示无隐藏无

Jan 2019-3-7

在 MATLAB Online 中打开

The attached file MagneticShootingMethod.m contains this line only:

((x-c)^2 - B^2)*(y''-ay)=(2*B^2)*(y-(x-c)y')/((x-c)^2)

This is neither a clear description nor running Matlab code. I cannot guess, what "i.e. Y_1=... and Y_2=..." could mean.

Please post, what you have tried so far and explain the problem clearly.

darova 2019-3-7

在 MATLAB Online 中打开

function dy = odefcn(x,y,a,B,c)
dy = zeros(2,1);
dy(1) = y(2);
dy(2) = (2*B^2)*(y(1)-(x-c)*dy(1))./((x-c).^2);
dy(2) = dy(2) ./ ((x-c).^2 - B^2) + a*y(1);
end

请先登录，再进行评论。

请先登录，再回答此问题。