Maple 18 Questions and Posts

These are Posts and Questions associated with the product, Maple 18

Hello,
When I try to put the Gcdex in a procedure and start maplemint, then there occurs an error.

Gcdex(x^2 - 1, x - 2, x ,'s','t') mod 3;

--> works

But:

restart;
a := proc()
    Gcdex(x^2 - 1, x - 2, x ,'s','t') mod 3;
end proc;
maplemint(a);

Then there is an error I don't understand.

Error, (in maplemint/expression) not implemented POLY

By the way I have a fundamental problem to understand, where the values s and t are saved after calling Gcdex (or Quo, Rem, etc.). Till now I thought, that variables s and t are created, but when I declare s, t at the beginning as local variables and start maplemint, then there is something like:

    These parameters have the same name as constants:
      3
    These local variables were used before they were assigned a value:
      r::name, (-x-1)::name, (x-1)::name

So the names of s and t changed, they don't assign a new value? I don't understand that.

 

I having a hard time with defining a vector, in order to store in it some data, then plot it and export it to a file, I copied all what's in the help instructures but it doesn't work everytime, please it's urgent for my PhD thesis !

Hi, I have a big system with 27 polynomial equations in 16 unknowns: f_1=...=f_27=0.  I can store these equations but I cannot calculate a Grobner basis of the ideal  J generated by my polynomials (allocation problem) - I use the library "with(FGb)"-  What interests me is whether my system is minimal in the following sense.

If, for example,  I remove f_1, is the ideal generated by (f_2,...f_27)  J again ? That is to say, is f_1 in the ideal generated by f_2,...,f_27 ? I would like to get an answer "yes" or "no" for each removed  f_i.

My question: can we solve the problem above  without calculating a Grobner basis of J?

Thanks in advance.

 

 

 

 

 

Hello every one,

I'm using Maple to create a compiled program of my model Maplesim. I had a error with "getcompiledProc" command, which is :

Error, (in GetCompiledProc) non-numeric initial condition for `Main.RAB.value`(t): Float(undefined)

the "Main.RAB.value(t)" represent a probe used to visualize an output of my system.

Does anybody had a such error ?
Thank you for your responses.

Larbi ANIBA

Hi, I need help with this little university assignment.

See the thumbnail below.

I have been away from the class because of a staph infection and missed the notes about vectors and Euler's formula.

 

It would be helpful if someone could give a rundown about how to solve the problems or give a general solution.

 

 

-Thanks

Hello, 

I have a PDE system. When I use pdsolve it gets me the messege " pdsolve->Warning: System is inconsistent". Is there a way I can see which equations breaks the system down? 
For this system, it's difficult to see from ayeball where the problem is. 
Thank you! 

test.mw

i copy maple code from notepad to maple in maple window,

there is no error

my function in the code

explicit define parameters are Local type

for example

appendto("...");

func1(aaa)

Local aaa;

 

but when i run cmaple to read the code text file in window 8

it return error

missing operator, syntax error

at Local aaa;

originally 

i have defined

Local aaa, ii;

for ii from 1 to nops(aaa) do

etc.

but it has error too,

then i change to one by one

Local aaa;

Local ii;

still have error at Local aaa;

 

I'm displaying a series of point plots as an animation, and would like to update a displayed parameter as well.  I have a nested list L[t] where there's a set of points for each t, and for each t there's also a numerical value M that I'd like to display.  (In my real problem, L[t] is the number of particles in each of several states, and M is the rms deviation from am algebraic probability distribution.)  

The closest I've gotten (for a simple L and M) is the following, but it displays all of the M values in the legend at once:

with(plots); with(Statistics);

L := [[1, 2, 3, 4, 5, 6], [2, 4, 6, 8, 10, 12]];

M := [1, 2];

display([seq(PointPlot(L[t], legend = M[t]), t = 1 .. 2)], insequence = true)

 

I don't need this to be in the legend.  Is there a way to display only the current value of M for each t?  Thanks very much.

Dear All

In following I tried to find symmetries of certain partial differential equation taken from paper "Group classification and exact solutions of generalized modified Boussinesq equation". But the determining equations are not matching with equations obtained in paper.


with(PDEtools)

DepVars := [f(u(x, t)), u(x, t)]; declare(f(u(x, t)), u(x, t))

[f(u(x, t)), u(x, t)]

 

f(u(x, t))*`will now be displayed as`*f

 

u(x, t)*`will now be displayed as`*u

(1)

PDE1 := diff(u(x, t), t, t)-delta*(diff(u(x, t), x, x, t, t))-(diff(f(u(x, t)), x, x))

diff(diff(u(x, t), t), t)-delta*(diff(diff(diff(diff(u(x, t), t), t), x), x))-((D@@2)(f))(u(x, t))*(diff(u(x, t), x))^2-(D(f))(u(x, t))*(diff(diff(u(x, t), x), x))

(2)

