Items tagged with constant constant Tagged Items Feed

When solving ODE with dsolve, I don't want to use "_C1" as the name of constant, I want to specify by myself, how can I do?


Anyone could help me in solving the following system of equations to get constants C1, C2, C3 and C4. MALPE give me this "soution may have been lost".  The MAPLE sheet is also attached.



Eq1:=simplify(C3*exp(-(1/4)*(C2*(x^2-2*0)+sqrt(C2*(x^2-2*0)^2+4*M*(x^2-2*0)*w1*(x^2-2*0)))/w1)+C4*exp((1/4)*(-C2*(x^2-2*0)+sqrt(C2*(x^2-2*0)^2+4*M*(x^2-2*0)*w1*(x^2-2*0)))/w1)-U) = 0;

C3*exp(-(1/4)*(C2*x^2+(x^4*(4*M*w1+C2))^(1/2))/w1)+C4*exp(-(1/4)*(C2*x^2-(x^4*(4*M*w1+C2))^(1/2))/w1)-U = 0


Eq2:=simplify(exp(-(1/4)*(C2+sqrt(C2^2+4*M*w1))*(x^2-2*0)/w1)*C3*x+exp((1/4)*(-C2+sqrt(C2^2+4*M*w1))*(x^2-2*0)/w1)*C4*x+C2-V-z) = 0;

exp(-(1/4)*(C2+(C2^2+4*M*w1)^(1/2))*x^2/w1)*C3*x+exp(-(1/4)*(C2-(C2^2+4*M*w1)^(1/2))*x^2/w1)*C4*x+C2-V-z = 0


Eq3:=simplify((-2*w2*w5*ln(C3*exp(-(1/2)*sqrt(w2*w4*(w2*w4+w3*w6))*C2*(x^2-2*0)/(w2*w4*w5))-C4)*sqrt(w2*w4*(w2*w4+w3*w6))+w2*w5*(-w2*w4+sqrt(w2*w4*(w2*w4+w3*w6)))*ln(exp(-(1/2)*sqrt(w2*w4*(w2*w4+w3*w6))*C2*(x^2-2*0)/(w2*w4*w5)))+C1*w3*w6*sqrt(w2*w4*(w2*w4+w3*w6)))/(sqrt(w2*w4*(w2*w4+w3*w6))*w3*w6)-1)= 0;

(-ln(exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5)))*w2^2*w4*w5+C1*w3*w6*(w2*w4*(w2*w4+w3*w6))^(1/2)-2*w2*w5*ln(C3*exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))-C4)*(w2*w4*(w2*w4+w3*w6))^(1/2)+ln(exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5)))*(w2*w4*(w2*w4+w3*w6))^(1/2)*w2*w5-(w2*w4*(w2*w4+w3*w6))^(1/2)*w3*w6)/((w2*w4*(w2*w4+w3*w6))^(1/2)*w3*w6) = 0


Eq4:= simplify((-C2*x^2*w2*w4-.50*C2*x^2*w3*w6+sqrt(w2*w4*(w2*w4+w3*w6))*C2*x^2+2.*w2*w4*w5*ln(w3^4*w6^2*(C3^2*exp(-1.0*sqrt(w2*w4*(w2*w4+w3*w6))*C2*x^2/(w2*w4*w5))-2*C3*exp(-.5*sqrt(w2^2*w4^2+w2*w3*w4*w6)*C2*x^2/(w2*w4*w5))*C4+C4^2)/(w2*w4*(w2*w4+w3*w6)*C2^2))-5.544000000*w2*w4*w5-w3^2*w6)/(w3^2*w6)) = 0;

(-C2*x^2*w2*w4-.5000000000*C2*x^2*w3*w6+(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2+2.*w2*w4*w5*ln(w3^4*w6^2*(C3^2*exp(-(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))-2.*C3*exp(-.5*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))*C4+C4^2)/(w2*w4*(w2*w4+w3*w6)*C2^2))-5.544000000*w2*w4*w5-w3^2*w6)/(w3^2*w6) = 0


