% confronto tra due formule per il calcolo delle radici
close all
clear all
clc

format long e

% radici esatte
% metto x1 la radice piu' grande
x1e = 10^5;
x2e = -1.0e-2;

% calcolo s e p
s = (x1e+x2e)/2;
p = x1e*x2e;

% calcolo le radici con la prima formula
x1 = s + sqrt(s*s-p);
x2 = s - sqrt(s*s-p);

% calcolo le radici con la seconda formula
x2n = p/(s + sqrt(s*s-p));
x1n = p/x2n;

% confronto risultati
x1e
x1n
x1

x2e
x2n
x2