G := [seq(xi[j](x, t, u), j = [x, t]), seq(eta[j](x, t, u), j = [u])]

[xi[x](x, t, u), xi[t](x, t, u), eta[u](x, t, u)]

(3)

declare(G)

eta(x, t, u)*`will now be displayed as`*eta

 

xi(x, t, u)*`will now be displayed as`*xi

(4)

DetSys := DeterminingPDE(PDE1, G, integrabilityconditions = false)

{diff(diff(xi[t](x, t, u), u), u)-(diff(diff(diff(diff(xi[t](x, t, u), u), u), x), x))*delta, (diff(diff(eta[u](x, t, u), x), x))*(diff(f(u), u))+(diff(diff(diff(diff(eta[u](x, t, u), t), t), x), x))*delta-(diff(diff(eta[u](x, t, u), t), t)), -(diff(diff(xi[x](x, t, u), u), u))*(diff(f(u), u))-(diff(diff(diff(diff(xi[x](x, t, u), t), t), u), u))*delta+(diff(diff(f(u), u), u))*(diff(xi[x](x, t, u), u)), 2*(diff(diff(diff(diff(xi[x](x, t, u), u), u), u), x))+2*(diff(diff(diff(diff(xi[t](x, t, u), t), u), u), u))-(diff(diff(diff(diff(eta[u](x, t, u), u), u), u), u)), 2*(diff(diff(diff(xi[x](x, t, u), u), u), x))+2*(diff(diff(diff(xi[t](x, t, u), t), u), u))-(diff(diff(diff(eta[u](x, t, u), u), u), u)), 4*(diff(diff(xi[x](x, t, u), t), x))-2*(diff(diff(eta[u](x, t, u), t), u))+diff(diff(xi[t](x, t, u), t), t), 2*(diff(diff(xi[x](x, t, u), u), x))+2*(diff(diff(xi[t](x, t, u), t), u))-(diff(diff(eta[u](x, t, u), u), u)), diff(diff(xi[x](x, t, u), x), x)-2*(diff(diff(eta[u](x, t, u), u), x))+4*(diff(diff(xi[t](x, t, u), t), x)), -2*(diff(xi[x](x, t, u), x))+(diff(diff(diff(eta[u](x, t, u), u), x), x))*delta-2*(diff(diff(diff(xi[t](x, t, u), t), x), x))*delta, -(diff(diff(diff(xi[x](x, t, u), u), x), x))*delta-4*(diff(diff(diff(xi[t](x, t, u), t), u), x))*delta+2*(diff(diff(diff(eta[u](x, t, u), u), u), x))*delta-2*(diff(xi[x](x, t, u), u)), (diff(diff(f(u), u), u))*eta[u](x, t, u)+(diff(diff(diff(eta[u](x, t, u), t), t), u))*delta-2*(diff(diff(diff(xi[x](x, t, u), t), t), x))*delta+2*(diff(xi[t](x, t, u), t))*(diff(f(u), u)), -2*(diff(xi[t](x, t, u), u))*(diff(f(u), u))+4*delta*(diff(diff(diff(eta[u](x, t, u), t), u), u))-2*(diff(diff(diff(xi[t](x, t, u), t), t), u))*delta-8*(diff(diff(diff(xi[x](x, t, u), t), u), x))*delta, 2*(diff(xi[t](x, t, u), u))*(diff(f(u), u))+2*delta*(diff(diff(diff(eta[u](x, t, u), t), u), u))-(diff(diff(diff(xi[t](x, t, u), t), t), u))*delta-4*(diff(diff(diff(xi[x](x, t, u), t), u), x))*delta, 2*(diff(diff(xi[t](x, t, u), t), u))-(diff(diff(eta[u](x, t, u), u), u))+(diff(diff(diff(diff(eta[u](x, t, u), u), u), x), x))*delta-2*(diff(diff(diff(diff(xi[t](x, t, u), t), u), x), x))*delta, diff(diff(xi[x](x, t, u), u), u)-4*(diff(diff(diff(diff(xi[t](x, t, u), t), u), u), x))*delta+2*(diff(diff(diff(diff(eta[u](x, t, u), u), u), u), x))*delta-(diff(diff(diff(diff(xi[x](x, t, u), u), u), x), x))*delta, 2*(diff(xi[x](x, t, u), u))-8*(diff(diff(diff(xi[t](x, t, u), t), u), x))*delta+4*(diff(diff(diff(eta[u](x, t, u), u), u), x))*delta-2*(diff(diff(diff(xi[x](x, t, u), u), x), x))*delta, -(diff(diff(xi[t](x, t, u), u), u))*(diff(f(u), u))-4*(diff(diff(diff(diff(xi[x](x, t, u), t), u), u), x))*delta-(diff(diff(diff(diff(xi[t](x, t, u), t), t), u), u))*delta+(diff(diff(f(u), u), u))*(diff(xi[t](x, t, u), u))+2*(diff(diff(diff(diff(eta[u](x, t, u), t), u), u), u))*delta, -(diff(diff(xi[t](x, t, u), x), x))*(diff(f(u), u))+diff(diff(xi[t](x, t, u), t), t)-2*(diff(diff(eta[u](x, t, u), t), u))+2*(diff(diff(diff(diff(eta[u](x, t, u), t), u), x), x))*delta-(diff(diff(diff(diff(xi[t](x, t, u), t), t), x), x))*delta, -2*(diff(xi[t](x, t, u), x))*(diff(f(u), u))+2*(diff(xi[x](x, t, u), t))-2*(diff(diff(diff(xi[x](x, t, u), t), x), x))*delta+4*(diff(diff(diff(eta[u](x, t, u), t), u), x))*delta-2*(diff(diff(diff(xi[t](x, t, u), t), t), x))*delta, (diff(diff(diff(diff(eta[u](x, t, u), t), t), u), u))*delta-2*(diff(diff(diff(diff(xi[x](x, t, u), t), t), u), x))*delta+(diff(diff(diff(f(u), u), u), u))*eta[u](x, t, u)+(diff(diff(eta[u](x, t, u), u), u))*(diff(f(u), u))-2*(diff(diff(xi[x](x, t, u), u), x))*(diff(f(u), u))+2*(diff(diff(f(u), u), u))*(diff(xi[t](x, t, u), t)+(1/2)*(diff(eta[u](x, t, u), u))), -(diff(diff(xi[x](x, t, u), x), x))*(diff(f(u), u))+2*(diff(diff(eta[u](x, t, u), u), x))*(diff(f(u), u))+2*(diff(diff(f(u), u), u))*(diff(eta[u](x, t, u), x))+diff(diff(xi[x](x, t, u), t), t)-delta*(diff(diff(diff(diff(xi[x](x, t, u), t), t), x), x))+2*delta*(diff(diff(diff(diff(eta[u](x, t, u), t), t), u), x)), 2*(diff(diff(xi[x](x, t, u), t), u))-2*(diff(diff(xi[t](x, t, u), u), x))*(diff(f(u), u))-2*(diff(xi[t](x, t, u), x))*(diff(diff(f(u), u), u))+4*(diff(diff(diff(diff(eta[u](x, t, u), t), u), u), x))*delta-2*(diff(diff(diff(diff(xi[t](x, t, u), t), t), u), x))*delta-2*(diff(diff(diff(diff(xi[x](x, t, u), t), u), x), x))*delta, diff(diff(diff(diff(xi[t](x, t, u), u), u), u), u), diff(diff(diff(diff(xi[t](x, t, u), u), u), u), x), diff(diff(diff(diff(xi[x](x, t, u), t), u), u), u), diff(diff(diff(diff(xi[x](x, t, u), u), u), u), u), diff(diff(diff(xi[t](x, t, u), u), u), u), diff(diff(diff(xi[t](x, t, u), u), u), x), diff(diff(diff(xi[t](x, t, u), u), x), x), diff(diff(diff(xi[x](x, t, u), t), t), u), diff(diff(diff(xi[x](x, t, u), t), u), u), diff(diff(diff(xi[x](x, t, u), u), u), u), diff(diff(xi[t](x, t, u), u), u), diff(diff(xi[t](x, t, u), u), x), diff(diff(xi[t](x, t, u), x), x), diff(diff(xi[x](x, t, u), t), t), diff(diff(xi[x](x, t, u), t), u), diff(diff(xi[x](x, t, u), u), u), diff(xi[t](x, t, u), u), diff(xi[t](x, t, u), x), diff(xi[x](x, t, u), t), diff(xi[x](x, t, u), u)}

