kill(all);

definisco le funzioni di forma per l'iso4

N : [
(1-xi)*(1-eta)/4,
(1+xi)*(1-eta)/4,
(1+xi)*(1+eta)/4,
(1-xi)*(1+eta)/4
];

creo una matrice le cui due righe sono costituite
dalle liste di derivate di N in xi ed eta

dN : matrix( 
    diff(N,xi,1),
    diff(N,eta,1)
);

converto in codice fortran, da copiaincollarsi

fortran(dN);