Уточнить корни уравнения в заданном диапазоне с заданной погрешностью методом Ньютона - C (СИ)

#include <stdio.h>
#include <math.h>
#define EPS 0.001
 
double f(double x)
{
 return (acos(x)-sqrt(1-0.3*x*x*x));
}
 
double f_prim(double x)
{
 return ((0.45*x*x/sqrt(1-0.3*x*x*x))-(1/sqrt(1-x*x)));
}
 
int main()
{
    int i=0;
    double mas[1000]={0.01,0};
 
    do
    {
        i++;
        mas[i]=mas[i-1]-f(mas[i-1])/f_prim(mas[i-1]);
        printf("xk= %lf, x=%lf, |xk-xk-1|=%lf\n",mas[i],mas[i-1],mas[i]-mas[i-1]);
    }
    while((mas[i]-mas[i-1]>EPS)||(mas[i]-mas[i-1]<-EPS));
 
printf("\nx=%lf\n",mas[i]);
return 0;
}