rm(list=ls());
x=matrix(c(1.3615,3.5616,-14.2411,-4.4950,2.3014,1.1066,-9.3409, 0.3779, 0.9386,-0.1838),1);
a=2;
print(a);
n=10;
for (i in 1:7)
{
  a=a+4*sum((x-a)/(1 + (x-a)^2))/n;
  print(a)
}