(5)

for EQ in sort([op(DetSys)], length) do EQ = 0 end do

diff(xi[t](x, t, u), u) = 0

 

diff(xi[t](x, t, u), x) = 0

 

diff(xi[x](x, t, u), t) = 0

 

diff(xi[x](x, t, u), u) = 0

 

diff(diff(xi[t](x, t, u), u), u) = 0

 

diff(diff(xi[t](x, t, u), u), x) = 0

 

diff(diff(xi[t](x, t, u), x), x) = 0

 

diff(diff(xi[x](x, t, u), t), t) = 0

 

diff(diff(xi[x](x, t, u), t), u) = 0

 

diff(diff(xi[x](x, t, u), u), u) = 0

 

diff(diff(diff(xi[t](x, t, u), u), u), u) = 0

 

diff(diff(diff(xi[t](x, t, u), u), u), x) = 0

 

diff(diff(diff(xi[t](x, t, u), u), x), x) = 0

 

diff(diff(diff(xi[x](x, t, u), t), t), u) = 0

 

diff(diff(diff(xi[x](x, t, u), t), u), u) = 0

 

diff(diff(diff(xi[x](x, t, u), u), u), u) = 0

 

diff(diff(diff(diff(xi[t](x, t, u), u), u), u), u) = 0

 

