use MATLAB(R) to solve a previously defined system, f, with ODE45 - Maple Help

Online Help

All Products    Maple    MapleSim

Home : Support : Online Help : Mathematics : Numerical Computations : Matlab : Matlab/ode45

Matlab[ode45] - use MATLAB(R) to solve a previously defined system, f, with ODE45

Matlab[ode15s] - use MATLAB(R) to solve a previously defined system, f, with ODE15s

Calling Sequence

ode45(f, Trange, IC)

ode45(f, Trange, IC, tol=Tol)

ode15s(f, Trange, IC)

ode15s(f, Trange, IC, tol=Tol)




string naming a function defined in MATLAB®



range, t_initial..t_final



initial conditions vector; for example, [0, 0, 0]



optional floating-point setting for tolerance (default is .001)



The ode45 command uses MATLAB® to compute the ODE45 solution of a differential system. The ode15s command uses MATLAB® to compute the ODE15S solution of a differential system.


To be valid, the call must name the function (f) defined in MATLAB®, and specify both the time range (Trange) and the initial condition vector (IC).


The function f must either be a built-in MATLAB® function, or a function defined in the file f.m in the active MATLAB® path


The time range must contain both a start time (t_initial) and an end time (t_final).


The initial condition (IC) must have a dimension consistent with the dimension of the system defined in the function, f.


The tolerance option (specified using 'tol'=Tol) must be a floating-point value.  It defaults to .001 if not specified.


Executing the ode45 command returns two Vectors: (T, X). T is the Vector of time steps where f was evaluated, and X is a Vector of values of f, evaluated at times in T.



Call ode45 to work on a built-in MATLAB® function.


Call ode15s to work on a built-in MATLAB® function.


Create a user-defined function, rocket, in Maple.  The file may also be predefined in the current path.

file:=openC:MATLABwork ocket.m,WRITE:

filecontents:= function yp = rocket(t,y) % time the rocket catches the target global INTERCEPT; % minimum distance representing intercept global DMIN; % ratio of pursuer speed to target speed k=1.3; % find speed and position of target % (other equations can be substituted here) if t < 10 % target changes direction after 10 seconds p = [ 10; t; 0 ]; vt = [ 0 ; 1; 0 ]; else p = [ 10; 10; t-10 ]; vt = [ 0 ; 0; 1 ]; end d = sqrt(sum((p-y).^2)); % calculate distance between P and T if d < DMIN % check if pursuer has caught targetN if t < INTERCEPT % if this is the first time, set the INTERCEPT = t; % interception time end k = 1; % slow down the pursuer end vp = k*sqrt(sum(vt.^2)); % set speed of pursuer yp = vp*(p-y)/d; % determine new position of pursuer &colon;



Use the ssystem command to set the permissions on UNIX systems.

ssystemchmod a+r rocket.m&colon;

Open the "Matlab Link".


Before accessing the user-defined function, the file needs to be opened (in MATLAB®).

evalMopen('C:MATLABwork ocket.m')&colon;

Set global variables (in MATLAB®).



Setup the problem.




Solve the problem.


Plot the results.




See Also

dsolve, Matlab, Matlab[evalM], MatlabMatrix

Download Help Document

Was this information helpful?

Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam