#include <cfloat>
#include <cmath>
#include <iomanip>
#include <iostream>

using namespace std;

/* Konstante kann nicht mit 4*atan(1.0) initialisiert werden
   wegen dem Funktionsaufruf
*/ 
const double pi = 4*atan(1);
const int N = 100;

/* Deklaration, Prototyp der Funktion f */
double f(double x);

int main(void)
{
    int i;

    // pi = 4*atan(1.0);
    cout << setiosflags(ios::scientific);
    for (i = 0; i <= N; i++)
    {
        cout << "i = " << setw(3) << i
             << ", f(" << setw(9) << setprecision(6) 
             << resetiosflags(ios::scientific) << setiosflags(ios::fixed) 
             << i*pi/N << ") = "
             << setw(18) << setprecision(10) 
             << resetiosflags(ios::fixed) << setiosflags(ios::scientific) 
             << f(i*pi/N) << endl;
    }
    return 0;
}

double f(double x)
{
    /* Ausnahmesituation: x nahe Pol pi/2 */
    if (fabs(x - pi/2) <= 100*DBL_EPSILON)
      return DBL_MAX;
    return tan(x);
}