
function S=area_ellisse(A,B)
% area dell'ellisse Ax^2+By^2=1 con A,B >0. 

if A <= 0 | B <= 0
    S=NaN;
else
    a=sqrt(1/A);
    b=sqrt(1/B);
    S=pi*a*b;
end