Digi Area Group - Math Tools for Professionals
   Maple and Mathematica packages - math tools for professionals

atlas™ - modern differential geometry for Maple™

> Features List & Examples
> Template Worksheets
> Screenshots
> Documentation & Downloads
> License & Pricing
> Buy Online

 
 
 
 
Google

Features List & Examples  |  Introduction  |  Dimension  |  Indexing  |  Forms  |  Metric  |  atlasWizard - Maplet™

atlas[div] - divergence operator

Calling Sequence:

     div(expr)

Parameters:

         expr - any vector expression.

Description:

  •  The div procedure calculates the divergence operator on a vector expression.

Examples:
restart:
with(atlas):

Declare constants:
Constants(Lambda);

{Catalan, _Z, Pi, I, Lambda, -I}

Declare functions:
Functions(f=f(x,y),h=h(f));

{h, f}

Declare p-forms:
Forms(e[k]=1,xi=1);

{e[k], xi}

Declare vectors:
Vectors(E[j],X,Y,Z);

{E[j], X, Y, Z}

By definition:
'div(X+Y)'=div(X+Y);
'div(F*X)'=div(F*X);

div(X+Y) = div(X)+div(Y)

div(F*X) = F*div(X)+iota[X](d(F))

Declare coframe:
Coframe(e[1]=d(x),e[2]=d(y));

[e[1] = d(x), e[2] = d(y)]

Declare frame:
Frame(E[k]);

[E[1] = Diff(``,x), E[2] = Diff(``,y)]

Declare metric:
Metric(g=4*(d(x)&.d(x)+d(y)&.d(y))/(1+Lambda*(x^2+y^2))^2);

g = 4*(`&.`(e[1],e[1])+`&.`(e[2],e[2]))/(1+Lambda*(x^2+y^2))^2

Calculate connection:
Connection(omega);

omega[i,j]

Divergence of frame vectors:
'div(E[1])'=div(E[1]);

div(E[1]) = -4*Lambda*x/(1+Lambda*x^2+Lambda*y^2)

Divergence of vector X:
'div(X)'=div(ToBasis(X));

div(X) = -4*iota[X](e[1])*Lambda*x/(1+Lambda*x^2+Lambda*y^2)+iota[E[1]](d(iota[X](e[1])))-4*iota[X](e[2])/(1+Lambda*x^2+Lambda*y^2)*Lambda*y+iota[E[2]](d(iota[X](e[2])))

Divergence of "rotation" vector:
'div(y*E[1]-x*E[2])'=div(y*E[1]-x*E[2]);

div(y*E[1]-x*E[2]) = 0

See Also:

atlas , atlas[d] , atlas[cov] , atlas[`&.`] , atlas[`&^`] , atlas[iota] .