diff(diff(diff(diff(xi[t](x, t, u), u), u), u), x) = 0

 

diff(diff(diff(diff(xi[x](x, t, u), t), u), u), u) = 0

 

diff(diff(diff(diff(xi[x](x, t, u), u), u), u), u) = 0

 

diff(diff(xi[t](x, t, u), u), u)-(diff(diff(diff(diff(xi[t](x, t, u), u), u), x), x))*delta = 0

 

4*(diff(diff(xi[x](x, t, u), t), x))-2*(diff(diff(eta[u](x, t, u), t), u))+diff(diff(xi[t](x, t, u), t), t) = 0

 

2*(diff(diff(xi[x](x, t, u), u), x))+2*(diff(diff(xi[t](x, t, u), t), u))-(diff(diff(eta[u](x, t, u), u), u)) = 0

 

diff(diff(xi[x](x, t, u), x), x)-2*(diff(diff(eta[u](x, t, u), u), x))+4*(diff(diff(xi[t](x, t, u), t), x)) = 0

 

2*(diff(diff(diff(xi[x](x, t, u), u), u), x))+2*(diff(diff(diff(xi[t](x, t, u), t), u), u))-(diff(diff(diff(eta[u](x, t, u), u), u), u)) = 0

 

-2*(diff(xi[x](x, t, u), x))+(diff(diff(diff(eta[u](x, t, u), u), x), x))*delta-2*(diff(diff(diff(xi[t](x, t, u), t), x), x))*delta = 0

 

(diff(diff(eta[u](x, t, u), x), x))*(diff(f(u), u))+(diff(diff(diff(diff(eta[u](x, t, u), t), t), x), x))*delta-(diff(diff(eta[u](x, t, u), t), t)) = 0

 

2*(diff(diff(diff(diff(xi[x](x, t, u), u), u), u), x))+2*(diff(diff(diff(diff(xi[t](x, t, u), t), u), u), u))-(diff(diff(diff(diff(eta[u](x, t, u), u), u), u), u)) = 0

 

-(diff(diff(xi[x](x, t, u), u), u))*(diff(f(u), u))-(diff(diff(diff(diff(xi[x](x, t, u), t), t), u), u))*delta+(diff(diff(f(u), u), u))*(diff(xi[x](x, t, u), u)) = 0

 

-(diff(diff(diff(xi[x](x, t, u), u), x), x))*delta-4*(diff(diff(diff(xi[t](x, t, u), t), u), x))*delta+2*(diff(diff(diff(eta[u](x, t, u), u), u), x))*delta-2*(diff(xi[x](x, t, u), u)) = 0

 

2*(diff(xi[x](x, t, u), u))-8*(diff(diff(diff(xi[t](x, t, u), t), u), x))*delta+4*(diff(diff(diff(eta[u](x, t, u), u), u), x))*delta-2*(diff(diff(diff(xi[x](x, t, u), u), x), x))*delta = 0

 

2*(diff(diff(xi[t](x, t, u), t), u))-(diff(diff(eta[u](x, t, u), u), u))+(diff(diff(diff(diff(eta[u](x, t, u), u), u), x), x))*delta-2*(diff(diff(diff(diff(xi[t](x, t, u), t), u), x), x))*delta = 0

 

(diff(diff(f(u), u), u))*eta[u](x, t, u)+(diff(diff(diff(eta[u](x, t, u), t), t), u))*delta-2*(diff(diff(diff(xi[x](x, t, u), t), t), x))*delta+2*(diff(xi[t](x, t, u), t))*(diff(f(u), u)) = 0

 

-2*(diff(xi[t](x, t, u), u))*(diff(f(u), u))+4*delta*(diff(diff(diff(eta[u](x, t, u), t), u), u))-2*(diff(diff(diff(xi[t](x, t, u), t), t), u))*delta-8*(diff(diff(diff(xi[x](x, t, u), t), u), x))*delta = 0

 

2*(diff(xi[t](x, t, u), u))*(diff(f(u), u))+2*delta*(diff(diff(diff(eta[u](x, t, u), t), u), u))-(diff(diff(diff(xi[t](x, t, u), t), t), u))*delta-4*(diff(diff(diff(xi[x](x, t, u), t), u), x))*delta = 0

 

diff(diff(xi[x](x, t, u), u), u)-4*(diff(diff(diff(diff(xi[t](x, t, u), t), u), u), x))*delta+2*(diff(diff(diff(diff(eta[u](x, t, u), u), u), u), x))*delta-(diff(diff(diff(diff(xi[x](x, t, u), u), u), x), x))*delta = 0

 

-2*(diff(xi[t](x, t, u), x))*(diff(f(u), u))+2*(diff(xi[x](x, t, u), t))-2*(diff(diff(diff(xi[x](x, t, u), t), x), x))*delta+4*(diff(diff(diff(eta[u](x, t, u), t), u), x))*delta-2*(diff(diff(diff(xi[t](x, t, u), t), t), x))*delta = 0

 

