## solving a PDE system with constant parameters ...

Hi guys!

I have a PDE system. The mayority of the equations are equal to zero, but two of them are:

where a, b, c, d are CONSTANT parameters. I know that if a=b=d=c=1 the system is inconsistent. But I also know that if a=-1, b=d=0 and c=1 the system is consistent and exist the solution. I wanna know if there's a way to ask maple to find another selections of my parameter that make my PDE consistent and what it's the solution for that selection of a,b,c,d.

Here's my PDE system (sys2).

test.mw

thank you so much for your time!

## Handling of scientific sonstants...

I understand that if I want to use a scientific constant in a Maple worksheet or document, I first have to declare it. Here is an example for the speed of light:

However if this worksheet is re-executed (e.g. using !!!) this happens:

Obviously there is a kind of recursion for c which is listed in the Variable Palette as a name for a variable but which is also the scientific constant. The problem can be solved e.g. by associating another name to the constant and unassigning the old definition in the Variable Palette:

This works but it makes the document harder to read since everybody is used to lowercase c for the speed of light. The same is of course true for other scientific constants.

Is there a way to use scientific constants in their "usual" notation in Maple, including using the unit of the constant?

## Problem with delta as parameter...

I want to ask., I put delta as my constant in maple program and I want the answer are in delta as well., but the thing is., when running., it let delta=0, delta=-1, and delta=delta.,
the condition is we cannot let delta=1 or delta=0 because it is just same for s5 and s7.,.(delta is refer to the s8). How can I get answer as delta? with the condition? here I attach my maple programme..

> derivation := proc (A, n)
local i, j, k, t, s5, s7, s8, m, D,
sols5, sols7, sols8, eqns5, eqns7, eqns8,
BChange5, BChange7, BChange8; eqns5 := {}; eqns7 := {}; eqns8 := {};
D := matrix(n, n);
BChange5 := matrix(n, n); BChange7 := matrix(n, n); BChange8 := matrix(n, n);
for i to n do for j to n do for m to n do
s5 := sum(0*A[i, j, k]*D[m, k], k = 1 .. n)-(sum(A[k, j, m]*D[k, i]+A[i, k, m]*D[k, j], k = 1 .. n));
s7 := sum(0*A[i, j, k]*D[m, k], k = 1 .. n)-(sum(A[k, j, m]*D[k, i]+0*A[i, k, m]*D[k, j], k = 1 .. n));
s8 := sum(0*A[i, j, k]*D[m, k], k = 1 .. n)-(sum(A[k, j, m]*D[k, i]+delta*A[i, k, m]*D[k, j], k = 1 .. n));
eqns5 := `union`(eqns5, {s5}); eqns7 := `union`(eqns7, {s7}); eqns8 := `union`(eqns8, {s8})
end do end do end do;
sols5 := [solve(eqns5)]; sols7 := [solve(eqns7)]; sols8 := [solve(eqns8)];
t := nops(sols5); t := nops(sols7); t := nops(sols8);
for i to t do for j to n do for k to n do
BChange5[k, j] := subs(sols5[i], D[k, j]);
BChange7[k, j] := subs(sols7[i], D[k, j]);
BChange8[k, j] := subs(sols8[i], D[k, j])
end do end do;
print("eqns&Assign;", eqns5); print("sols:=", sols5); print("BChange5:=", BChange5);
print("eqns&Assign;", eqns7); print("sols:=", sols7); print("BChange8:=", BChange7);
print("eqns&Assign;", eqns8); print("sols:=", sols8); print("BChange8:=", BChange8)
end do end proc;

> AS1 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 1, 2) = 1]);
> derivation(AS1, 2);

> AS2 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 1, 1) = 1, (1, 2, 2) = 1]);
> derivation(AS2, 2);

> AS3 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 1, 1) = 1, (2, 1, 2) = 1]);
> derivation(AS3, 2);

> AS4 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 1, 1) = 1, (2, 2, 2) = 1]);
> derivation(AS4, 2);

> AS5 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 1, 1) = 1, (1, 2, 2) = 1, (2, 1, 2) = 1]);
> derivation(AS5, 2);

> AS1 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 3, 2) = 1, (3, 1, 2) = 1]);
> derivation(AS1, 3);

> AS2 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 3, 2) = 1, (3, 1, 2) = alpha]);
> derivation(AS2, 3);

> AS3 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 1, 2) = 1, (1, 2, 3) = 1, (2, 1, 3) = 1]);
> derivation(AS3, 3);

