Integrating a 2nd order ODE
조회 수: 3 (최근 30일)
이전 댓글 표시
I am given an equation,
d^2y/dx^2 + q(x) = 0
x ranges from 0 to 1 and y(0) = 1 and y(L) = 1.5
where L = 1
and q(x) = 2*cos((pi*x)/L)
here is what I have
function dydx = yfunc (x,y)
x = linspace(0,1,30);
L = 1;
for i = 1:length(x)
qx(i)= 2*cos((pi*x(i))/L);
end
dydx= -qx;
end
%main script
[x,y] = ode45(@yfunc,x,[1 1.5]);
I am running into an error stating Dimensions of arrays being concatenated are not consistent. Could you please help me?
댓글 수: 0
채택된 답변
Star Strider
2021년 10월 17일
syms y(x) x L Y
q(x) = 2*cos(pi*x/L);
Dy = diff(y);
D2y = diff(Dy);
ODE = D2y + q
[VF,Subs] = odeToVectorField(ODE)
bvpfcn = matlabFunction(VF, 'Vars',{x,Y,L})
I solved it completely, however I do not want to deprive you of the same feeling of accomplishment, so I leave the rest to you. It is a straightforward solution. Follow the examples in the documentation I linked to.
.
댓글 수: 0
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Ordinary Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!