Warning, new definition for norm

Warning, new definition for trace

> Lagrange:=proc(X::vector,fX::vector)
local i,j,poli,L,n;
n:=vectdim(X);
poli:=0;
for i from 1 to n do
L:=1;
for j from 1 to n do
if i<>j then
L:=L*evalf((x-X[j])/(X[i]-X[j]))
fi;
od;
poli:=poli+fX[i]*L;
od;
RETURN(simplify(poli));
end:

> LagrangePt:=proc(pt::numeric,X::vector,fX::vector)
local polinomio, res;
polinomio:=Lagrange(X,fX);
res:=subs(x=pt,polinomio);
RETURN(evalf(res));
end: