Вычислить значение функции

#include <iostream>
#include <cmath>
 
using namespace std;
 
double fM(double a, double x)
{
   return ((((-1.)*a)*pow(x,2))/2.);
}
 
main()
{
   double x, a, b;
   cout<<"X = ";
   cin>>x;
   cout<<"A = ";
   cin>>a;
   cout<<"B = ";
   cin>>b;
 
 cout<<sqrt( (exp((-1.)*x) + exp((1.)*b)) / (exp(x)+exp((-1.)*a)) + log(fabs(x)) + pow(fM(a,x),2 )/2.)<<endl;
 system("pause");
}

double fM(double a, double x)
{
   return -a*x*x/2.;
}

#include <iostream>
#include <cmath> 
using namespace std;
 
int main()
{
    double a, b, x, W, eps=1e-9;
    cout<<"a="; cin>>a;
    cout<<"b="; cin>>b;    
    do
    {
    cout<<"x="; cin>>x;
    if ((exp(-x)+exp(-b))/(exp(x)+exp(-a))+log(fabs(x)+a*x*x/2.0*a*x*x/2.0/2.0) <-eps || x==0.0)
    cout <<"x does not belong to the domain !   Repeat please...\n";
    }
    while ((exp(-x)+exp(-b))/(exp(x)+exp(-a))+log(fabs(x)+a*x*x/2.0*a*x*x/2.0/2.0) <-eps || x==0.0);    
       
    W=sqrt((exp(-x)+exp(-b))/(exp(x)+exp(-a))+log(fabs(x)+a*x*x/2.0*a*x*x/2.0/2.0));
 
    cout<<"W="<<W<<"\n";
    
system("pause");
return 0;
}