Итерационный цикл

#include <iostream>
#include <math.h>
 
using namespace std;
 
void main() {
    long double x, r, a, b, p, realP;
    const long double pi = 3.14159;
    cout << "Input accuracy: ";
    cin >> x;
    cout << "Input radius: ";
    cin >> r;
    a = (2 * r)/sqrt(3);
    b = 100;
    int n = 6;
    realP = 2 * pi * r;
    cout << "Real perimeter is " << realP << '\n';
 
    while ((x < b) & (a > 0))
    {
        cout << "n= " << n << '\t';
        cout << "a= " << a <<"\t";
        p = n*a;
        a = sqrt(2 * r*r - r*sqrt(4 * r*r - a*a));
        n*=2;
        b = fabs(p - realP);
        cout << "p= " << p << '\t';
        cout << "b= " << b << '\t' << '\n';
    }
    cout << '\n' << "p= " << p << '\n';
    system("pause");
}

#include <iostream>
#include <math.h>
 
using namespace std;
 
void main() {
    long double x, r, a, b, p, realP;
    const long double pi = 3.14159;
    cout << "Input accuracy: ";
    cin >> x; //Вводим точность (предел)
    cout << "Input radius: ";
    cin >> r; //Вводим радиус
    int n = 6;
    a = (2 * r)/sqrt(3); //По формуле вычисляем сторону правильного шестиугольника
    b = 100; //Вводим произвольное значение >x, чтобы войти в цикл
    realP = 2 * pi * r; //Высчитываем действительный периметр окружности
    cout << "Real perimeter is " << realP << '\n';
 
    while ((a > x) && (x < b))
    {
        cout << "n= " << n << '\t';  
        cout << "a= " << a <<"\t";
        p = n*a; //Cчитаем периметр. В первой итерации всегда правильно посчитается периметр шестиугольника
        n*=2; //Удваиваем число сторон
        b = fabs(p - realP); //Находим абсолютную величину разницы между фактически посчитаным программой периметром и реальным значением периметра
        a = sqrt(2 * r*r - r*sqrt(4 * r*r - a*a)); // Применяем формулу удвоения стороны n-угольника
        cout << "p= " << p << '\t';
        cout << "b= " << b << '\t' << '\n';
    }
    cout << '\n' << "p= " << p << '\n'; //Выводим окончательный результат работы программы
    system("pause");
}

#include <iostream>
#define _USE_MATH_DEFINES 1
#include <math.h>
 
using namespace std;
 
void main() {
    double x, r, a, b, p, realP;
    const double pi = M_PI;
    cout << "Input accuracy: ";
    cin >> x;
    cout << "Input radius: ";
    cin >> r;
    a = r;          // а никак не  a = (2 * r)/sqrt(3); - сторона вписанного 6-угольника равна радиусу описанной окружности
    b = 100;
    int n = 6;
    realP = 2 * pi * r;
    cout << "Real perimeter is " << realP << '\n';
    p = n*a;
 
    double step_p = x;
 
    while (step_p >= x)
    {
        a = sqrt(2 * r*r - r*sqrt(4 * r*r - a*a));
        n*=2;
        cout << "n= " << n << '\t';
        cout << "a= " << a <<"\t";
        step_p = -p;
        p = n*a;
        step_p += p;
        cout << "step_p= " << step_p <<"\t";
 
        b = fabs(p - realP);
        cout << "p= " << p << '\t';
        cout << "b= " << b << '\n';
    }
    cout << '\n' << "p= " << p << '\n';
    cout << '\n' << "pi= " << (p/(2*r)) << '\n';
 
    system("pause");
}