Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Topic: ode45 help
Replies: 2   Last Post: Dec 7, 2012 9:13 AM

 Search Thread: Advanced Search

 Messages: [ Previous | Next ]
 mojtaBa malekpour Posts: 1 Registered: 10/16/12
ode45 help
Posted: Oct 16, 2012 5:06 AM
 Plain Text Reply

why matlab can't solve this?

function dy = ffffd(t,y)
dy = zeros(4,1); % a column vector
dy(1) =y(2);
dy(2) =-81276055.762*cos(atan((0.05*sin(523.5988*t)+0.12*sin(y(3)))/(-y(1)+0.05*cos(523.5988*t)+0.12*cos(y(3)))))*([[y(1)-0.05*cos(523.5988*t)-0.12*cos(y(3))]^2+[0.05*sin(523.5988*t)+0.12*sin(y(3))]^2]^0.5-0.0005)^1.5;
dy(3) =y(4);
dy(4) =171687.934*sin(y(3)-523.5988*t)+1357293887.572*sin(y(3)-atan((0.05*sin(523.5988*t)+0.12*sin(y(3)))/(-y(1)+0.05*cos(523.5988*t)+0.12*cos(y(3)))))*([[y(1)-0.05*cos(523.5988*t)-0.12*cos(y(3))]^2+[0.05*sin(523.5988*t)+0.12*sin(y(3))]^2]^0.5-0.0005)^1.5+122.868*cos(y(3));

[T,Y] = ode45(@ffffd,[0 2.784],[0.1705 0 0 0]);

Date Subject Author
10/16/12 mojtaBa malekpour
10/16/12 Torsten
12/7/12 Kristin

© The Math Forum at NCTM 1994-2016. All Rights Reserved.