Составить программу, вычисляющую значения функции f от двух переменных а и b на C++

 #include <iostream>
#include <cmath>

using namespace std;

double g(double p, double q) {
  return sqrt(p + q + 1) * (exp(p) + exp(-q)) /
         (pow(p + pow(tan(q / M_PI), 2), 2) * sqrt(3 * p + 2 * q));
}

double f(double a, double b) {
  return g(a * b + 1, a + 2 * b) + g(2 * a + b, sqrt(a) + b);
}

int main() {
  double a, b;
  cin >> a >> b;
  cout << f(a, b);
}