double Lagrange(double *X, double *Y, int NumberOfKnots, double X_input){
	double solution = 0, *temp;
	int i, j;
	
	temp = new double[NumberOfKnots];
	
	for(i = 0;i < NumberOfKnots;i++){
		temp[i] = 1;
		for(j = 0;j < NumberOfKnots;j++){
			if(i != j){
				temp[i] = temp[i] * (X_input - X[j]) / (X[i] - X[j]);
			}
		}
		solution = solution + Y[i] * temp[i];
	}
	
	delete []temp;
	return solution;
}