/*
wxMaxima 0.7.5 http://wxmaxima.sourceforge.net
Maxima 5.15.0 http://maxima.sourceforge.net
Using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (aka GCL)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) load(”vect”);
(%o1) d:/PROGRA~1/MAXIMA~1.0/share/maxima/5.15.0/share/vector/vect.mac
(%i2) demo(”vect”);
batching #pd:/PROGRA~1/MAXIMA~1.0/share/maxima/5.15.0/share/vector/vect.dem
At the _ prompt, type ‘;’ followed by enter to get next demo
(%i3) if get(’vect,’version)=false then load(vect)
_ ;
(%i4) declare([p,q,f,g],nonscalar)
_ ;
(%i5) f~p.q~g+q~f.p~g+p~q.f~g=0
(%o5) f.(-(-g)~q)~p-f.(-(-g)~p)~q+f.g~p~q=0
_ ;
(%i6) ev(vectorsimp(%),expandall)
(%o6) 0=0
_ ;
(%i7) example:laplacian(%pi*(h+s))=div(3*s*p)
(%o7) %pi*laplacian(s+h)=3*div(p*s)
_ ;
(%i8) ev(vectorsimp(example),expandall:true)
(%o8) %pi*laplacian(s)+%pi*laplacian(h)=3*div(p)*s+3*(p.grad(s))
_ ;
(%i9) parabolic:[[(u^2-v^2)/2,u*v],u,v]
(%o9) [[(u^2-v^2)/2,u*v],u,v]
_ ;
(%i10) scalefactors(parabolic)
_ ;
(%i11) example:express(example)
(%o11)(%pi*(’diff((s+h),v,2)+’diff((s+h),u,2)))/(v^2+u^2)
=(3*(’diff((s*sqrt(v^2+u^2)*p[v]),v,1)+’diff((s*p[u]*sqrt(v^2+u^2)),u,1)))/(v^2+u^2)
_ ;
(%i12) depends([s,h],u,[s,p],v)
_ ;
(%i13) ev(example,diff,factor)
(%o13) (%pi*(’diff(s,v,2)+’diff(s,u,2)+’diff(h,u,2)))/(v^2+u^2)=
(3*(s*v^2*(’diff(p[v],v,1))+s*u^2*(’diff(p[v],v,1))+(’diff(s,v,1))*v^2*p[v]+s*v*p[v]
+(’diff(s,v,1))*u^2*p[v]+(’diff(s,u,1))*p[u]*v^2+(’diff(s,u,1))*u^2*p[u]+s*u*p[u]))/(v^2+u^2)^(3/2)
_ ;
(%i14) example:[(3*u^3*v+2*u*v^3)/(u^2+v^2),(u^4+2*u^2*v^2)/(v^2+u^2)]
(%o14) [(2*u*v^3+3*u^3*v)/(v^2+u^2),(2*u^2*v^2+u^4)/(v^2+u^2)]
_ ;
(%i15) potential(example)
(%o15) u^2*v*sqrt(v^2+u^2)
_ ;
(%o16) d:/PROGRA~1/MAXIMA~1.0/share/maxima/5.15.0/share/vector/vect.dem