Пересечение кругов

#include <math.h>
#include <iostream>
#include <cmath>
using namespace std;
const double pi = 3.1415926535897932384626433832795;
int main() {
    double x0, y0, r0, x1, y1, r1, d, a, h, x2, y2,sc,sc1;
    double x3, y3, kut, s, x3l, y3l, leng,  kut1, s1;
    cin >> x0 >> y0 >> r0;
    cin >> x1 >> y1 >> r1;
    d = sqrt(pow((x0 - x1), 2) + pow((y0 - y1), 2));//відстань між центрами кіл
    
    
        a = (pow(r0, 2) - pow(r1, 2) + pow(d, 2)) / (2 * d);//on picture
        h = sqrt(pow(r0, 2) - pow(a, 2));//on picture
 
        x2 = x0 + (a*(x1 - x0)) / d;
        y2 = y0 + (a*(y1 - y0)) / d;//координати точки п2
 
        x3 = x2 - (h*(y1 - y0)) / d;
        y3 = y2 - (h*(x1 - x0)) / d;
        x3l = x2 + (h*(y1 - y0)) / d;
        y3l = y2 + (h*(x1 - x0)) / d;
        //перший трикутник
        kut = (180 / pi) * (h / r0);//градусна міра кута першого трикутника
 
 
        leng = h * 2;//відстань між точками перетину кіл
        s = (a * leng) / 2;//площа першого трикутника
 
        //другий трикутник
        kut1 = ((180 / pi) * (h / r1));
        s1 = ((d - a) * leng) / 2;//площа другого трикутника
 
        //сигменти
        sc = ((pi*pow(r0, 2)*kut) / 360) - s;
        sc1 = ((pi*pow(r1, 2)*kut) / 360) - s1;
        cout << sc + sc1;
    
    system("pause");
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main() 
{
    double x1, y1, r1, x2, y2, r2, d, a1, a2, h, S1, S2, S;    
    cin >> x1 >> y1 >> r1;
    cin >> x2 >> y2 >> r2;
    d = hypot(x1-x2,y1-y2);   
    a1 = (d*d+r1*r1-r2*r2)/(2*d);    
    a2 = (d*d+r2*r2-r1*r1)/(2*d);    
    h = sqrt(r1*r1-a1*a1);    
    S1 = r1*r1*asin(h/r1)-a1*h;
    S2 = r2*r2*asin(h/r2)-a2*h;
    S=S1+S2;    
    cout <<setprecision(9)<<fixed<<S<<endl;      
system("pause");
return 0;
}