I am encountering problems solving a system of differential equations.

In the attached file, on the first try the boundary conditions are defined in H: = He, resulting in "Error, (in DEtools / convertsys) unable to compute coeff". On the second try, the boundary conditions are defined in H: = He + 0.00000000000000001 and this works.

What is the possible cause of the problem?

Thank you

problem.mw

Hello,

I am attempting to check solutions to a system of ODEs using odetest. However, odetest only appears to partially substitute the provided solution. Furthermore, it appears to be related to the alphabetical order of the functions.

For instance, here I have two functions, phiL and phiM, that satisfy Laplace's equation and are coupled through the boundary conditions, BCs.

laplace := {-phiL(z) + diff(phiL(z),z$2)=0, -phiM(z) + diff(phiM(z),z$2)=0}:
BCs := {phiL(d1)=0,phiM(-d1)=0,phiL(0)=phiM(0), D(phiL)(0)-D(phiM)(0)=-n}:
sol := {
phiM(z) = n/2/coth(d1)*(cosh(z)+coth(d1)*sinh(z)),
phiL(z) = n/2/coth(d1)*(cosh(z)-coth(d1)*sinh(z))
}:

odetest(sol, laplace union BCs,{phiL(z),phiM(z)});

This returns

{0, 1/2*(2*phiL(0)*coth(d1)-n)/coth(d1), D(phiM)(0)-1/2*n}

Here, phiL(0) and phiM(0) are unevaluated even though the provided solutions are valid there.

Furthermore, while renaming phiL to an alphabetically earlier name (eg, phiJ) causes the corresponding change in the output. However, renaming it to something alphabetically after phiM (eg, phiN) causes the terms in the output to switch. That is, changing phiL to phiN in the above code results in

{0, 1/2*(2*phiM(0)*coth(d1)-n)/coth(d1), D(phiN)(0)+1/2*n}

Therefore, it seems to be related to the way Maple internally stores the list of variables.

Is this a bug? Or is there something I'm missing?

Thanks!

My code below works but I an looking to speed it up to apply to large arrays. I have large integers and I want to store the exponent and first and last few digits.

nelems := 10;  ~ takes about 50mins when n=100,000,000
n := 374894756873546859847556;
op(n);
A := Array(1 .. 4, 1 .. nelems);
length(n);
st := time();
for i to nelems do
    A[1, i] := i^10*n;
    A[2, i] := length(A[1, i]) - 1;
    b := convert(A[1, i], string);
    A[3, i] := parse(b[1 .. 3]);
    A[4, i] := parse(b[-3 .. -1]);
end do;
time() - st;
A;
A[1, -2];
A[2, -2];
A[3, -2];
A[4, -2];

Hello,
I am trying to implement a Squarefree Factorization Algo in Maple. 
For several algorithms (like the one in Victor Shoup- A computational Introduction to Number Theory and Algebra) one is in need of computing the pth root (F_q is the finite field with q=p^k) of a polynomial. I tried stuff like pol^(1/p), root[p](pol), proot(pol,p), but either I get NOROOT as output or just the polynomial with the exponent 1/p, like e.g.

(x^6+(2*tau+1)*x^3+2)^(1/3)

(tau my RootOf for the field extension, F_q with q=3^2)
But that doesnt help me, because I will need a polynomial without any exponent so that i can get the derivate pol'=dpol/dx, and that doenst work if there is any exponent...
cTest := `mod`(root[p](c), p);
cStrich := `mod`(diff(cTest, x), p);
          (x^6+(2*tau+1)*x^3+2)^(1/3)                    
          Error, the modular inverse does not exist

Does anybody know how I can calculate the pth root of a polynomial? (At the moment I am following the Algo that is mentioned on wiki for Squarefree Factorization...)

Or, does anybody know a better Algo for a squarefree Factorization? 

Thanks for any help and any answers!

In dark matter cosmology the following integral in x  does not yet have an analytical solution. I report the integral in

Maple 2018 notation

int(1/sqrt((1+x)^3*a+(1-a)*(1+x)^(3+3*b)), x)

If someone has some smart idea is welcome

I require one of these Weyl + W or Weyl - W (where W is Weyl tensor with hodge star operator, the equation may be found in my code) to be identically equal to zero. The computer will not recognise that it is identally equal unless i input the original p.d.e which I have, and I have also added some partial differentials up to the third order, maybe I need to include up to the fourth order? I am unsure what to add can someone please help me make this identically equal to zero! Thanks in advance for any help!

(I have put colon's after the equations Weyl+W and Weyl-W as the output was too long, but rest assured its not looking identically equal to zero!)
 

Download NEED_11.mw

restart;
ddesys := {diff(pred(t), t) = pred(t)*[prey(t)/(prey(t)+10)-2/3], diff(prey(t), t) = prey(t)*[2-(1/20)*prey(t-tau)]-prey(t)*[pred(t)/(prey(t)+10)]-10, pred(0) = 1, prey(0) = 1};
 

dsn := dsolve[eval(ddesys, tau = 0)];
plots[odeplot](dsn, [[t, prey(t), color = red], [t, pred(t), color = blue]], 0 .. 300, legend = [prey, pred], labels = [t, ""]);
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution 

I do not understand what i mistake, please help me regarding this issue, thanks in advance

This is a question out of curiosity only ...

I can easily import a matrix from an Excel spreadsheet / source. Maple then displays matrix dimensions, data type, storage and order ... which is fine. Now - suppose I wished to import the standard matrix in question, with all entries displayed - can this be done?

I am working on a problem in geometry where I have ended up with a system of nonlinear ODEs in F and G where F and G are functions of a coordinate y, and A and B are both real constants.  I have included the worksheet where I am working on the equations.

Ric_Equations.mw
 

 I would like to solve these in as explicit a way as possible by obtaining expressions for F and G and was wondering how best to do this with Maple (it doesn't really matter if the explicit solutions are messy, as I just need to know that they exist and can be written down in some form).  I have tried using the dsolve command directly but I receive the error message which you can see above.

Sam := Sample(('Uniform')(10, 20), 30)

My question has two steps:

STEP 1:  The multiplication  of is defined as follows

if n<>l, then

.

if n=l and m<=s,

Question 1: I wrote a code for calculating the multiplication  of. Is it right?

The code for Step 1  

restart;

multiply:=proc(n,m,l,s) local g,a: 
a:=unapply(doublefactorial(2*j-1)/factorial(j),j):
g:=unapply((a(m-j)*a(j)*a(s-j)/a(m+s-j))*(2*m+2*s-4*j+1)/(2*m+2*s-2*j+1),j):

if n<>l then 0 else
sqrt((2*m+1)*(2*s+1))*2^(K/2-1).add((g(j)/sqrt(m+s-2*j+1/2))*phi[n, m+s-2*j],j=0..m) 
end if
end proc:
 
n:=2:
l:=2:
m:=1:
s:=1:
multiply(n,m,l,s);

when I compared the results which I got and the results which is given in the book as follows, I think it is right.

Step 2:

We know that the outer product matrix is calculated as follows 

We found the elements of the outer product matrix in Step 1. 

Question 2 : I want to write the elements which are derived in step 1 to the outer product matrix in step 2. In here, the outer product matrix is NxN matrix. N=(M+1).2^(K-1) where K, M are any integers.

Hello,

How I can take variation from left-hand side of  5, and reach to right-hand side of  5. After by using integral by part obtained  7?

Thank you

I use Maple 2015, is it possible to animate a standing wave as shown below?

 A working animation can be found on https://oceanservice.noaa.gov/facts/seiche.html

Download dummy_1.mw