/* benchmrk.f -- translated by f2c (version of 16 May 1991  13:06:06).
   You must link the resulting object file with the libraries:
	-link <S|C|M|L>f2c.lib   (in that order)
*/

#include <stdlib.h>
#include <math.h>
#include <stdio.h>
#include "/home/software/f2c/f2c.h"

/* Table of constant values */

static integer c__9 = 9;
static integer c__1 = 1;
static integer c__3 = 3;
static doublereal c_b6 = 2.2;
static integer c__5 = 5;

/* Das Ergebnis sollte  476.33479743066 sein. */
/* Die Anzahl der outer loops kannst Du zu Anfang */
/* auf 10 setzen, um in den Bereich von Sekunden */
/* zu kommen. Die HP-730 benoetigt ca. 0.02 sec. */
/* Mein alter ATARI-ST (jedoch mit 68881 Co-Proz) */
/* benoetigt 0.7 bis 1.0 sec. mit gutem Compiler */
/* (in Pascal sogar etwas schneller als in Fortran). */
/* Mit einem schlechten Compiler habe ich auch */
/* schon mal 20 sec. benoetigt! */
/* Besser als mit der Stoppuhr geht es natuerlich */
/* mit entspr. Systemroutinen, die aber nicht */
/* auf allen Rechnern verfuegbar sind. */

/* Viel Spass, */
/*  Norbert */

/* Main program */ main()
{
    /* System generated locals */
    integer i__1;

    /* Builtin functions */
    integer s_wsle(), do_lio(), e_wsle(), s_rsle(), e_rsle();
    double pow_dd(), sqrt(), exp();

    /* Local variables */
    static doublereal a;
    static integer i, j, n;
    static doublereal z;

    /* Fortran I/O blocks */
    static cilist io___1 = { 0, 6, 0, 0, 0 };
    static cilist io___2 = { 0, 5, 0, 0, 0 };
    static cilist io___8 = { 0, 6, 0, 0, 0 };


/*  benchmark test for all computers, */
/*  adopted for  HP-workstation */
/*  930526 NZ */
    s_wsle(&io___1);
    do_lio(&c__9, &c__1, "input number of outer loops:", 28L);
    e_wsle();
    s_rsle(&io___2);
    do_lio(&c__3, &c__1, (char *)&n, (ftnlen)sizeof(integer));
    e_rsle();
    i__1 = n;
    for (j = 1; j <= i__1; ++j) {
	a = (float)0.;
	z = (float)0.;
	for (i = 1; i <= 2000; ++i) {
	    z += (float)1.;
	    a += pow_dd(&z, &c_b6) / (exp(sqrt(z)) + (float)1.);
	}
    }
    s_wsle(&io___8);
    do_lio(&c__9, &c__1, "result =  ", 10L);
    do_lio(&c__5, &c__1, (char *)&a, (ftnlen)sizeof(doublereal));
    e_wsle();
} /* MAIN__ */

/* Main program alias */ int savage_ () { main(); }