function demo_II
% demo sul calcolo dell'area dell'ellisse Ax^2+By^2=1 con A,B >0.


fid=fopen('risultati.txt','w');

for esempio=1:10
    A=0.5-rand(1);
    B=0.5-rand(1);
    S=area_ellisse(A,B);
    fprintf(fid,'\n \n Esempio %3.0f',esempio);
    fprintf(fid,'\n A: %1.15e B: %1.15e',A,B);
    if isnan(S) == 0
        fprintf(fid,'\n S: %1.15e',S);
    else
        fprintf(fid,'\n L''equazione Ax^2+By^2=1 non e'' quella di un ellisse');
    end
end

fclose(fid);