Items tagged with ode ode Tagged Items Feed

i am solving 3 ODE question with boundary condition. when i running the programm i got this error.. any one could help me please.. :)


restart; with(plots); k := .1; E := 1.0; Pr := 7.0; Ec := 1.0; p := 2.0; blt := 11.5

Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))+Gr*theta(eta)-k*(diff(f(eta), eta))+2*E*g(eta) = 0;

diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+Gr*theta(eta)-.1*(diff(f(eta), eta))+2.0*g(eta) = 0


Eq2 := diff(g(eta), eta, eta)+f(eta)*(diff(g(eta), eta))-k*g(eta)-2*E*(diff(f(eta), eta)) = 0;

diff(diff(g(eta), eta), eta)+f(eta)*(diff(g(eta), eta))-.1*g(eta)-2.0*(diff(f(eta), eta)) = 0


Eq3 := diff(theta(eta), eta, eta)+Pr*(diff(theta(eta), eta))*f(eta)+Pr*Ec*((diff(f(eta), eta, eta))^2+(diff(g(eta), eta))^2) = 0;

diff(diff(theta(eta), eta), eta)+7.0*(diff(theta(eta), eta))*f(eta)+7.00*(diff(diff(f(eta), eta), eta))^2+7.00*(diff(g(eta), eta))^2 = 0


bcs1 := f(0) = p, (D(f))(0) = 1, g(0) = 0, theta(0) = 1, theta(blt) = 0, (D(f))(blt) = 0, g(blt) = 0;

f(0) = 2.0, (D(f))(0) = 1, g(0) = 0, theta(0) = 1, theta(11.5) = 0, (D(f))(11.5) = 0, g(11.5) = 0


L := [10, 11, 12];

[10, 11, 12]


for k to 3 do R := dsolve(eval({Eq1, Eq2, Eq3, bcs1}, Gr = L[k]), [f(eta), g(eta), theta(eta)], numeric, output = listprocedure); Y || k := rhs(R[3]) end do

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging





plot([Y || (1 .. 3)], 0 .. 10, labels = [eta, (D(f))(eta)]);

Warning, unable to evaluate the functions to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct





Download tyera(a).mw

Dear all;

Thank you very much is you can make the computation using maple.

I have this system:



I would like to make to change the system in another system using polar coordinate: like x=r*cos(theta) and y=r*sin(theta), and then derive a simple system: r'=F(r, theta) and theta'=G(r, theta), where here F and G are two unknowns functions. 

Many thinks if someone can help me using a Maple code. 


Hi every body:

i have a second order ode and will convert to two ode of first order with maple,how do this work???

eq := diff(y(x), x, x)+2*y(x)+y(x)^2 = 0

Let N be an integer. 


For each pair of integers (n,m) where 1<= n,m <= N, we have a variable f_{n,m}(t). 


Then for these we have a system of ODEs 


d/dt f_{n,m}(t) = \sum_{n', m'} f_m'n' * f_m''n'' * (m'n'' - m''n') 


where m''=m-m', n''=n=n', and the sum is simply over for all pairs (n',m'). 


I simply do not know how to put these set of equations into Maple in a nice way. 


I will really appreciate any help!

I'm trying to solve this system of ODEs by Laplace transform. 

> de1 := d^2*y(t)/dt^2 = y(t)+3*x(t)

> de2 := d^2*x(t)/dt^2 = 4*y(t)-4*exp(t)

with initial conditions 

> ICs := y(0) = 2, (D(y))(0) = 3, x(0) = 1, (D(x))(0) = 2



> deqns := de1, de2


> var := y(t), x(t)


I need to solve it for both y(t) and x(t), I have tried this by:

> dsolve({ICs, deqns}, var, method = laplace)


> dsolve({ICs, deqns}, y(t), method = laplace)

> dsolve({ICs, deqns}, x(t), method = laplace)


However I get this error message:

Error, (in dsolve/process_input) invalid initial condition


Any help is appreciated



could you help me solve this error ? I don't understand what it means.


> eq3:=diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*(x(t)-(diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*x(t)+omega[0]^2*X[0])/omega[0]^2) = -omega[0]^2*X[0]:
> dsolve(eq3);
Warning, it is required that the numerator of the given ODE depends on the highest derivative. Returning NULL.




I have a non linear ode with sinosoial term, (sin(x)).

How can we Analyse the system and plot the bifurcation diagram:


Thank you very much for your help.



I have three system of ODE and i would like to solve it using Homotopy perturbation method. Could you please provide to me the code in Maple or the Maple pachage that used to solve it by Homotopy perturbation method ?

I hope to hear you soon




I need you to answer to thos question using maple,

I have x'(t)=e^{2*x(t)}  + alpha -2*x(t),  where alpha is a constant in R.

I would like to compute the equilibrium points and if there is a saddle-node bifurcation.

Thank you very much.


Good day, can any one help in writing maple programme for the finite difference (FD) formulae define to solve this coupled non-linear  ODEs. See it here Thank you

NOTE: please disregard the earlier link.

Hi I am simply trying to rearange equation X, the steoes I want to do are differentiate X w.r.t t then insert equation Y into the new differentiated form but I am stuck, I am not really sure what to try I have used algsubs but it comes up with an error.


here is my code: 


X := diff(x(t), t)+(k1+k2+k3)*x(t)-(k4-k1)*y(t)-k1 = 0;
print(`output redirected...`); # input placeholder
/ d \
|--- x(t)| + (k1 + k2 + k3) x(t) - (k4 - k1) y(t) - k1 = 0
\ dt /
Y := diff(y(t), t)+(k4+k5+k6)*y(t)-(k3-k6)*x(t)-k6 = 0;
print(`output redirected...`); # input placeholder
/ d \
|--- y(t)| + (k4 + k5 + k6) y(t) - (k3 - k6) x(t) - k6 = 0
\ dt /
A := diff(X, t);
print(`output redirected...`); # input placeholder
/ d / d \\ / d \
|--- |--- x(t)|| + (k1 + k2 + k3) |--- x(t)|
\ dt \ dt // \ dt /

/ d \
- (k4 - k1) |--- y(t)| = 0
\ dt /
algsubs(Y, A);
Error, (in algsubs) cannot compute degree of pattern in y(t)


Hi all

I have written the following code in maple to approximate arbitrary functions by hybrid of block-pulse and bernstein functions but it doesn't work properly especially for f(t)=1.0, so what is the matter?



best wishes

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department


I'm dealing with 2nd order ODE on Maple. By using " infolevel 5" Maple tell me that it use Kovacic's algorithm to find the solution. Could anybody tell me how or at least some idea so that I can go on this my self. Following here my ODE

Thank you so much


Hello! How can I find extremes of numeric solution of ODE system obtained using "dsolve"? Can I use something like "extrema" function?

I failed to solve the ODE system shown as follows, where y1(x) and y2(x) are functions of x, ranging from -L/2 to L/2. All the other parameters are constants (A,B,C,F,G). The analytic or numeric solution of y1(x) and y2(x) are wanted.Really appreciate for you experts' help and time!!!


boundary conditions:y1(0)=0, diff(y2(L/2),x$2)=0, D(y2)(0)=0, y2(L/2)=0



2 3 4 5 6 7 8 Last Page 4 of 21