/* quadr.f -- translated by f2c (version 19950314).
   You must link the resulting object file with the libraries:
	-lF77 -lI77 -lm   (in that order)
*/

#include "f2c.h"

/* Subroutine */ int quadr_(doublereal *f0, doublereal *f1, doublereal *f2, 
	doublereal *x1, doublereal *x2, doublereal *a, doublereal *b, 
	doublereal *c)
{
    /* System generated locals */
    doublereal d__1, d__2;

/* *************************************************** */
/*                                                  * */
/*    QUADR CALCULATES THE A, B AND C IN THE EQUNS. * */
/*                                                  * */
/*     A                   =   F0                   * */
/*     A + B.X0 + C.X0**2  =   F1                   * */
/*     A + B.X2 + C.X2**2  =   F2                   * */
/*                                                  * */
/* *************************************************** */
/* Computing 2nd power */
    d__1 = *x1;
/* Computing 2nd power */
    d__2 = *x2;
    *c = (*x2 * (*f1 - *f0) - *x1 * (*f2 - *f0)) / (*x2 * (d__1 * d__1) - *x1 
	    * (d__2 * d__2));
/* Computing 2nd power */
    d__1 = *x1;
    *b = (*f1 - *f0 - *c * (d__1 * d__1)) / *x1;
    *a = *f0;
    return 0;
} /* quadr_ */