C语言 计算一元三次方程ax³+bx²+cx+d=0中x在1附近的一个实根

#include <stdio.h>
#include <math.h> 
float solut(float a,float b,float c,float d){
	float x=1,x0,f,f1;
	do{
		x0=x;
		f=((a*x0+b)*x0+c)*x0+d;
		f1=(3*a*x0+2*b)*x0+c;
		x=x0-f/f1;
	}while(fabs(x-x0)>=1e-3);
	return x;
}
	
int main(){
	float a,b,c,d;
	printf("input a,b,c,d:");
	scanf("%f,%f,%f,%f",&a,&b,&c,&d);
	printf("x=%10.7f",solut(a,b,c,d));
	return 0;
}