System of 2nd order ODE with Euler.
More Answers (1)
- y0 must be a 4x1 vector, not a scalar.
- ye = zeros(4,n) instead of ye=zeros(1,n)
- ye(:,1) = y0 instead of ye(1) = y0
- ye(:,i+1) = ye(:,i) + h*f(x(i),ye(:,i)) instead of the expression in your loop
- dy = [Y(3);Y(4);Y3;Y4] instead of the row vector in your code