c++计算导数

#include<iostream>
#include<cmath>
using namespace std;

int main(void){
	double epsilon = 0.0001;  // 10^-4
	// y = x^2+x+3; x=2  d=5
	double x1 = 2+epsilon, x2 = 2-epsilon;
	double f1 = pow(x1,2)+x1+3, f2 = pow(x2,2)+x2+3;

	double derivative = (f1-f2)/(2*epsilon);
	cout<<derivative<<endl;
	return 0;
}