-(diff(diff(xi[t](x, t, u), x), x))*(diff(f(u), u))+diff(diff(xi[t](x, t, u), t), t)-2*(diff(diff(eta[u](x, t, u), t), u))+2*(diff(diff(diff(diff(eta[u](x, t, u), t), u), x), x))*delta-(diff(diff(diff(diff(xi[t](x, t, u), t), t), x), x))*delta = 0

 

-(diff(diff(xi[t](x, t, u), u), u))*(diff(f(u), u))-4*(diff(diff(diff(diff(xi[x](x, t, u), t), u), u), x))*delta-(diff(diff(diff(diff(xi[t](x, t, u), t), t), u), u))*delta+(diff(diff(f(u), u), u))*(diff(xi[t](x, t, u), u))+2*(diff(diff(diff(diff(eta[u](x, t, u), t), u), u), u))*delta = 0

 

-(diff(diff(xi[x](x, t, u), x), x))*(diff(f(u), u))+2*(diff(diff(eta[u](x, t, u), u), x))*(diff(f(u), u))+2*(diff(diff(f(u), u), u))*(diff(eta[u](x, t, u), x))+diff(diff(xi[x](x, t, u), t), t)-delta*(diff(diff(diff(diff(xi[x](x, t, u), t), t), x), x))+2*delta*(diff(diff(diff(diff(eta[u](x, t, u), t), t), u), x)) = 0

 

2*(diff(diff(xi[x](x, t, u), t), u))-2*(diff(diff(xi[t](x, t, u), u), x))*(diff(f(u), u))-2*(diff(xi[t](x, t, u), x))*(diff(diff(f(u), u), u))+4*(diff(diff(diff(diff(eta[u](x, t, u), t), u), u), x))*delta-2*(diff(diff(diff(diff(xi[t](x, t, u), t), t), u), x))*delta-2*(diff(diff(diff(diff(xi[x](x, t, u), t), u), x), x))*delta = 0

 

(diff(diff(diff(diff(eta[u](x, t, u), t), t), u), u))*delta-2*(diff(diff(diff(diff(xi[x](x, t, u), t), t), u), x))*delta+(diff(diff(diff(f(u), u), u), u))*eta[u](x, t, u)+(diff(diff(eta[u](x, t, u), u), u))*(diff(f(u), u))-2*(diff(diff(xi[x](x, t, u), u), x))*(diff(f(u), u))+2*(diff(diff(f(u), u), u))*(diff(xi[t](x, t, u), t)+(1/2)*(diff(eta[u](x, t, u), u))) = 0

(6)

DetSys1 := dsubs(diff(xi[t](x, t, u), u) = 0, diff(xi[t](x, t, u), x) = 0, diff(xi[x](x, t, u), t) = 0, diff(xi[x](x, t, u), u) = 0, diff(eta[u](x, t, u), u, u) = 0, diff(eta[u](x, t, u), x, u, t) = 0, DetSys)

{0, diff(diff(xi[t](x, t, u), t), t)-2*(diff(diff(eta[u](x, t, u), t), u)), diff(diff(xi[x](x, t, u), x), x)-2*(diff(diff(eta[u](x, t, u), u), x)), -2*(diff(xi[x](x, t, u), x))+(diff(diff(diff(eta[u](x, t, u), u), x), x))*delta, (diff(diff(diff(f(u), u), u), u))*eta[u](x, t, u)+2*(diff(diff(f(u), u), u))*(diff(xi[t](x, t, u), t))+(diff(diff(f(u), u), u))*(diff(eta[u](x, t, u), u)), (diff(diff(f(u), u), u))*eta[u](x, t, u)+(diff(diff(diff(eta[u](x, t, u), t), t), u))*delta+2*(diff(xi[t](x, t, u), t))*(diff(f(u), u)), (diff(diff(eta[u](x, t, u), x), x))*(diff(f(u), u))+(diff(diff(diff(diff(eta[u](x, t, u), t), t), x), x))*delta-(diff(diff(eta[u](x, t, u), t), t)), -(diff(diff(xi[x](x, t, u), x), x))*(diff(f(u), u))+2*(diff(diff(eta[u](x, t, u), u), x))*(diff(f(u), u))+2*(diff(diff(f(u), u), u))*(diff(eta[u](x, t, u), x))}

(7)

for EQ in sort([op(DetSys1)], length) do EQ = 0 end do

0 = 0

 

diff(diff(xi[t](x, t, u), t), t)-2*(diff(diff(eta[u](x, t, u), t), u)) = 0

 

diff(diff(xi[x](x, t, u), x), x)-2*(diff(diff(eta[u](x, t, u), u), x)) = 0

 

