Items tagged with equation


Hi all


I am having a very complicated equation that has the form of 



Of course, the actual equation is more complicate than above. It is just an example. I want to solve the equation in terms of x. And I know that both x and y are real, and they are positive (greater than 0). My question is, how should I specify this when solving the equation?



PS: I try to run the program to the solve the equation (without specifying that they are real and positive), and at the output, it gave me something like:

"RootOf(y^2+2y+.......)".   What is that "RootOf" means? Square-root or what?

Hy all.

I want to solve this equation, with„dd” as numerical result. What do I do wrong? Thanks. Nico

TTot := 70;
TC := 17;
GM := .26;
QMax := 870;
V := 3600*GM*QMax*TTot;
eq := V = int(QMax*exp((-t+TC)/dd)*(1+(t-TC)/TC)^(TC/dd), t = 0 .. TTot);
fsolve(eq, dd);


I have a system of diff equations (see below). I am trying to obtain analytical solution. when I assume that z=wN, I receive such solution. Do anybody have idea if I know that z>wN, does this system has an analytical solution?

diff(K(t), t) = -(1/2)*(Q(t)^2*alpha^2*eta*upsilon-2*eta*alpha*(N*upsilon*w*C[max]-z*alpha*K(t))*Q(t)+N*w*(-2*C[max]*z*eta*alpha*K(t)+upsilon*((-N*w+z)*alpha+N*C[max]^2*w*eta)))*K(t)/((C[max]*w*N-alpha*Q(t))*upsilon*N*w)

diff(Q(t), t) = (1/2)*(-z*(Q(t)^2*alpha^2*eta-2*N*Q(t)*alpha*eta*w*C[max]+w*(w*(eta*C[max]^2-alpha)*N+z*alpha)*N)*K(t)-2*N*upsilon*w*(N*w-z)*(C[max]*w*N-alpha*Q(t)))/((C[max]*w*N-alpha*Q(t))*upsilon*N*w)

K(0) = K0, Q(0) = Q0




sin(xy) = x + y

subs( y(x)=y, solve( diff(subs( y=y(x), (1) ),x), diff(y(x),x) ) );


subs( x(xy)=x, solve( diff(subs( x=x(xy), (1) ),xy), diff(x(xy),xy) ) );


I can not get this answer to come out correctly when using this software please help. the correct answer in the back of the book is 


1- y cos(xy)


x cos(xy) - 1


And is there a way that this program can tutor me on how to get this answer intead of spitting out the answer 

the diff tutor only allows for a one sided equation to be entered.

I have a system of first order diff. equations which I would like to solve symbolically. Unfortunately, Maple does not solve the system. Do anybody have suggestions how can I solve this system (please see below):

diff(S(t), t) = -eta*(C[max]*w*N-alpha*Q(t))*K(t)*S(t)/(w*N*(S(t)+K(t))),

diff(K(t), t) = S(t)*((z*eta*alpha*(C[max]*w*N-alpha*Q(t))*S(t)-upsilon*(eta*alpha^2*Q(t)^2-2*C[max]*w*N*eta*alpha*Q(t)+((-N*w+z)*alpha+N*C[max]^2*w*eta)*N*w))*K(t)^2+(2*((1/2)*z*eta*(C[max]*w*N-alpha*Q(t))*S(t)+N*w*upsilon*(N*w-z)))*S(t)*alpha*K(t)+N*S(t)^2*w*alpha*upsilon*(N*w-z))/((K(t)^2*alpha*z+3*S(t)*K(t)*alpha*z+S(t)*(2*S(t)*z*alpha+upsilon*(C[max]*w*N-alpha*Q(t))))*(S(t)+K(t))*N*w),

diff(Q(t), t) = (-alpha*z*(z*eta*(C[max]*w*N-alpha*Q(t))*K(t)+N*w*upsilon*(N*w-z))*S(t)^2+(-z^2*eta*alpha*(C[max]*w*N-alpha*Q(t))*K(t)^2-(eta*alpha^2*Q(t)^2-2*C[max]*w*N*eta*alpha*Q(t)+N*w*((2*N*w-2*z)*alpha+N*C[max]^2*w*eta))*z*upsilon*K(t)-N*w*upsilon^2*(N*w-z)*(C[max]*w*N-alpha*Q(t)))*S(t)-N*w*z*alpha*upsilon*K(t)^2*(N*w-z))/((2*S(t)^2*alpha*z+(3*z*alpha*K(t)+upsilon*(C[max]*w*N-alpha*Q(t)))*S(t)+K(t)^2*alpha*z)*N*w*upsilon)

where initials conditions are:

S(0) = S0, K(0) = K0, Q(0) = Q0






Hello guys ...

I used a numerically method to solve couple differential equation that it has some boundary conditions. My problem is that some range of answers has 50% error . Do you know things for improving our answers in maple ?

my problem is :



suggestion method by preben Alsholm:

a,b,c,d,e,j,h are constants.suppose some numbers for these constants . I used this code:



I used digits for my code at the first of writting.

please help me ... what should i do?

