Drexel dragonThe Math ForumDonate to the Math Forum


Math Forum
 Ask Dr. Math
 Discussions
 Internet Newsletter
 MathTools
 Problems of the Week
 Teacher2Teacher
 Teacher Exchange
 Workshops

Search All of the Math Forum:
Browse our Internet Mathematics Library

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


Math Forum » Discussions » Software » comp.soft-sys.matlab

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

Advanced Search
Reply to this Topic Reply to this Topic

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
mojtaBa malekpour

Posts: 1
Registered: 10/16/12
ode45 help
Posted: Oct 16, 2012 5:06 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic 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
Read ode45 help
mojtaBa malekpour
10/16/12
Read Re: ode45 help
Torsten
12/7/12
Read Re: ode45 help
Kristin

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2013. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.