[問題] 二階微分方程畫圖的問題.已回收
小弟是使用function的初學者.
我的老師出了一個問題是
y'' - y = sinx
initial condition是
y(0) - y'(0) = 1
要我畫出所有(y(pi),y'(pi))的點
我在用ode45 function時.tspan是用[0,pi]
但是我要怎麼畫才能只畫出 (y(pi),y'(pi))
我的function.m file是
function w = rhs5(x,y)
w = [y(2); y(1)+sin(x)];
end
程式內容是
y = [0;0];
for j = -5:1:5
y(1) = j;
y(2) = y(1)-1;
[T,Y] = ode45(@rhs5,[0,pi],[y(1),y(2)]);
plot(Y(:,1),Y(:,2));
hold on
axis([-10 10 -10 10]);
end
我知道我這樣做能化出很多的(y(x),y'(x)),
但是接下來就不知道怎麼做下去了
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 173.70.6.63