#include <iostream>
#include <math.h>
using namespace std;
 
int main() {
 
   float a1, a2, h, h1, b, b2;
   a1 = 1;
   a2 = 5;
   h = 1;
 
   h1 = 0.05;
 
     float eps = 1.0e-7;
     //cout << "\nEnter step b:";
     //cin >> h1;
 
    do {
        cout << " " << a1 << endl;;
 
        b = max(a1-1.0,(-(5.0+4.0*a1)+sqrt(25.0+48.0*a1))/4.0);
        b2 = min(a1-0.5,(-1.0+sqrt(1.0+4.0*a1*(1.0-a1)))/2.0);
 
	cout << "\t\tmax(a1-1,(-(5+4*a1)+sqrt(25+48*a1))/4)=" << b << endl;
	cout << "\t\tmin(a1-0.5,(-1+sqrt(1+4*a1*(1-a1)))/2)=" << b2 << endl;
 
        do{
            cout << "\t\t" << b << endl;
            b += h1;
        } while (b <= b2+eps); 
 
        a1 += h;
    } while (a1 <= a2+eps);   
 
	return 0;
}

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8bWF0aC5oPgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKaW50IG1haW4oKSB7CgkKICAgZmxvYXQgYTEsIGEyLCBoLCBoMSwgYiwgYjI7CiAgIGExID0gMTsKICAgYTIgPSA1OwogICBoID0gMTsKICAgCiAgIGgxID0gMC4wNTsKCiAgICAgZmxvYXQgZXBzID0gMS4wZS03OwogICAgIC8vY291dCA8PCAiXG5FbnRlciBzdGVwIGI6IjsKICAgICAvL2NpbiA+PiBoMTsKCiAgICBkbyB7CiAgICAgICAgY291dCA8PCAiICIgPDwgYTEgPDwgZW5kbDs7CgogICAgICAgIGIgPSBtYXgoYTEtMS4wLCgtKDUuMCs0LjAqYTEpK3NxcnQoMjUuMCs0OC4wKmExKSkvNC4wKTsKICAgICAgICBiMiA9IG1pbihhMS0wLjUsKC0xLjArc3FydCgxLjArNC4wKmExKigxLjAtYTEpKSkvMi4wKTsKCgljb3V0IDw8ICJcdFx0bWF4KGExLTEsKC0oNSs0KmExKStzcXJ0KDI1KzQ4KmExKSkvNCk9IiA8PCBiIDw8IGVuZGw7Cgljb3V0IDw8ICJcdFx0bWluKGExLTAuNSwoLTErc3FydCgxKzQqYTEqKDEtYTEpKSkvMik9IiA8PCBiMiA8PCBlbmRsOwoKICAgICAgICBkb3sKICAgICAgICAgICAgY291dCA8PCAiXHRcdCIgPDwgYiA8PCBlbmRsOwogICAgICAgICAgICBiICs9IGgxOwogICAgICAgIH0gd2hpbGUgKGIgPD0gYjIrZXBzKTsgCgogICAgICAgIGExICs9IGg7CiAgICB9IHdoaWxlIChhMSA8PSBhMitlcHMpOyAgIAoJCglyZXR1cm4gMDsKfQ==