Line integral, jacobian...

[salvomic]

here my hypothesis for Jacobian Matrix function...
It should works, I hope without errors (please, see if any).

Input: [f(1), f(2), ...], [x,y,...]

Salvo

Code:


#cas
jacob(args):=
// Matrice jacobiana by Salvo MiccichÈ
// input expr vettoriale, vettore variabili
BEGIN
local argv, argc, mat, f, var, fn, j, k,  gr, vd;
argv:=[args];
argc:=size(argv);
IF argc !=2 THEN
return "Input:[f1(x),f1(y),f1(z)...], [x,y,z,...]"; 
ELSE
f:=argv(1);
var:=argv(2);
fn:=size(f);
vd:=size(var);
mat:=makemat(0,fn,vd);
FOR j FROM 1 TO fn DO // gradienti
gr:=grad(f(j),var);
FOR k FROM 1 TO vd DO // elementi
mat[j,k]:=gr(k);
END; // for k
END; // for j
return mat;
END; // if-else

END;
#end

Now, please, help me to have a short program to calculate surface integrals and flux integral (Int(f.n)), thank you