It looks like you're new here. If you want to get involved, click one of these buttons!

- 141K All Categories
- 103.8K Programming Languages
- 6.5K Assembler Developer
- 1.9K Basic
- 40K C and C++
- 2.9K C#
- 7.9K Delphi and Kylix
- 4 Haskell
- 9.7K Java
- 4.1K Pascal
- 1.3K Perl
- 2K PHP
- 551 Python
- 37 Ruby
- 4.4K VB.NET
- 1.6K VBA
- 20.9K Visual Basic
- 2.6K Game programming
- 317 Console programming
- 92 DirectX Game dev
- 1 Minecraft
- 112 Newbie Game Programmers
- 2 Oculus Rift
- 9K Applications
- 1.8K Computer Graphics
- 744 Computer Hardware
- 3.5K Database & SQL
- 535 Electronics development
- 1.6K Matlab
- 628 Sound & Music
- 258 XML Development
- 3.3K Classifieds
- 199 Co-operative Projects
- 198 For sale
- 190 FreeLance Software City
- 1.9K Jobs Available
- 604 Jobs Wanted
- 210 Wanted
- 2.9K Microsoft .NET
- 1.8K ASP.NET
- 1.1K .NET General
- 3.4K Miscellaneous
- 7 Join the Team
- 355 Comments on this site
- 70 Computer Emulators
- 2.1K General programming
- 187 New programming languages
- 629 Off topic board
- 217 Mobile & Wireless
- 89 Android
- 126 Palm Pilot
- 340 Multimedia
- 156 Demo programming
- 184 MP3 programming
- Bash scripts
- 27 Cloud Computing
- 53 FreeBSD
- 1.7K LINUX programming
- 370 MS-DOS
- Shell scripting
- 321 Windows CE & Pocket PC
- 4.1K Windows programming
- 942 Software Development
- 417 Algorithms
- 68 Object Orientation
- 92 Project Management
- 95 Quality & Testing
- 269 Security
- 7.7K WEB-Development
- 1.8K Active Server Pages
- 62 AJAX
- 5 Bootstrap Themes
- 55 CGI Development
- 28 ColdFusion
- 224 Flash development
- 1.4K HTML & WEB-Design
- 1.4K Internet Development
- 2.2K JavaScript
- 37 JQuery
- 308 WEB Servers
- 151 WEB-Services / SOAP

airwalkery2k
Member Posts: **2**

in Matlab

Hi,

I am using Matlab for a class, and unfortunately, I don't have much experience in it.

For my problem, I need to solve an ordinary differential equation using Matlab's built-in ODE45 function. Using some tutorials, I came up with the following code to solve and plot my function.

[b]function jtd

[t,x] = ode45(@dfile,[0,20],[0;0]);

plot(t,x(:,1))

title('nonlinear')

xlabel('t'), ylabel('y'), grid

function xprime = dfile(t,x)

F=13.4;

a=0.1;

xprime = zeros(2,1);

xprime(1) = x(2);

xprime(2) = F*cos(t) - a*x(2) - x(1)^3;[/b]

It works well, but I need to be able to vary F and a using a loop to get multiple plots as the two variables change. Yet I can't for the life of me figure out how to get a variable from the main function JTD into the function xprime when I call the ode45 function.

If I could figure that out, it would help me immensely in continueing my program.

I am using Matlab for a class, and unfortunately, I don't have much experience in it.

For my problem, I need to solve an ordinary differential equation using Matlab's built-in ODE45 function. Using some tutorials, I came up with the following code to solve and plot my function.

[b]function jtd

[t,x] = ode45(@dfile,[0,20],[0;0]);

plot(t,x(:,1))

title('nonlinear')

xlabel('t'), ylabel('y'), grid

function xprime = dfile(t,x)

F=13.4;

a=0.1;

xprime = zeros(2,1);

xprime(1) = x(2);

xprime(2) = F*cos(t) - a*x(2) - x(1)^3;[/b]

It works well, but I need to be able to vary F and a using a loop to get multiple plots as the two variables change. Yet I can't for the life of me figure out how to get a variable from the main function JTD into the function xprime when I call the ode45 function.

If I could figure that out, it would help me immensely in continueing my program.

Terms of use / Privacy statement / Publisher: Lars Hagelin

Programmers Heaven articles / Programmers Heaven files / Programmers Heaven uploaded content / Programmers Heaven C Sharp ebook / Operated by CommunityHeaven

© 1997-2017 Programmersheaven.com - All rights reserved.

## Comments

217[b]F=13.4;

a=0.1;

fun=@(x,t) dfile(x,t,F,a);[/b]

[t,x] = ode45([b]fun[/b],[0,20],[0;0]);

plot(t,x(:,1))

title('nonlinear')

xlabel('t'), ylabel('y'), grid

function xprime = dfile(t,x,[b]F,a[/b])

xprime = zeros(2,1);

xprime(1) = x(2);

xprime(2) = F*cos(t) - a*x(2) - x(1)^3;

I define the function handle of dfile out of the ode45 function. In this way you can pass the parameter to the function.

Now you can make a for cicle to make different plots, for example:

F=[12.3 13.4];

a=[0.2 0.1];

for i=1:length(F)

figure();

fun=@(x,t) dfile(x,t,F(i),a(i));

...

end

I add the function figure() so that plots on different figures.

2From here on out, it should be smooth sailing. [b]Thanks, Giug![/b]

2171