Viết hộ mình cái code ra

#include <bits/stdc++.h>
using namespace std;

const double PI = 3.141592653589793238463;

int main(){
 double a,p,s,r1,r2;
 cin >> a;
 p = a*1.5;
 s = sqrt(p*pow(p-a,3));
 r1 = s/p;
 r2 = pow(a,3)/(4*s);
 float k = (r2*r2-r1*r1)*PI;
 cout << fixed << setprecision(6) << k;
}