% Resevanje kvadratne enacbe 
% (zgled 1.7.3 iz knjige Razsirjen uvod v numericne metode)
%
% Eno niclo moramo izracunati preko Vietove formule, da bo stabilno

% Bor Plestenjak 2010

format long e

% Imamo kvadratno enacbo ax^2+bx+c=0, kjer so

a = 1.2345678;
b = 76543210.5;
c = 0.1122334455;

% Klasicen nacin resevanja je preko diskriminante
% in znane formule x_{1,2}=(-b+-sqrt(D))/(2a)
%
% Tako dobimo

disp('Klasicen nacin preko formule x_{1,2}=(-b+-sqrt(D))/(2a)')

D = b^2-4*a*c;
x1 = (-b+sqrt(D))/(2*a)
x2 = (-b-sqrt(D))/(2*a)

% Vendar, tocni nicli sta

disp('Tocni nicli za primerjavo')
tocno = roots([a b c])

% in ena nicla je tocna, druga pa nima niti ene tocne decimalke!
pause

disp('Drugo niclo izracunano po Vietovi formuli')
x1t = c/a/x2