solve({Eq1, Eq2, Eq3,Eq4}, {C1, C2, C3,C4});

Warning, solutions may have been lost







There is an heuristic-free algorithm, designed by Greg Reid to find structure constant of finite dimension Lie algebra of symmetries, admitted by differential equation. Details could be found in his paper:


The main problem - I couldn't find implementation of it in Maple...

Hi all

How can I construct following triangular block matrix?



P is (N*M)*(N*M) and E and H are M*M matrices. tf is known constant.

thanks in advance


Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

I would like to announce a new unofficial record computation of the MRB constant that was finished on Sun 21 Sep 2014 18:35:06.

I really would like to see someone beat it with Maple!

It took 1 month 27 days 2 hours 45 minutes 15 seconds. I computed 3,014,991 digits of the MRB constant, (confirming my previous 2,00,000 or more digit computation was actually accurate to 2,009,993 digits), with Mathematica 10.0. I Used my version of Richard Crandall's code:



(*Fastest (at MRB's end) as of 25 Jul 2014.*)


prec = 3000000;(*Number of required decimals.*)ClearSystemCache[];

T0 = SessionTime[];

expM[pre_] := 

  Module[{a, d, s, k, bb, c, n, end, iprec, xvals, x, pc, cores = 12, 

    tsize = 2^7, chunksize, start = 1, ll, ctab, 

    pr = Floor[1.005 pre]}, chunksize = cores*tsize;

   n = Floor[1.32 pr];

   end = Ceiling[n/chunksize];

   Print["Iterations required: ", n];

   Print["end ", end];

   Print[end*chunksize]; d = ChebyshevT[n, 3];

   {b, c, s} = {SetPrecision[-1, 1.1*n], -d, 0};

   iprec = Ceiling[pr/27];

   Do[xvals = Flatten[ParallelTable[Table[ll = start + j*tsize + l;

        x = N[E^(Log[ll]/(ll)), iprec];

        pc = iprec;

        While[pc < pr, pc = Min[3 pc, pr];

         x = SetPrecision[x, pc];

         y = x^ll - ll;

         x = x (1 - 2 y/((ll + 1) y + 2 ll ll));];(*N[Exp[Log[ll]/ll],

        pr]*)x, {l, 0, tsize - 1}], {j, 0, cores - 1}, 

       Method -> "EvaluationsPerKernel" -> 4]];

    ctab = ParallelTable[Table[c = b - c;

       ll = start + l - 2;

       b *= 2 (ll + n) (ll - n)/((ll + 1) (2 ll + 1));

       c, {l, chunksize}], Method -> "EvaluationsPerKernel" -> 2];

    s += ctab.(xvals - 1);

    start += chunksize;

    Print["done iter ", k*chunksize, " ", SessionTime[] - T0];, {k, 0,

      end - 1}];

   N[-s/d, pr]];

t2 = Timing[MRBtest2 = expM[prec];]; DateString[]


MRBtest2 - MRBtest2M



I used a six core Intel(R) Core(TM) i7-3930K CPU @ 3.20 GHz 3.20 GHz with 64 GB of RAM of which only 16 GB was used.

t2 From the computation was {1.961004112059*10^6, Null}.




For all real a, the partial sums sn= sum((-1)^k (k^(1/k) -a), k=1..n) are bounded so that their limit points form an interval [-1.+  the MRB constant +a, MRB constant] of length 1-a, where the MRB constant is limit(sum((-1)^k*(k^(1/k)), k = 1 ..2*N),N=infinity).

For all complex z, the upper limit point of  sn= sum((-1)^k (k^(1/k) -z), k=1..n) is the  the MRB constant.

We see that maple knows the basics of this because when we enter sum((-1)^k*(k^(1/k)-z), k = 1 .. n) 

maple gives

sum((-1)^k*(k^(1/k)-z), k = 1 .. n)

I have some data for a model in MapleSim that I would like to use a time look up table with.  I've found that the two options for interpolation are linear and 1st derivative, but the data was intended to be interpretted as piecewise constant.  Is there any way to acheive this option in MapleSim?

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?

I need to get rid of the type "constant" for gamma. In Maple, gamma is defined as Euler's constant by default. While it is easy enough to unprotect(gamma) and then get rid of its value, Maple will refuse to solve an equation for gamma, as it remains of type constant even after deassigning it. So I need to regain gamma as a variable.

Some may feel this is an unwise thing to do. But it actually is not: I am writing a document involving physics, and gamma is the accepted symbol for the relativistic energy. I cannot avoid using that, lest mass confusion ensues (this involves students). I really don't want to write gammar instead. Euler's constant, otoh, does not figure at all in my document.

Note that I need a solution that works in Maple 15 and later as I am working in a heterogeneous environment as far as Maple versions are concerned.

Thanks in advance,

Mac Dude




Could not evaluate numerical integral with constant in it. I use method = _cuhre.  Maple print solution like this:


Int(Int(Int(max(0., (0.9483573506e-3*(-1.*sin(a)*cos(w)-1.*cos(a)*sin(w)*sin(b)))*cos(a)*cos(b)^2*(-58.5*signum(cos(b)*sin(w)*sin(b))*kk+200.*cos(b)*sin(w)*sin(b))*Heaviside(-58.5+200.*cos(b)*sin(w)*sin(b)*signum(cos(b)*sin(w)*sin(b))*kk)*Heaviside(1.-.98*cos(b)^2)/sqrt(1.-.98*cos(b)^2)), a = 0. .. 6.283185308), b = 0. .. 1.570796327), w = 0. .. 6.283185308)+Int(Int(Int(min(0., (0.9483573506e-3*(-1.*sin(a)*cos(w)-1.*cos(a)*sin(w)*sin(b)))*cos(a)*cos(b)^2*(-58.5*signum(cos(b)*sin(w)*sin(b))*kk+200.*cos(b)*sin(w)*sin(b))*Heaviside(-58.5+200.*cos(b)*sin(w)*sin(b)*signum(cos(b)*sin(w)*sin(b))*kk)*Heaviside(1.-.98*cos(b)^2)/sqrt(1.-.98*cos(b)^2)), a = 0. .. 6.283185308), b = 0. .. 1.570796327), w = 0. .. 6.283185308)

How could it be taken.




Say, I got an expression that depends on two variables, x and y. How can I tell Maple, that y is actually just a (real) constant, so y does not depend on x?

Because when I apply a differentiation with the "D" - command, it would always also write out expressions, where y is differentiated w.r.t. x.


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...

Hello, I have a problem with plotting graphs with a constant.. I have one system of equations and one another equation:






I need to compare theit curves.

k and C are constants. I thought it is alright if I would remove constants, BUT there is one in the power, so I have no idea what to do.. 

I just tried to plot the system, but it doesn't work either:

Thanks a lot Markiyan, and I am sorry I may have overlooked this already answered question. But, now I have just one more thing to know. On using AllSolutions = true in solve I get the answer as x = _z1~ ∏ / k, I then try to subs my fav symbol N instead of _Z1~, but Maple doesn't take it! It keeps it as it is. What kind of constant is this _Z1~, I had thought it to be the same as _C1, which is used by Maple while solving any ODE.

Reason why I want to...

As many of you know now the MRB constant = sum((-1)^n*(n^(1/n)-1),n=1..infinity).

Here are some equations involving various forms of that summation.

The first one involves convergent series and is too obvious. The others involve divergent series.

The last two, however, are new!


Let c=MRB constant and a, c~, x, and y = any number.


sum((-1)^n*(c~*n^(1/n)-c~),n=1..infinity)= c*c~.

evalf(sum((-1)^n*(n^(1/n)-a),n=1..infinity)) gives c-1/2*(1-a).

evalf(sum((-1)^n*(x*n^(1/n)+y*n),n=1..infinity)) gives (c-1/2)*x-1/4*y.

And it appears that

evalf(sum((-1)^n*(x*n^(1/n)-a),n=1..infinity)) gives (c - 1/2)*x + 1/2*a.

1 2 3 4 Page 1 of 4