программирование циклических процессов

#include <iostream>
#include <math.h>
using namespace std;
int main()
{
    setlocale(LC_ALL, "RUS");
    double eps, yn = 1, yn1, e, p;
    float x = 0, xEnd = 1, step = 0.1;
    cout << "Введите значение точности и степеня" << endl;
    cin >> eps >> p;
    for (x; x <= xEnd; x = x + step)
    {
        e = (1 - x);
        yn = 1;
        yn1 = ((p - 1) * yn + e / pow(yn, p - 1)) / p;
        do
        {
            yn = yn1;
            yn1 = ((p - 1) * yn + e / pow(yn, p - 1)) / p;
        } while (fabs(yn1 - yn) > eps);
        cout << x << " " << 1 - x << " " << sqrt(x)/ yn1 << endl;
    }
    return 0;
}

#include <iostream>
#include <math.h>
using namespace std;
int main()
{
    setlocale(LC_ALL, "RUS");
    double eps, yn = 1, yn1, e, p;
    float x = 0, xEnd = 1, step = 0.1;
    cout << "Введите значение точности и степеня" << endl;
    cin >> eps >> p;
    for (x; x <= xEnd; x = x + step)
    {
        e = (1 - x);
        yn = 1;
        yn1 = ((p - 1) * yn + e / pow(yn, p - 1)) / p;
        do
        {
            yn = yn1;
            yn1 = ((p - 1) * yn + e / pow(yn, p - 1)) / p;
        } while (fabs(yn1 - yn) > eps);
        cout << x << " " << 1 - x << " " << sqrt(x)/ yn1 << endl;
    }
    return 0;
}

cout << x << " " << yn1<< " " << sqrt(x)/ (1.0-x) << endl;

#include <iostream>
#include <math.h>
using namespace std;
int main()
{
setlocale(LC_ALL, "RUS");
double eps, yn = 1, yn1, e, p;
float x = 0, xEnd = 1, step = 0.1;
cout << "Введите значение точности и степеня" << endl;
cin >> eps >> p;
for (x; x <= xEnd; x = x + step)
{
e = (1 - x);
yn = 1;
yn1 = ((p - 1) * yn + e / pow(yn, p - 1)) / p;
do
{
yn = yn1;
yn1 = ((p - 1) * yn + e / pow(yn, p - 1)) / p;
} while (fabs(yn1 - yn) > eps);
cout << x << " " << 1 - x << " " << sqrt(x)/ yn1 << endl;
}
return 0;
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{        
    double x1=0., x2=0.9, dx=0.1, yn, prev=0., eps=0.0001;      
    
    cout << " x   result     y(x)\n";
    for (double x=x1; x<=x2+dx/2.; x+=dx)
    {
        yn=1;
        while(fabs(yn-prev)>eps)
        {          
           prev=yn;           
           yn=0.5*(yn+x/yn);                      
        } 
        
    cout << setprecision(1) << fixed<< x << setprecision(5) << "  " << yn/(1-x) << "  " << sqrt(x)/(1-x) << "\n";       
    }
    
system("pause");
return 0;
}