Items tagged with phase



I have the ODEs: x'=x



sketch the phase space and extended phase space of previous ode.


hello guys,


i have a system of autonomous equations which i want to plot its 3D phase space with directional field,

i have some problem with it , and i dont know how to command for add some directional field for 3D phase space .


thank you guys


Hi there,
I have a set of differential equations whose solution, Jacobian matrix and its eigenvalues, direction field, phase portrait and nullclines, need to be computed.

Each of the equations has a varying parameter.

I know how to get the above for a single parameter value, but when I set a range of values for the parameters, Maple is not able to handle all cases as I would expect: solving the differential equation system:

eq1 := x*(1.6*(1-(1/100)*x)-phi*y)
eq2 := (x/(15+x)-0.3e-1*x-.4)*y+.6+theta
desys := [eq1, eq2];
vars := [x, y];
steadyStates := map2(eval, vars, [solve(desys)])

already yields an error:
Error, (in unknown) invalid input: Utilities:-SetEquations expects its 2nd argument, equations, to be of type set({boolean, algebraic, relation}), but received {-600*y+(Array(1..2, {(1) = 8400, (2) = 15900})), Array(1..5, {(1) = 0, (2) = 0, (3) = 0, (4) = 0, (5) = 0})}

The equations are the following:
de1 := diff(x(t), t) = x(t)*(1.6*(1-(1/100)*x(t))-phi*y(t));
de2 := diff(y(t), t) = (x(t)/(15+x(t))-0.3e-1*x(t)-.4)*y(t)+.6+theta

the parameters being:
phi:=[0 0.5 1 1.5 2]
theta:=[5. 10.]

How can I handle the situation so that Maple computes each of the above for each combination of the parameters?

I would like to avoid using two for loops and having to store all results in increasingly bigger and complicated arrays.

The worksheet at issue is this:


Hi there,

I would like to compute and display the nullclines of a set of ordinary differential equations.

AFAIK, I can compute the nullclines in Maple by defining the equations and solving the system


# Define the equations
eq1 := u(t)*(1-u(t)/kappa)-u(t)*v(t) = 0;
eq2 := g*(u(t)-1)*v(t) = 0;

# Solve the system (i.e. compute the nullclines)
sol := solve({eq1, eq2}, {u(t), v(t)});

However, I am not quite able to imagine how to display them over a dfieldplot or a phaseportrait.

Attached is an example with some differential equations, and their vector field and trajectories:

It can be use to illustrate how to (compute and) display the nullclines.


Thank you,


Hi all,

I have a simple transfer function like



I want to calculate the amplitude ratio and the phase shift.

Without maple I would set s=I*omega and solve for real and imaginary part to obtain the phase shift and amplitude ratio.



Is there any possibility to do that in Maple more easily?



> with(DETools);
> eq1 := diff(y(t), t) = v(t);
--- y(t) = v(t)
> eq2 := diff(v(t), t) = -3*v(t)+10*y(t);
--- v(t) = -3 v(t) + 10 y(t)
> phaseportrait([eq1, eq2], [y, v], t = 0 .. 10, [[y(0) = 0, v(0) = 1]], y = -10 .. 10, v = -10 .. 10, linecolor = blue);



why my phase portrait not working

Below is my code for solving a solution of three first order ODEs using the 4th-order Runge-Kutta method.

I have been able to successfully plot the solutions of each of the ODEs (x,y,z) against time t, however I am struggling to produce a plot in the three dimensional phase domain of x,y,z. Could anybody suggest what commands to use as everything I have tried (plot, plot3d, implicitplot3d etc) has produced an error. 

h:= 0.01:
N:= 200:

Hey, I have a system of ODE and I can't draw it's phase curve. I tried to use DEplot and phaseportrait, but it doesn't work. Here is my system:


dy/dt=ky              (k is a constant)


Here is my piece of code:

DE := [diff(x(t), t) = x(t)];

DF := [diff(y(t), t) = k*y(t)];


phaseportrait([DE, DF], [y, x], t = -5 .. 5, y = -5 .. 5, x = -5 .. 5, k...

Hey, I'm trying to plot the phase curve of the simple system of ODE, BUT it shows me an error which one I don't understand at all :)

Here is my code:

DE := diff(x(t), t) = x(t);
DF := diff(y(t), t) = k*y(t);
phaseportrait([DE, D], [y, x], t = -5 .. 5, [[y(0) = 1, x(0)], [y(0) = 0, x(0) = 2], [y(0) = 0, z(0) = -2]], y = -5 .. 5, x = -5 .. 5, color = black, linecolor = red);

Error, (in DEtools/phaseportrait) invalid initial...


I am trying to do a phase plot of an autonomous system using DEplot command. However, no phase plot appears. Could you tell me what am I doing wrong? The same code worked for my professor in class, but it's not working for me. I am using Maple 16. The code is posted below.


> with DEtools

> sys := diff(x(t), t) = y(t), diff(y(t), t) = x(t)*(1-x(t)*x(t))+y(t);

> DEplot([sys], [x(t), y(t)], t = 0 .. 0.1e-8, x = -3 .. 3, y = -3 .. 3, color = black)


I want  to make the phase plot of d(Z(t))/dt (along y) and Z(t) (along x axis). The initial conditiona are arbitrary. I only have one equation:

d(dz/dt)/dz = - (sqrt(d(z^2)/dt + 10^10 * z^2)*dz/dt + 10^(10)*z)/dz/dt

Is it still possible to make a phase space plot or do i need another equation?


I have the following problem:

My function is defined by the determinant of 2 Heun functions

If I plot the phase I get something which looks quite what I'm looking for.

To get a better result I thought I would manually carry out the Wronskian as far as possible...

Doing some manipulations I get another form of the Wronskian which in fact should give the same result...

the problem is it doesnt :-(

I've added the spreadsheet....

Hi everybody,

Does somebody has an idea about how to model the temperature of water during freezing ?

I mean, if you consider a certain mass of water and you put it, for example on a fridge. I don't consider the convection but only the diffusion, and I would like to model, if possible with an animation the temperature of the water (as a cube for the shape).

Normally, during the freezing I should observe a decrease of the temperature, then a freezing...

Hi everyone,


I have this differential equations:


2 2
l := t -> sqrt((x1(t) - x2(t)) + (y1(t) - y2(t)) )


I want to draw the phase portraits (x1(t...



I'm trying to plot a few phase portraits for this four equation system but maple wont let me...

> restart;

> s := 10; d[T] := 0.2e-1; k1 := 2.4*10^(-5); k2 := 3*10^(-3);
> N := 5; delt := .24; c := 2.4; T0 := 1; 


Can someone please help me? Thanks
1 2 Page 1 of 2