way # 1 -- with Gauss ie, int int_S F dot dS = int int int_solid S [ div F ] dV

> restart : with(linalg):

Warning, new definition for norm

Warning, new definition for trace

> F_1 := x;

> F_2 := y;

> F_3 := 3 ;

> F := [ F_1, F_2, F_3 ];

> div_of_F := diverge(F, [x,y,z]);

> volume_of_unit_sphere := 4*Pi / 3 ;

> flux := div_of_F * volume_of_unit_sphere ;