#include <iostream>
#include <iomanip>

using namespace std;

typedef float PopulationT;

const PopulationT R = 3.97;
const int MAX_POPULATION = 10000;

PopulationT NewPopulation(PopulationT oldPopulation);

int main() {
    int i;

    PopulationT currentGen,
           nextGen;

   currentGen = .5;

   cout << setprecision(10);
   cout << "Generation    % of Max    Population" << endl;
   for (i=0; i < 20; i++) {
       cout <<setw(3) << i  << setw(20) << currentGen*100 << "% or " 
            << int(MAX_POPULATION * currentGen)  << endl;
       nextGen = NewPopulation(currentGen);
       currentGen = nextGen; 
   }
   cout <<setw(3) << i << setw(20) << currentGen*100 << "% or " 
        << int(MAX_POPULATION * currentGen)  << endl;


   return 0;
}

PopulationT NewPopulation(PopulationT oldPop){
     PopulationT newPop;

     newPop = R * oldPop * (1.0 -  oldPop);

     return newPop;
}