-2*(diff(xi[x](x, t, u), x))+(diff(diff(diff(eta[u](x, t, u), u), x), x))*delta = 0

 

(diff(diff(f(u), u), u))*eta[u](x, t, u)+(diff(diff(diff(eta[u](x, t, u), t), t), u))*delta+2*(diff(xi[t](x, t, u), t))*(diff(f(u), u)) = 0

 

(diff(diff(eta[u](x, t, u), x), x))*(diff(f(u), u))+(diff(diff(diff(diff(eta[u](x, t, u), t), t), x), x))*delta-(diff(diff(eta[u](x, t, u), t), t)) = 0

 

(diff(diff(diff(f(u), u), u), u))*eta[u](x, t, u)+2*(diff(diff(f(u), u), u))*(diff(xi[t](x, t, u), t))+(diff(diff(f(u), u), u))*(diff(eta[u](x, t, u), u)) = 0

 

-(diff(diff(xi[x](x, t, u), x), x))*(diff(f(u), u))+2*(diff(diff(eta[u](x, t, u), u), x))*(diff(f(u), u))+2*(diff(diff(f(u), u), u))*(diff(eta[u](x, t, u), x)) = 0

(8)

The third equation in (8) can simplify last equation. This will give us eta[u][x] = 0as f[u, u] is non zero.

NULL


Download [1116]_Symmetries_Determination.mw[1116]_Group_classification_and_exact_solutions_of_generalized_modified_Boussinesq_equation.pdf

Regards

Good morning, I'm a student and I installed Maple 18 on my PC, as a "Single User Profile". Unfortunately, I had to buy another computer, so now I'm trying to install the program on my new one. Even if I enter my Purchase Code and my details, Maple says that no more activations are left in my account, even if the expiration date is on March 2020. Is it possible to activate again the program? 

Thank you,

Laura

hi...please help me for solve this nonlinear equations with pdsolve

thanksoffcenter2.mw

La := .25; Lb := 0.1e-1

h := 0.4e-2

rho := 7900

E := 0.200e12

nu := .3

ve := 5

g := 9.8

M := .5

Z0 := 0.1e-2

K := 5/6

C := sqrt(E/rho)

NULL

 

PDE[1] := diff(u(x, t), x, x)+(diff(w(x, t), x))*(diff(w(x, t), x, x)) = (diff(u(x, t), t, t))/C^2

diff(diff(u(x, t), x), x)+(diff(w(x, t), x))*(diff(diff(w(x, t), x), x)) = 0.3949999999e-7*(diff(diff(u(x, t), t), t))

(1)

PDE[2] := K*(diff(phi(x, t), x)+diff(w(x, t), x, x))/(2*(1+nu))+(diff(w(x, t), x))*(diff(u(x, t), x, x))+(diff(u(x, t), x))*(diff(w(x, t), x, x))+(3/2)*(diff(w(x, t), x, x))*(diff(w(x, t), x))^2 = (diff(w(x, t), t, t))/C^2

.3205128205*(diff(phi(x, t), x))+.3205128205*(diff(diff(w(x, t), x), x))+(diff(w(x, t), x))*(diff(diff(u(x, t), x), x))+(diff(u(x, t), x))*(diff(diff(w(x, t), x), x))+(3/2)*(diff(diff(w(x, t), x), x))*(diff(w(x, t), x))^2 = 0.3949999999e-7*(diff(diff(w(x, t), t), t))

(2)

 

PDE[3] := diff(phi(x, t), x, x)-6*K*(diff(w(x, t), x)+phi(x, t))/(h^2*(1+nu)) = (diff(phi(x, t), t, t))/C^2

diff(diff(phi(x, t), x), x)-240384.6154*(diff(w(x, t), x))-240384.6154*phi(x, t) = 0.3949999999e-7*(diff(diff(phi(x, t), t), t))

(3)

 

 

#####################################

(4)

at x= La

PDE[a1] := diff(u(x, t), x)+(1/2)*(diff(w(x, t), x))^2-M*(g-(diff(u(x, t), t, t))-Z0*(diff(phi(x, t), t, t)))/(E*Lb*h) = 0

diff(u(x, t), x)+(1/2)*(diff(w(x, t), x))^2-0.6125000000e-6+0.6250000000e-7*(diff(diff(u(x, t), t), t))+0.6250000000e-10*(diff(diff(phi(x, t), t), t)) = 0

(5)

PDE[a2] := diff(phi(x, t), x)-12*M*Z0*(g-(diff(u(x, t), t, t))-Z0*(diff(phi(x, t), t, t)))/(E*Lb*h^3) = 0

diff(phi(x, t), x)-0.4593750000e-3+0.4687500000e-4*(diff(diff(u(x, t), t), t))+0.4687500000e-7*(diff(diff(phi(x, t), t), t)) = 0

(6)

PDE[a3] := w(x, t) = 0

w(x, t) = 0

(7)

NULL

############################################

``

at x=0 NULL

