Info
此问题已关闭。 请重新打开它进行编辑或回答。
Subscript indices must either be real positive integers or logicals.
1 次查看(过去 30 天)
显示 更早的评论
I am trying to run eulers method for a first order differential equation. However recieve the above error when I try to plot the results.
See my code blow:
x=zeros(400,1);
y=zeros(400,1);
x(1)=2;
y(1)=4;
for n=1:400
x(n+1)=x(n)+0.005;
y(n+1)=x(n)+0.005.*(((y).^2)+(7*y)+2)/(3*x((y).^2)+((1-x)*(3*((y)^2)+7)));
end
plot(x,y)
title('eulers')
grid on
What am I doing wrong?
0 个评论
回答(1 个)
KL
2017-12-1
when you write,
y(n+1)=x(n)+0.005.*(((y).^2)+(7*y)+2)/(3*x((y).^2)+((1-x)*(3*((y)^2)+7)));
^
you're trying to access the zeroth element of x (since y is vector with all zeroes except its first element). You need to rewrite that equation. How does the equation in your book look like?
Anyway, if you want to create a vector with 0.005 increaments use linspace or colon(:) operator.
2 个评论
KL
2017-12-1
if you have access to symbolic math toolbox, it's a lot easier. Go through the example here: https://www.mathworks.com/help/symbolic/solve-a-single-differential-equation.html#f1-7555
此问题已关闭。
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!