d^2y/dt^2=(1/126)*t-1/18

이런 2계 미분 방정식이 있는데

dsolve로 풀고 ode45로 풀어야하는데

dsolve는 했습니다...

ode45로 코드를 짜봤는데

이상하게 나오네요..

function ydot=rhss(t,y)
ydot(1)=y(2)
ydot(2)=(1/126)*t-1/18
ydot=[ydot(1);ydot(2)];
end


명령창

 [t,y]=ode45('rhss',[0,7],[0,0])

 plot(t,y)


알려주세요제발....내일까진데 모르겠어요