oggsait

80 Reputation

6 Badges

13 years, 62 days

MaplePrimes Activity


These are questions asked by oggsait

i wrote a procedure to  find some functions but at the end  of the procedure gives unexpected string error. i didn't find any way out.

Best regards. 

Download procedureA.mws

Hello i want to sort according to u derivatives (k) system.  And finding determining equations system and solving this system. Thank you very much.  

restart

with(PDEtools)

[CanonicalCoordinates, ChangeSymmetry, CharacteristicQ, CharacteristicQInvariants, ConservedCurrentTest, ConservedCurrents, ConsistencyTest, D_Dx, DeterminingPDE, Eta_k, Euler, FromJet, InfinitesimalGenerator, Infinitesimals, IntegratingFactorTest, IntegratingFactors, InvariantEquation, InvariantSolutions, InvariantTransformation, Invariants, Laplace, Library, PDEplot, PolynomialSolutions, ReducedForm, SimilaritySolutions, SimilarityTransformation, Solve, SymmetrySolutions, SymmetryTest, SymmetryTransformation, TWSolutions, ToJet, build, casesplit, charstrip, dchange, dcoeffs, declare, diff_table, difforder, dpolyform, dsubs, mapde, separability, splitstrip, splitsys, undeclare]

(1)

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

table( [(  ) = u(x, y, t) ] )

(2)

declare(U[])

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

(3)

pde := diff(U[t]-(3/2)*U[x]-6*U[]^2*U[x]+U[x, x, x], x)+U[y, y] = 0

diff(diff(u(x, y, t), t), x)-(3/2)*(diff(diff(u(x, y, t), x), x))-12*u(x, y, t)*(diff(u(x, y, t), x))^2-6*u(x, y, t)^2*(diff(diff(u(x, y, t), x), x))+diff(diff(diff(diff(u(x, y, t), x), x), x), x)+diff(diff(u(x, y, t), y), y) = 0

(4)

NULL

w := phi(x, y, t, U[])

phi(x, y, t, u(x, y, t))

(5)

w*(-12*U[x]^2-12*U[]*U[x, x])+12*w*U[x]^2+12*U[]*w*U[x, x]+(diff(w, x, x))*(-3/2-6*U[]^2)+diff(diff(w, t), x)+diff(w, y, y)+diff(w, x, x, x, x)-lambda*(diff(U[t]-(3/2)*U[x]-6*U[]^2*U[x]+U[x, x, x], x)+U[y, y])

-lambda*(diff(diff(u(x, y, t), t), x)-(3/2)*(diff(diff(u(x, y, t), x), x))-12*u(x, y, t)*(diff(u(x, y, t), x))^2-6*u(x, y, t)^2*(diff(diff(u(x, y, t), x), x))+diff(diff(diff(diff(u(x, y, t), x), x), x), x)+diff(diff(u(x, y, t), y), y))+(D[1, 1, 1, 1](phi))(x, y, t, u(x, y, t))+(D[1, 1, 1, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))+((D[1, 1, 1, 4](phi))(x, y, t, u(x, y, t))+(D[1, 1, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(u(x, y, t), x))+(D[1, 1, 4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), x), x))+((D[1, 1, 1, 4](phi))(x, y, t, u(x, y, t))+(D[1, 1, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))+((D[1, 1, 4, 4](phi))(x, y, t, u(x, y, t))+(D[1, 4, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(u(x, y, t), x))+(D[1, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), x), x)))*(diff(u(x, y, t), x))+2*((D[1, 1, 4](phi))(x, y, t, u(x, y, t))+(D[1, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(diff(u(x, y, t), x), x))+(D[1, 4](phi))(x, y, t, u(x, y, t))*(diff(diff(diff(u(x, y, t), x), x), x))+((D[1, 1, 1, 4](phi))(x, y, t, u(x, y, t))+(D[1, 1, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))+((D[1, 1, 4, 4](phi))(x, y, t, u(x, y, t))+(D[1, 4, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(u(x, y, t), x))+(D[1, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), x), x))+((D[1, 1, 4, 4](phi))(x, y, t, u(x, y, t))+(D[1, 4, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))+((D[1, 4, 4, 4](phi))(x, y, t, u(x, y, t))+(D[4, 4, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(u(x, y, t), x))+(D[4, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), x), x)))*(diff(u(x, y, t), x))+2*((D[1, 4, 4](phi))(x, y, t, u(x, y, t))+(D[4, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(diff(u(x, y, t), x), x))+(D[4, 4](phi))(x, y, t, u(x, y, t))*(diff(diff(diff(u(x, y, t), x), x), x)))*(diff(u(x, y, t), x))+3*((D[1, 1, 4](phi))(x, y, t, u(x, y, t))+(D[1, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))+((D[1, 4, 4](phi))(x, y, t, u(x, y, t))+(D[4, 4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(u(x, y, t), x))+(D[4, 4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), x), x)))*(diff(diff(u(x, y, t), x), x))+3*((D[1, 4](phi))(x, y, t, u(x, y, t))+(D[4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(diff(diff(u(x, y, t), x), x), x))+(D[4](phi))(x, y, t, u(x, y, t))*(diff(diff(diff(diff(u(x, y, t), x), x), x), x))+(D[2, 2](phi))(x, y, t, u(x, y, t))+(D[2, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), y))+((D[2, 4](phi))(x, y, t, u(x, y, t))+(D[4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), y)))*(diff(u(x, y, t), y))+(D[4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), y), y))+(D[1, 3](phi))(x, y, t, u(x, y, t))+(D[3, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))+((D[1, 4](phi))(x, y, t, u(x, y, t))+(D[4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(u(x, y, t), t))+(D[4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), t), x))+((D[1, 1](phi))(x, y, t, u(x, y, t))+(D[1, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))+((D[1, 4](phi))(x, y, t, u(x, y, t))+(D[4, 4](phi))(x, y, t, u(x, y, t))*(diff(u(x, y, t), x)))*(diff(u(x, y, t), x))+(D[4](phi))(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), x), x)))*(-3/2-6*u(x, y, t)^2)+12*u(x, y, t)*phi(x, y, t, u(x, y, t))*(diff(diff(u(x, y, t), x), x))+12*phi(x, y, t, u(x, y, t))*(diff(u(x, y, t), x))^2+phi(x, y, t, u(x, y, t))*(-12*(diff(u(x, y, t), x))^2-12*u(x, y, t)*(diff(diff(u(x, y, t), x), x)))