> AS4 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 3, 2) = 1, (2, 3, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS4, 3);

> AS5 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(2, 3, 2) = 1, (3, 1, 1) = 1, (3, 3, 3) = 1]);
> derivation(AS5, 3);

> AS6 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(3, 1, 2) = 1, (3, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS6, 3);

> AS7 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 2, 1) = 1, (2, 2, 2) = 1, (3, 1, 1) = 1, (3, 3, 3) = 1]);
> derivation(AS7, 3);

> AS8 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 3, 1) = 1, (2, 3, 2) = 1, (3, 1, 1) = 1, (3, 3, 3) = 1]);
> derivation(AS8, 3);

> AS9 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(2, 3, 2) = 1, (3, 1, 1) = 1, (3, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS9, 3);

> AS10 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 3, 1) = 1, (2, 3, 2) = 1, (3, 1, 1) = 1, (3, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS10, 3);

> AS11 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 3, 2) = 1, (2, 3, 2) = 1, (3, 1, 2) = 1, (3, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS11, 3);

> AS12 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 1, 2) = 1, (1, 3, 1) = 1, (2, 3, 2) = 1, (3, 1, 1) = 1, (3, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS12, 3);

> AS13 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 1, 1) = 1, (2, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS13, 3);

> AS14 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 2, 1) = 1, (2, 1, 1) = 1, (2, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS14, 3);

> AS15 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 2, 1) = 1, (2, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS15, 3);

> AS16 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(2, 1, 1) = 1, (2, 2, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS16, 3);

> AS17 := array(sparse, 1 .. 3, 1 .. 3, 1 .. 3, [(1, 1, 2) = 1, (3, 3, 3) = 1]);
> derivation(AS17, 3);
>

## How to make constants?...

_C1, _C2, _C3 are constant, how to set them constant, to make diff(_C2) = 0 etc

eval(simplify(subs(a=_C1,subs(b=1/(diff(c(t), t)),subs(c=_C2+_C3*exp(-t),eq2)))));

(diff(_C2(t), t)+(diff(_C3(t), t))*(exp(-t))(t)+_C3(t)*(diff((exp(-t))(t), t)))*(-(diff(_C1(t), t))*(diff(diff((c(t))(t), t), t))/(diff((c(t))(t), t))^2+(diff(s(t), t))/(diff((c(t))(t), t)))
a

## Constants in dsolve ...

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?

## Problem with system of equations...

Hi

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.

 > restart:
 > 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;
 (1)
 > 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;
 (2)
 > 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;
 (3)
 > 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;
 (4)
 > solve({Eq1, Eq2, Eq3,Eq4}, {C1, C2, C3,C4});
 >

## Structure constant...

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:

http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=2316480

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

## triangular block matrix...

Hi all

How can I construct following triangular block matrix?

where

and

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

Ph.D Candidate

Applied Mathematics Department

## New record computation

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

DateString[]

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[]

Print[MRBtest2]

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

## Special bounded partial sums

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(,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

marvinrayburns.com

## Piecewise Constant Interpolation Table...

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?

## How can i improve my nummerical solution?...

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 :

a*Φ''''(x)+b*Φ''(x)+c*Φ(x)+d*Ψ''(x)+e*Ψ(x):=0

d*Φ''(x)+e*Φ(x)+j*Ψ''(x)+h*Ψ(x):=0

suggestion method by preben Alsholm:

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

VR22:=0.1178*diff(phi(x),x,x,x,x)-0.2167*diff(phi(x),x,x)+0.0156*diff(psi(x),x,x)+0.2852*phi(x)+0.0804*psi(x);
VS22:=0.3668*diff(psi(x),x,x)-0.0156*diff(phi(x),x,x)-0.8043*psi(x)-0.80400*phi(x);
bok:=evalf(dsolve({VR22=0,VS22=0}));

PHI,PSI:=op(subs(bok,[phi(x),psi(x)]));
Eqs:={eval(PHI,x=1.366)=1,eval(diff(PHI,x),x=1.366)=0,eval(PHI,x=-1.366)=1,eval(diff(PHI,x),x=-1.366)=0,
eval(PSI,x=1.366)=1,eval(PSI,x=1.366)=1};
C:=fsolve(Eqs,indets(%,name));
eval(bok,C);
SOL:=fnormal(evalc(%));

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

## How to get rid of a constant?...

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.

Mac Dude

## Numerical Integration...

Hello.

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.

Thank!!!