вторая производная в С

#include <iostream.h>
#include <math.h>
 
float f( float x)
    {
      return((cos(x)+6)/sin(x));
    }
 
void main()
    {
      float eps = 0.001;
      float d1,d2;//d1 - первая производная, d2 - вторая производная
      //по классическому определению производной (f'(x) = (f(x) - f(x0))/(x - x0)) определяем первую производную
      //с точностью eps, т.е. eps = x - x0
      //для удобства введем переменную x0 = x + eps
      float x;//аргумент функции
      cin >> x;
      float x0 = x + eps;
      //находим первую производную
      d1 = (f(x) - f(x0))/(x - x0);
      //таким же образом опять преобразовываем полученное выражение, т.к. f''(x) = (f'(x))' 
      //также определяем производную с точностью eps
      float Numerator = (f(x + eps) - f(x0 + eps))/((x + eps) - (x0 + eps)) - (f(x) - f(x0))/(x - x0);
      //считаем вторую производную, поделив числитель дроби на (x+eps) - x, т.е. на eps
      d2 = Numerator/eps;
      //Выводим результат
      cout << "f''(x) = " << d2 << endl;
    }