﻿function [roots] = DannyAllPoly(c, x0, TOL, N)
roots=[];
n=length(c);
k=1;
while (k<=(n - 3))
    [x,q] = DannyNewtonHorner(c, x0, TOL, N);
    % disp(x)
    % disp(roots)
    roots = [roots, x];
    c=q;
    k++; 
endwhile

a=c(1);
b=c(2);
d=c(3);

[r1, r2] = DannyQuadraticRoots(a,b,d);
roots=[roots, r1, r2];
end