Newton's sqrt in c++

#include<bits/stdc++.h> 
#define EPSILON 0.001 
using namespace std; 

double func(double x) { 
    return x*x*x - x*x + 2; 
} 
  
double derivFunc(double x) { 
    return 3*x*x - 2*x; 
} 
   
void sqrt(double x) { 
    double h = func(x) / derivFunc(x); 
    while (abs(h) >= EPSILON) { 
        h = func(x)/derivFunc(x);    
        x = x - h; 
    } 
    cout << "The value of the root is : " << x; 
} 
  
int main() { 
  	ios::sync_with_stdio(0);
  	cin.tie(0);
    double x0;
  	cin >> x0;
    sqrt(x0); 
    return 0; 
}