/* erot2.f -- translated by f2c (version of 23 April 1993  18:34:30).
   You must link the resulting object file with the libraries:
	-lf2c -lm   (in that order)
*/

#include "f2c.h"

/* Table of constant values */

static integer c__1 = 1;

/* 		Iwannou Spyros, Kontopodis Dhmhtrhs */
/* Main program */ MAIN__()
{
    /* Format strings */
    static char fmt_11[] = "(3f8.5)";
    static char fmt_20[] = "(\002COND=\002,d20.15)";

    /* System generated locals */
    integer i__1, i__2;

    /* Builtin functions */
    integer s_wsfe(), do_fio(), e_wsfe();
    /* Subroutine */ int s_stop();

    /* Local variables */
    static doublereal cond;
    static integer ndim;
    static doublereal work[3];
    static integer ipvt[3];
    static doublereal h[9]	/* was [3][3] */;
    static integer i, j, n;
    extern /* Subroutine */ int decomp_();

    /* Fortran I/O blocks */
    static cilist io___6 = { 0, 6, 0, fmt_11, 0 };
    static cilist io___10 = { 0, 6, 0, fmt_20, 0 };


    n = 3;
    ndim = 3;
    i__1 = n;
    for (j = 1; j <= i__1; ++j) {
	i__2 = n;
	for (i = 1; i <= i__2; ++i) {
	    h[i + j * 3 - 4] = (float)1. / (real) (i + j - 1);
/* L10: */
	}
    }
/* L11: */
    i__2 = n;
    for (i = 1; i <= i__2; ++i) {
	s_wsfe(&io___6);
	i__1 = n;
	for (j = 1; j <= i__1; ++j) {
	    do_fio(&c__1, (char *)&h[i + j * 3 - 4], (ftnlen)sizeof(
		    doublereal));
	}
	e_wsfe();
/* L15: */
    }
    decomp_(&ndim, &n, h, &cond, ipvt, work);
/* L20: */
    s_wsfe(&io___10);
    do_fio(&c__1, (char *)&cond, (ftnlen)sizeof(doublereal));
    e_wsfe();
    s_stop("", 0L);
} /* MAIN__ */

/* Main program alias */ int erotisi2_ () { MAIN__ (); }