double fgcd(double a, double b)
{
    if (b > a)
    {
        return fgcd(b, a);
    }
    if (fabs(b) < 0.001)
    {
        return a;
    }
    else
    {
        return (fgcd(b, a - floor(a / b) * b));
    }
}