2012-06-18 55 views
1

這將是它: Y'(T)= Y(t)的/ TT^2/Y^2 * T Y(1)= 1用matlab解這個微分方程?

我曾嘗試: 功能hazi3b()

[T,Y] = ode45(@bfugveny,[112],1); 情節(T,Y, 'GX')

和:

功能DY = bfugveny(T,Y)

DY = Y(T)/噸 - 噸^ 2 /(Y^2 * T);

+0

出了什麼問題? –

+0

試圖訪問y(1.22);索引必須是正整數或邏輯。 bfugveny中的錯誤(第3行) dy = y(t)/ t-t^2 /(y^2 * t); (1),y + f * hB(:1),odeArgs {:});在ode45中的錯誤(行262) f(:,2)= feval(odeFcn,t + hA hazi3b(第2行)中的錯誤 [T,Y] = ode45(@bfugveny,[112],1); –

回答

2

您不必在公式中寫出y(t)

傳入你的oracle的y已經是y -evaluated-at-time-t的猜測。

所以儘量

dy = y/t - t^2/(y^2*t); 
+0

謝謝,它工作:) –

+0

所以接受答案....感恩是好的,但接受對你和幫助你的人是好的。也投了票。 – duffymo