Hello guys, i have a system of equations ( dynamical system ) which i have its critical points but when i compute its critical points with maple i get different points , i dont know what is wrong . thank you for your time.

Hello all,

I'm an engineering student working on bicycle shifting. After an dynamic study i encounter the equation:


where C is a real constant

More exactly, i would need to find a formule for theta12 using the succesive derivatives of theta10.

I heard that i need to search a numerical solution but,up to now, don't succeed to solve it.

I'm open to all idea.


Hi dear users:

i will plot the equation below abs(y) in terms of x,(note:abs(y) and x is real values),can every body help me?


-32.46753247/(Pi*x^2)+1.053598444*10^8*Pi^2*y/x^2-5.342210338*10^14*Pi^2*y*(2.574000000*10^8*Pi^2-.7700000000*x^2)/((-3.904240733*10^6*x^2+1.305131902*10^15*Pi^2-159.8797200*Pi^2*x^2+2.672275320*10^10*Pi^4+2.391363333*10^(-7)*x^4)*x^2)+1.504285714*10^9*Pi^4*y^3/x^2 = y

i have a non linear equation that depends on three variables e, theta and z.

i have done calculations to calculate e while varying theta and z. theta varied among the vector [0, Pi/4, Pi/3, Pi/2] and z was varying between 1 and 20

when plotting my data it gives the following plot where z is represented on the x-axis and each curve correspond to one theta


i am currently able of fitting one plot to one equation i would like to fit the data points using the nonlinearfit function and to only get one equation for all the plots. is that possible in maple or not


How can I solve a differential equation set of the type,

dy(x)/dx +y^2 =P(x); dP(x)/dx = R(P) numerically

This is the first presentation of updates for the DE and Mathematical Functions programs of Maple 18. It includes several improvements, all in the Mathematical Functions sector, as well as some fixes. The update and instructions for its installation are available on the Maplesoft R&D webpage for DEs and mathematical functions. Some of the items below were mentioned here in Mapleprimes - you are welcome to present suggestions or issues; if possible they will be addressed right away in the next update.

  • Filling gaps in the FunctionAdvisor regarding all the 6 complex components: abs, argument, conjugate, Im, Re, signum, as well as regarding Heaviside (step function), Dirac, min and max.
  • Fix the simplification and differentation rule for doublefactorial
  • Make convert(..., hypergeometric) work the same way as convert(blabla, hypergeom)
  • Implement integral forms for Heaviside(z) and JacobiAM(z, k) via convert(..., Int)
  • Implement appropriate display for the inert %intat function as well as its conversion to the inert Int
  • Make the FunctionAdvisor/DE return not just the PDE system satisfied by f(z, k) = JacobiAM(z, k)and also (new) the ODE satisfied by f(z) = JacobiAM(z, k)
  • Fix conversion rule from Heaviside(z) to Sum
  • Fix unexpected error interruption when differentiating min(...) and max(...) containing more than three arguments
  • Fix issue in simplify/conjugate
  • Improvement in expand/int: factors in disguise are put outside the integration sign
  • Various improvements in the case of multiple integrals involving the Dirac function
  • Make Intat fully inert (before it was evaluating its arguments)
  • Make value of inert indexed objects work

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

I am trying to solve an equation by continuation method (one vairable and one parameter) containing more than 1000 terms and also hight exponents like 800. I want to get omega' values for different 'p' values.The maple file  is appended below. the link is

Many thanks in advance.

Hi,   I want to substitute every  U[jj+1,-1] by U[jj+1,1] in these lines. Many thanks to send your remarks.

for jj from 1 to M do:
  sys[jj] := eval(BTCS_general,j=jj);

subs(U[jj,M+2]=U[jj,M]; U[jj+1,-1]=U[jj+1,1], sys);

  Why in this equation U[jj+1,-1] doesnt changed by U[jj+1,1];
                    s u[i + 1, -1]   4 s u[i + 1, 0]   6 s u[i + 1, 1]
      u[i + 1, 1] + -------------- - --------------- + ---------------
                           4                4                 4       
                          h                h                 h        

           4 s u[i + 1, 2]   s u[i + 1, 3]          
         - --------------- + ------------- = u[i, 1]
                  4                4                
                 h                h                 

Hi ;

I need your help to write the system contains all these equation:

# system 1
# system 2
for j from 1 to N do
u[0, j] := (1/4)*(u[1,j]+u[1,j]+u[0, j-1]+u[0,j+1]-f[0,j]*h^2);
end do;
# System 3
for j from 1 to N do
u[N+1, j] := (1/4)*(u[N, j]+u[N-1,j]+u[N+1, j-1]+u[N+1, j+1]-f[N+1,j]*h^2);
end do;
 eqs := [ seq(seq(Stencil[1](h,i,j,u,f),i=1..N),j=1..N)];

How can collect these system 1 system2, system3 ans system4 in a set with one name.

sys:=[eqs,system3,system2,system2]; ( sys: here contains all the equation).

Thanks you.


First 7 8 9 10 11 12 13 Last Page 9 of 24