Items tagged with ode


I am solving a large system of ODEs, using the following command,

> Sol := dsolve({seq(ode[j], j = 0 .. 21), seq(v[j](0) = 0, j = 1 .. 21), v[0](0) = 1}, [seq(v[j](t), j = 0 .. 21)]);



and then plot the quantities I want by something like

> plots:-odeplot(Sol, [t, v[3](t)+v[5](t)], t = 0 .. 1.5);

My problem is that, I do not know a priori which quantity I want to plot, and plotting using above method requires solving the ODEs each time separately, which takes a long time.


So I was curious if there is a scheme that I can solve my system for once and for all, and then plot any quantities that I would like to see.


The idea: parameter "a" will have a new random value each 10 days.

The way I did it is working but it can get very long especially if I do it for a system of equations and for long time more than a year (365 days).

The code:

with(DEtools); with(plots);
n := 5;

for i to n do Ra[i] := RandomTools:-Generate(distribution(Uniform(0.1e-1, .5))); a[[i]] := Ra[i] end do;

b := 0.1e-2;

T := 10;

 eq := diff(L(t), t) = a*L(t)-b;

init[1] := L(0) = 100;
 sol[1] := dsolve({init[1], subs(a = a[[1]], eq)}, L(t), range = 0 .. T, numeric);

init[2] := L(T) = rhs(sol[1](T)[2]);

sol[2] := dsolve({init[2], subs(a = a[[2]], eq)}, L(t), range = T .. 2*T, numeric);


init[3] := L(2*T) = rhs(sol[2](2*T)[2]);
sol[3] := dsolve({init[3], subs(a = a[[3]], eq)}, L(t), range = 2*T .. 3*T, numeric);

p[1] := odeplot(sol[1], [t, L(t)], t = 0 .. T);

p[2] := odeplot(sol[2], [t, L(t)], t = T .. 2*T);

p[3] := odeplot(sol[3], [t, L(t)], t = 2*T .. 3*T);

p := display([p[1], p[2], p[3]]);


Thank you

Hi Guys,

I am trying to solve the folling ODE in maple but am struggling to get to the correct solution. I first of all have the following equations;

Diff(V(x), [x]) = q(x)

Diff(M(x), [x]) = V(x)

theta(x) = Diff(v(x),[x])


Which give me the following;

q(x) = Diff(EI*(Diff(v(x), [`$`(x, 2)])), [`$`(x, 2)])

With thw boundary conditions;

M(l) = 0, M(0) = 0, v(0) = 0, v(l) = 0

The given solution is;


Anything that might point me in the right direction would be great!





Attached is a photo with the code I am working for.  

On the top is practice code with a simpler ODE to help with trouble shooting, on the bottom is the ODE I am working with.

I was hoping to gain insight about the _z1 symbol in the solution, I haven't been able to find much help on other threads.  I would like to know how I can go about working with it - if it is something on my end or if it is the nature of the equation I am working with.


Thank you for any help,


This is the simplest method to explain numerically solving an ODE, more precisely, an IVP.

Using the method, to get a fell for numerics as well as for the nature of IVPs, solve the IVP numerically with a PC, 10steps.

Graph the computed values and the solution curve on the same coordinate axes.


y'=(y-x)^2, y(0)=0 , h=0.1

Sol. y=x-tanh(x)


I don't know well maple. 

I study Advanced Engineering Math and using maple, but i am stopped in this test.

I want to know how solve this problem.

please teach me~ 

IT IS EULER's method

Hi there. I'm Student

i want to know how solve this problem.

please teach me! 

y'=(y-x)^2, y(0)=0, h=0.1


how solve this problem for maple? 

please teach me~

Hi, I am an student and I am currently working on a system that sketches the relation of predator and prey of yellowstone's gray wolf and elk. I tried using the Lotka-Volterra model, but I wanted to add more parameters and add a carrying capacity for the system. Unfortunatley I cannot find a way to edit the Lotka model to my needs, and because I am new i do not know how to create my own model. This is the two equations I want to use: (D(x))(t) = alpha*x(t)-ax^2/k-b*x(t)*y(t)-gx(t), (D(y))(t) = -beta*y(t)+c*x(t)*y(t)-gy(t)

