om = logspace(0,5,100)';
for indice = 1:length(om),
    w = om(indice);
    exec('/home/maurizio/Documenti/Scienza/Laboratorio E 21-22/Esercitazione 7/CircEsZNow.sce', -1);
    [E,V,I] = SolReteMetodoFasoriale(Circ1ZNow);
    out(indice) = E(3);
end
clf();
plot(om,abs(out));
qq = gca();
qq.log_flags = 'lnn';
xgrid();
xlabel('pulsazione (rad)');
qq.x_label.font_size = 2;
ylabel('modulo');
qq.y_label.font_size = 2;
qq.title.text = 'Risposta in frequenza';
qq.title.font_size = 2;
qq.title.font_style = 8;
qq.children.children.thickness = 2;
