#include #include #include #include #include using namespace std; #include #include #include #include #include #include #include #include double zed(double* s, double *p){ double ess = s[0] ; double dv = p[0]; return sin(ess)-dv; // function you wish to invert minus dummy variable } double root_int(double* dummy_variable, double* p) { double dv = dummy_variable[0] ; // dummy variable double B = p[0]; TF1 small_z("zed", zed, 0.0, 0.7, 1); // numerator function small_z.SetParameter(0,dv); // return B*small_z.GetX(0,0,0.7) ; // y of inverted function } int main() { double B = 1 ; // constant TF1* alpha= new TF1("alpha", root_int, 0., 0.7, 1); // numerator function alpha -> SetParameter(0, B); // Set y0, mu, R, E0 alpha-> Draw(); return 0; }