#include #include #include using namespace std; double powerlaw(double x) { return 5*pow(x,-1.1); } double bpl(double x) { double x_b=5; double f_b=10; double gamma1=-0.8; double gamma2=-1.5; if(x>x>>x1>>x2>>y>>y1>>y2) { cout<