/*x^2-3 x+1=0の解*/

#include<stdio.h>
#include<math.h>
void main(void){

double x1,x2,cc;

x1=1.50;
cc=1.0;
while(c>=pow(10,-5)){
x2=(x1*x1+1.0)/3.0;
cc=fabs((x2-x1)/x1);

x1=x2;
}
printf("x=%f",x1);
}

で、x1=1.50の辺りから、どうしてそうなるのかわからないので、
教えていただける方、教えてください。