(6)

k := simplify(%)

-(3/2)*(D[1, 1](phi))(x, y, t, u(x, y, t))+(D[1, 3](phi))(x, y, t, u(x, y, t))+(D[2, 2](phi))(x, y, t, u(x, y, t))+(D[1, 1, 1, 1](phi))(x, y, t, u(x, y, t))+4*(D[1, 1, 1, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)+6*(D[1, 1, 4](phi))(x, y, t, u(x, y, t))*(D[1, 1](u))(x, y, t)+4*(D[1, 4](phi))(x, y, t, u(x, y, t))*(D[1, 1, 1](u))(x, y, t)+(D[4](phi))(x, y, t, u(x, y, t))*(D[1, 1, 1, 1](u))(x, y, t)+2*(D[2, 4](phi))(x, y, t, u(x, y, t))*(D[2](u))(x, y, t)+(D[4](phi))(x, y, t, u(x, y, t))*(D[2, 2](u))(x, y, t)+(D[3, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)+(D[4](phi))(x, y, t, u(x, y, t))*(D[1, 3](u))(x, y, t)-3*(D[1, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)-(3/2)*(D[4](phi))(x, y, t, u(x, y, t))*(D[1, 1](u))(x, y, t)-lambda*(D[1, 3](u))(x, y, t)+(3/2)*lambda*(D[1, 1](u))(x, y, t)-lambda*(D[1, 1, 1, 1](u))(x, y, t)-lambda*(D[2, 2](u))(x, y, t)+6*(D[1, 1, 4, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)^2+4*(D[1, 4, 4, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)^3+(D[4, 4, 4, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)^4+3*(D[4, 4](phi))(x, y, t, u(x, y, t))*(D[1, 1](u))(x, y, t)^2+(D[4, 4](phi))(x, y, t, u(x, y, t))*(D[2](u))(x, y, t)^2+(D[3](u))(x, y, t)*(D[1, 4](phi))(x, y, t, u(x, y, t))-(3/2)*(D[4, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)^2-6*(D[1, 1](phi))(x, y, t, u(x, y, t))*u(x, y, t)^2+12*lambda*u(x, y, t)*(D[1](u))(x, y, t)^2+6*lambda*u(x, y, t)^2*(D[1, 1](u))(x, y, t)+12*(D[1](u))(x, y, t)*(D[1, 4, 4](phi))(x, y, t, u(x, y, t))*(D[1, 1](u))(x, y, t)+6*(D[1](u))(x, y, t)^2*(D[4, 4, 4](phi))(x, y, t, u(x, y, t))*(D[1, 1](u))(x, y, t)+4*(D[1](u))(x, y, t)*(D[4, 4](phi))(x, y, t, u(x, y, t))*(D[1, 1, 1](u))(x, y, t)+(D[3](u))(x, y, t)*(D[4, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)-12*(D[1, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)*u(x, y, t)^2-6*(D[4, 4](phi))(x, y, t, u(x, y, t))*(D[1](u))(x, y, t)^2*u(x, y, t)^2-6*(D[4](phi))(x, y, t, u(x, y, t))*(D[1, 1](u))(x, y, t)*u(x, y, t)^2

(7)

frontend(coeff, [k, U[x]^2]);

0

(8)

frontend(coeff, [k, U[x]*U[x, x]])

Error, invalid input: coeff received O*O, which is not valid for its 2nd argument, x

 

NULL


Download det.eq..mw

 

i have got alot of mixed and high degree derivatives. For example:

u[x]*u[x,t]*eta[x,t]+u[]^2*u[x]*eta[x]+kis(x,y)u[x,t]^2*u[]+eta(x, y)*u[]*u[x]^2+ksi[x,t]*u[x]^2*u[x,t]+......

like this alot of terms

my question is how can i solve divided by the derivative of the u(x,t) partial differential equations system and so  how can i find eta(x,t,u) and ksi(x,t,u) 

How can i sort like below? i want to sort x dependent y is fixed 

 

Original case          sorted case

x   y                         x(i)   y[]

2   1                         1      4

3   2                         2      1

5   3                         3      2

1   4                         4      5

4   5                         5      3

How can I draw following piecewise closed function t in  3-dimensional cartesian space

D124:=max(-x+y+z,x-y+z,x+y-z)=1; 
D134:=max(-x+y+z,x-y+z,-x-y-z)=1;
D123:=max(-x+y+z,x+y-z,-x-y-z)=1;
D234:=max(x+y-z,x-y+z,-x-y-z)=1;

 

t:=piecewise(x>0 and y>0 and z>0,D124,x<0 and y<0 and z>0,D134,x<0 and y>0 and z<0,D123,x>0 and y<0 and z<0,D234);

1 2 3 4 Page 3 of 4