liesymm
makeforms
construct a set of differential forms from a PDE
Calling Sequence
Parameters
Description
Examples
makeforms(eqns, fncs, rootname)
makeforms(eqns, fncs, [Extd])
eqns
-
PDE, or a list or set of PDEs
fcns
list of functions, for example, u(t,x)
rootname
name used to construct the names for the extended coordinates
Extd
name or list of names for constructing the extended variable names
This routine is part of the liesymm package and is loaded via with(liesymm).
Given one or more partial differential equations this routine constructs a set of differential forms which after closure is equivalent to the original system of equations in the sense of Cartan. The forms are obtained by first reducing the system of PDEs to a system of first order equations. It will handle systems which are quasi-linear.
If more than one equation is given, they must be given as a list or a set.
The second argument specifies the dependent and independent variables as in h⁡t,x. If there is more than one dependent variable, these must be specified as a set or list. The Maple expression h,u⁡t,x,y can be used to construct the required list when more than one dependent variable is involved.
The third argument is used to construct the names of the extended variables. If it is the name k then the names k1, k2, k3, ..., kn are constructed as needed. If a list of names is provided then these names are used in place of k1, k2, ..., kn.
The dependent variables (as defined by h,u⁡t,x) are processed in the order given. For each dependent variable we introduce extended variables corresponding to the partial derivatives with respect to t and x in that order. These extended variables are in turn treated as dependent variables until the reduction to first order is achieved. The resulting coordinate list is t,x,h,k1,k2,u,k3,w4 with k1 and k2 being the partials of h with respect to t and x and w3 and w4 being the partials of u with respect to t and x.
with⁡liesymm:
eq1≔Diff⁡h⁡t,x,x,x=Diff⁡h⁡t,x,t
eq1≔∂2∂x2h⁡t,x=∂∂th⁡t,x
makeforms⁡eq1,h⁡t,x,w
d⁡h−w1⁢d⁡t−w2⁢d⁡x,−d⁡w2&ˆd⁡t−w1⁢d⁡t&ˆd⁡x
makeforms⁡eq1,h⁡t,x,p,q
d⁡h−p⁢d⁡t−q⁢d⁡x,−d⁡q&ˆd⁡t−p⁢d⁡t&ˆd⁡x
annul⁡,t,x
∂∂th⁡t,x−p⁡t,x=0,∂∂xh⁡t,x−q⁡t,x=0,∂∂xq⁡t,x−p⁡t,x=0
eq3≔Diff⁡U⁡t,r,r,r+1r⁢Diff⁡U⁡t,r,r−Diff⁡U⁡t,r,t,t=exp⁡−2⁢U⁡t,r⁢Diff⁡C⁡t,r,t2−Diff⁡C⁡t,r,r2
eq3≔∂2∂r2U⁡t,r+∂∂rU⁡t,rr−∂2∂t2U⁡t,r=ⅇ−2⁢U⁡t,r⁢∂∂tC⁡t,r2−∂∂rC⁡t,r2
eq4≔Diff⁡C⁡t,r,r,r+1r⁢Diff⁡C⁡t,r,r−Diff⁡C⁡t,r,t,t=2⁢Diff⁡C⁡t,r,r⁢Diff⁡U⁡t,r,r−Diff⁡C⁡t,r,t⁢Diff⁡U⁡t,r,t
eq4≔∂2∂r2C⁡t,r+∂∂rC⁡t,rr−∂2∂t2C⁡t,r=2⁢∂∂rC⁡t,r⁢∂∂rU⁡t,r−2⁢∂∂tC⁡t,r⁢∂∂tU⁡t,r
makeforms⁡eq3,eq4,U⁡t,r,C⁡t,r,A,B,F,G
d⁡C−A⁢d⁡t−B⁢d⁡r,d⁡U−F⁢d⁡t−G⁢d⁡r,d⁡G&ˆd⁡t+d⁡F&ˆd⁡r+Gr−A2ⅇU2+B2ⅇU2⁢d⁡r&ˆd⁡t,d⁡B&ˆd⁡t+d⁡A&ˆd⁡r+Br−2⁢B⁢G+2⁢A⁢F⁢d⁡r&ˆd⁡t
annul⁡,t,r
∂∂rC⁡t,r−B⁡t,r=0,∂∂tC⁡t,r−A⁡t,r=0,∂∂rU⁡t,r−G⁡t,r=0,∂∂tU⁡t,r−F⁡t,r=0,∂∂rG⁡t,r−∂∂tF⁡t,r+G⁡t,rr−A⁡t,r2ⅇU⁡t,r2+B⁡t,r2ⅇU⁡t,r2=0,∂∂rB⁡t,r−∂∂tA⁡t,r+B⁡t,rr−2⁢B⁡t,r⁢G⁡t,r+2⁢A⁡t,r⁢F⁡t,r=0
See Also
liesymm[annul]
liesymm[close]
liesymm[determine]
liesymm[setup]
with
Download Help Document