were k is carrying capacity.

Basically what I am asking is that if someone can help make the system workable on Maple and some steps of how to do it. 

Hello guys,

I am stuck with the following code. I cannot get a solution. Mapple doesn't even run. Please help, I am in urgent need on help!!!



I've got the following four differential equations :

d2v_x:=-((C_d)*rho*Pi*(r^2)*(v_x)*sqrt((v_x)^2 +(v_y)^2))/(2*m);
d2v_y:=-((C_d)*rho*Pi*(r^2)*(v_y)*sqrt((v_x)^2 +(v_y)^2))/(2*m)-g;

and the following initial value conditions:

x(0)=0,y(0)=0,v_x(0)=v0/sqrt(2),v_y(0)=v0/sqrt(2) given v0=65 

I need to solve these using the numeric type and then draw overlaid plots

(i) setting C_d=0

(ii) leaving C_d as a variable

before plotting y(t) vs x(t). The hint for this last part is that the path can be seeing using [x(t),y(t)] instead of [t,y(t)]

I've tried to do it but seemed to have several syntax errors.



I've got the following diff.eq

y'(x)=sin(x*y(x)) given y(0)=1 

and need to solve it numerically which is why I've used:


This code doesn't return a value though and in fact, ans3 is being displayed as a procedure

"ans3:=proc(x_rkf45) ... end proc"

I don't quite understand why and what I need to do to get the required numerical solution


In the attached Maple worksheet I attempt to plot the solution of an initial value problem for a first order ODE.  DEplot fails with a cryptic message.  Strangely enough, if I give the "arrows=none" option to DEplot, it produces the correct plot!

I see this behavior in Maple 17 and 18.

Maple 11, however, works fine with or without the "arrows=none" option.

Is there an explanation for this or is it a bug?

hai everyone. i am currently trying to solve an integration of the following ∫g(η)dη . integrate from 0 to 10.

from the following odes.

f ''' +1-(f ')2 +ff ''=0,


with boundary conditions f(0)=0, f'(0)=λ, f'(∞)=1, g(0)=1,g(∞)=0

First, i solve the odes using the shooting method. then i used the trapezoidal rule to solve for the integration of g(eta) using the following codes

> with(student);
> trapezoid(g(eta), eta = 0 .. 10, 10);
> evalf(%);

it seems that it can not read the data from the shooting method. can anyone suggest why it is happening?

thank you verymuch for your concern :)

So I am trying to solve a given ODE using calculated christoffel symbols found by maple, and in order to get the correct christoffel symbols, I need my function to be r=x-x_s, where x is not a function of t. However, I then have to solve the ODE where x is a function of t. Maple used the r value to find the christoffel symbols which has x in it, and now I want to find the origonal function of x(t), but I can't have x and x(t) in the same ODE. If I change r=x-x_s to r=x(t)-x_s, I get the wrong christoffel symbols. How can I solve my ODE?

I try to solve a set of differential equations with boundary conditions by dsolve command but I got this error:


odes := D2*(diff(P(x), x))/((3*D1*a+4*D2)*P(x))-(diff(S(x), x))/(q*S(x)-1) = 0, diff(S(x), `$`(x, 2))+(diff(S(x), x))*cotx+4*pi*(3*D1*a+4*D1)*P(x)/((q*S(x)-1)*D2) = 0

ics := P((1/2)*pi) = 1, S((1/2)*pi) = -1, (D(S))((1/2)*pi) = 0, (D(P))((1/2)*pi) = 0

sol := dsolve({ics, odes}, numeric);
Error, (in dsolve/numeric/process_input) input system must be an ODE system, got independent variables {x, (1/2)*pi}


Does any body knows what the problem is??



i want to solve the 2nd order ODE y''=0.8y'arctan(x)+xy    -2<x<1

the maple code used is:

>maplesol := dsolve({ (D(D(y)))(x) =0.8*arctan(x)*(D(y))(x)+x*y(x),y(-2)=1.3,D(y)(-2)=0.7});

and maple only gives out maplesol:=


I think its something to do with the arctan(x) because when i remove it, i get an answer.

First 9 10 11 12 13 14 15 Last Page 11 of 31