DIVERGENCE VISUALIZATION PROJECT (flowbox.mws)

Matt Miller, USC

(updated to Release 5 by Doug Meade, 28 Sep 1998)

> with(linalg): with(plots):

Warning, new definition for norm

Warning, new definition for trace

First we compute the divergence of a simple vector field by using the coordinate-dependent formula. You are going to enter other vector fields here later on.

> F:= ( x , y ) -> ( x - y , x + y + 0.001*x^2) ;

> vars:= x , y :

> div_F:= diverge( [ F(x,y) ] , [ vars ] );

> eval( div_F, [ x = 1.5, y = 2.0 ] );

Part I: Local Behavior

Part II: Global Behavior