(8)

PDE[b1] := u(x, t) = 0 

PDE[b2] := w(x, t) = 0

PDE[b3] := diff(phi(x, t), x) = 0

diff(phi(x, t), x) = 0

(9)

################################################

at t=0 for x= [0,La]

u(x, t) = 0

u(x, t) = 0

(10)

w(x, t) = 0

w(x, t) = 0

(11)

phi(x, t) = 0

phi(x, t) = 0

(12)

diff(phi(x, t), t) = 0

diff(phi(x, t), t) = 0

(13)

diff(w(x, t), t) = 0

diff(w(x, t), t) = 0

(14)

diff(phi(x, t), t, t) = 0

diff(diff(phi(x, t), t), t) = 0

(15)

diff(w(x, t), t, t) = 0

diff(diff(w(x, t), t), t) = 0

(16)

######################################################

at t=0 for x= [0,La)

diff(u(x, t), t) = 0

diff(u(x, t), t) = 0

(17)

diff(u(x, t), t, t) = 0

diff(diff(u(x, t), t), t) = 0

(18)

###################################################

at t=0 for x=La

NULL

diff(u(x, t), t) = -ve

diff(u(x, t), t) = -5

(19)

diff(u(x, t), t, t) = g

diff(diff(u(x, t), t), t) = 9.8

(20)

NULL

NULL

 

Download offcenter2.mw

Hello guys, i would like to do parallel computation in my code written in the Maple18. The question that can help me is:

Given a procedure that compute an function g, where g = f1+f2+f3+f4+f5+f6+f7+f8, i would like to compute all fi at same time.
Now, i´m using " grid:-seq('f[i]',[i=1,2,3,4,5,6,7,8])" and it works very well. However, i think that for my case an better solution should be;
Calculate the f1 in core 1, f2 in core 2, f3 in core 3 ... f8 in core 8 at same time, and after this, to sum all results(f1+f2+f3+..+f8). How i can do this?

Att,

Griffith.

Dear All

It is well known that the package "PDEtools" is helpful in finding infinitesimal transformations for PDEs which I illustrate as follow:


with(PDEtools):

DepVars := [u(x, y, t)]

[u(x, y, t)]

(1)

declare(u(x, y, t)):

u(x, y, t)*`will now be displayed as`*u

(2)

U := diff_table(u(x, y, t)):

PDE1 := U[t, x]+(3/2)*u(x, y, t)*U[x, x]+(3/2)*U[x]^2+(1/4)*U[x, x, x, x]+(3/4)*sigma*U[y, y] = 0:

G := [seq(xi[j](x, y, t, u), j = [x, y, t]), seq(eta[j](x, y, t, u), j = [u])]:

declare(G):

eta(x, y, t, u)*`will now be displayed as`*eta

 

xi(x, y, t, u)*`will now be displayed as`*xi

(3)

DetSys := DeterminingPDE(PDE1, G, integrabilityconditions = false):

pdsolve(DetSys)

{eta[u](x, y, t, u) = (1/9)*(-2*(diff(diff(diff(_F1(t), t), t), t))*y^2-4*(diff(diff(_F2(t), t), t))*y+6*sigma*(-(3/2)*(diff(_F1(t), t))*u+(1/2)*(diff(diff(_F1(t), t), t))*x+diff(_F3(t), t)))/sigma, xi[t](x, y, t, u) = (3/2)*_F1(t)+_C1, xi[x](x, y, t, u) = (1/6)*(-2*(diff(diff(_F1(t), t), t))*y^2-4*(diff(_F2(t), t))*y+3*sigma*((diff(_F1(t), t))*x+2*_F3(t)))/sigma, xi[y](x, y, t, u) = (diff(_F1(t), t))*y+_F2(t)}

(4)

The set (4) gives infinitesimal transformations. How we can write  vector fields corresponding to arbitrary constant C1and arbitrary functions "F1(t), F2(t), F3(t) "?"" 

``


Download Writing_Vector_fields.mw

Regards


Here, I attached my maple code. I need to find root. I am using fsolve. But I am not geting the root. Please any one help me... to find the root.

reatart:NULL``

m1 := 0.3e-1;

0.3e-1

(1)

m2 := .4;

.4

(2)

m3 := 2.5;

2.5

(3)

m4 := .3;

.3

(4)

be := .1;

.1

(5)

rho := .1;

.1

(6)

ga := 25;

25

(7)

a := 3.142;

3.142

(8)

q := .5;

.5

(9)

z[0] := 3;

3

(10)

x[0] := 1.5152;

1.5152

(11)

w[0] := 1.1152;

1.1152

(12)

a1 := be*z[0];

.3

(13)

a2 := be*x[0];

.15152

(14)

a3 := rho*w[0];

.11152

(15)

a4 := rho*z[0];

.3

(16)

a5 := rho*w[0];

.11152

(17)

a6 := rho*z[0];

