clear;
// Inserire il nome completo del file ReteEs.sce
File = '/home/maurizio/Documents/Scienza/Scilab/Laboratorio E/ReteEs.sce';
// Inserire il nome completo del file Esercitazione_7.sce
File1 = '/home/maurizio/Documents/Scienza/Scilab/Laboratorio E/Esercitazione_7.sce';
exec(File1);
pert = 0;
exec(File);
[E0,V0,I0] = SoluzioneRete(ReteEs);
pert = 1;
ERR_E = [];
N = 100;
for k = 1:N,
    exec(File);
    [E,V,I] = SoluzioneRete(ReteEs);
    ERR_E($+1,1) = norm(E - E0, %inf) / norm(E0, %inf);
end;
clf(); plot2d('nl',[1:N]',ERR_E, style = 5); xgrid();
xlabel('k'); ylabel('ERR_E(k)');
//
// Modificare in modo da calcolare, ad ogni iterazione,
// il massimo del vettore degli errori relativi delle tensioni di nodo
//
//      |E(i) - E0(i)| / |E0(i)|  ,   i = 1,...,N
//
// e poi disegnare tale massimo in funzione di k.
