how to plot by varying initial conditions
2 次查看(过去 30 天)
显示 更早的评论
%plot n vs x n vs y
for n=1:10
[t,r]=ode45(@fn,[0 140],[0.01+n 0.01 ])
plot (n,r(:,1))
end
function dr=fn(t,r)
x=r(1); y=r(2);
dr=zeros(2,1)
dr(1)=x+y;
dr(2)=x-y;
end
pl plot n vs x
0 个评论
回答(1 个)
KSSV
2021-7-25
function myfun()
%plot n vs x n vs y
R = cell(10,1) ;
figure
hold on
for n=1:10
[t,r]=ode45(@fn,[0 140],[0.01+n 0.01 ]) ;
R{n} = r(:,1) ;
plot(t,R{n})
end
end
function dr=fn(t,r)
x=r(1); y=r(2);
dr=zeros(2,1) ;
dr(1)=x+y;
dr(2)=x-y;
end
4 个评论
KSSV
2021-7-25

This is the plot obtained. n is a single number, how you can plot? Plot wrt to t and go for legend. Check is code giving correct results.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!