.3

(18)

b1 := a1*a4*ga+a4*ga*m1;

2.475

(19)

D1 := a1+m1+m2+m3+m4;

3.53

(20)

D2 := a1*m2+a1*m3+a1*m4-a2*ga+a3*ga+m1*m2+m1*m3+m1*m4+m2*m3+m2*m4+m3*m4;

1.92600

(21)

D3 := a1*a3*ga+a1*m2*m3+a1*m2*m4+a1*m3*m4-a2*ga*m1-a2*ga*m4+a3*ga*m1+a3*ga*m4+m1*m2*m3+m1*m3*m4+m2*m3*m4+m1*m2*m3;

1.4499000

(22)

D4 := a1*a3*a4*ga+a1*m2*m3*m4-a2*ga*m1*m4+a3*ga*m1*m4+m1*m2*m3*m4;

.3409200

(23)

G1 := -a1*a6-a6*m1-a6*m2-a6*m3;

-.969

(24)

G2 := -a1*a6*m2-a1*a6*m3+a2*a6*ga-a3*a6*ga+a4*a5*ga-a6*m1*m2-a6*m1*m3-a6*m2*m3;

.549300

(25)

G3 := -a1*a3*a6*ga-a1*a6*m2*m3+a2*a6*ga*m1-a3*a6*ga*m1-a6*m1*m2*m3;

-.3409200

(26)

A1 := w^(4*q)*cos(4*q*a*(1/2))+D1*w^(3*q)*cos(3*q*a*(1/2))+D2*w^(2*q)*cos(2*q*a*(1/2))+D3*w^q*cos((1/2)*q*a)+D4;

-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200

(27)

B1 := w^(4*q)*sin(4*q*a*(1/2))+D1*w^(3*q)*sin(3*q*a*(1/2))+D2*w^(2*q)*sin(2*q*a*(1/2))+D3*w^q*sin((1/2)*q*a);

-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5

(28)

A2 := -w^(3*q)*a6*cos(3*q*a*(1/2))+G1*w^(2*q)*cos(2*q*a*(1/2))+G2*w^q*cos((1/2)*q*a)+G3;

.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200

(29)

B2 := -w^(3*q)*a6*sin(3*q*a*(1/2))+G1*w^(2*q)*sin(2*q*a*(1/2))+G2*w^q*sin((1/2)*q*a);

-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5

(30)

C := .27601200;

.27601200

(31)

Q1 := 4*C^2*(A2^2+B2^2);

.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2

(32)

Q2 := -4*C*A2*(A1^2-A2^2+B1^2-B2^2-C^2);

-1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)

(33)

Q3 := (A1^2-A2^2+B1^2-B2^2-C^2)^2-4*C^2*B2^2;

((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)^2-.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2

(34)

V := simplify(-4*Q1*Q3+Q2^2);

-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2)

(35)

x := (-Q2+sqrt(V))/(2*Q1);

(1/2)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)

(36)

E := -2*A1*C*x-A1^2+A2^2-B1^2+B2^2-C^2;

-.2760120000*(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)-(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2+(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2-(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2+(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1

(37)

y := -E/(2*C*B1);

-1.811515442*(-.2760120000*(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)-(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2+(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2-(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2+(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)/(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)

(38)

``

fsolve(x^2+y^2 = 1, w)

fsolve((1/4)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))^2/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)^2+3.281588197*(-.2760120000*(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)*(1.10404800*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)*((-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2-(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2-(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)+(-0.1e-12-0.5481797400e-1*w^11-40.93358002*w^(19/2)-212.0102604*w^(17/2)-1.048226159*w^(21/2)-8.667039897*w^10-119.4464160*w^9-208.1803245*w^8-54.3436016*w^7-38.4722894*w^6+2.67061391*w^5-2.29413863*w^4-.136247212*w^2+.899997750*w^3+0.1e-10*w^(1/2)-0.150073928e-1*w^(3/2)+0.54469063e-2*w-2.53869438*w^(11/2)-2.40374793*w^(9/2)-84.14780373*w^(15/2)-86.62603442*w^(13/2)+2.023073705*w^(7/2)-0.6906749e-2*w^(5/2))^(1/2))/(.3047304966*(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2+.3047304966*(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2)-(-.9999999170*w^2.0-2.496849400*w^1.5-0.3922745903e-3*w^1.0+1.025129710*w^.5+.3409200)^2+(.2121968329*w^1.5+0.1973593344e-3*w^1.0+.3883741982*w^.5-.3409200)^2-(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2+(-.2120672160*w^1.5-.9689999799*w^1.0+.3884533076*w^.5)^2-0.7618262414e-1)^2/(-0.4073463989e-3*w^2.0+2.495324242*w^1.5+1.925999960*w^1.0+1.025338523*w^.5)^2 = 1, w)

(39)

``

 

Download root.mw

1 2 3 4 5 6 7 Last Page 3 of 61