blob: eb4c6330f04315b62c877ffdda4da460fdaa5bb2 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
|
//Example 2.13
clc
clear
a = [-6 11 -6 1];
maxit = 3;
for iter = 1:maxit
a = [a(4)^2 -(a(3)^2 -2*a(2)*a(4)) (a(2)^2 - 2*a(1)*a(3)) -a(1)^2];
root = abs([a(4)/a(3) a(3)/a(2) a(2)/a(1)])^(1/(2^iter));
end
root = round(root*10^5) / 10^5;
disp(root,"Estimated roots for the polynomial are: ")
|
