const double mu=0.891;
double func1(double x)
{
return mu*x+2*(1-mu)*x*x/(1+x*x);
}