Программа - решение кубического уравнения методом хорд.

[Гена1]

Нужно написать программу решения кубического уравнения на Си методом хорд! Это цело численный метод решения уравнений( вот описание метода (http://ru.wikipedia.org/wiki/Метод_хорд)) У меня есть написанная мной программа , если кто может , помогите исправить ошибки или напишите заново . Вот мой код: #include <stdio.h>
#include <stdlib.h>
#include <math.h>

int main()
{
double a,b,c,d;//заводим коэффициенты кубического уравнения
double x[10];//заводим массив с ответами(мы получаем всего 1 ответ для этого уравнения , но массив нам нужен так как каждый ответ получается с определенной точностью,п)
double k[10];
double b1[10];
double q[10];
double y1[10];
scanf("%lf",&a);
scanf("%lf",&b);
scanf("%lf",&c);
scanf("%lf",&d);
q[0]=0;
double t=1;
double y2;
k[0]=((y2-y1[0])/(t-q[0]));
b1[0]=y1[0];
y1[0]=a*0+b*0+c*0+d;
y2=a*1+b*1+c*1+d;
int i;
for(i=0;i<10;i++)
{
q[i]=-b1[i]*k[i];
y1[i]=a*(q[i]*q[i]*q[i])+b*(q[i]*q[i])+c*q[i]+d;
k[i]=((y2-y1[i])/(t-q[i]));
b1[i]=y1[i]-(k[i]*q[i]);
x[i+1]=x[i]-(((x[i]-x[i-1])*(k[i]*x[i]+b1[i]))/((k[i]*x[i]+b1[i])-(k[i-1]*x[i-1]+b1[i-1])));
printf("%lf",x[i+1]);
}
}