/* nag_quartic_roots (c02alc) Example Program. * * Copyright 2000 Numerical Algorithms Group. * * NAG C Library * * Mark 6, 2000. * Mark 7, revised, 2001. */ #include #include #include #include int main(void) { double a, b, c, d, e; double *errest=0, *zeroi=0, *zeror=0; Integer i; Integer exit_status=0; NagError fail; INIT_FAIL(fail); Vprintf("nag_quartic_roots (c02alc) Example Program Results\n\n"); if ( ! (errest = NAG_ALLOC(4,double)) || ! (zeroi = NAG_ALLOC(4,double)) || ! (zeror = NAG_ALLOC(4,double)) ) { Vprintf ("Allocation failure\n"); exit_status=-1; goto END; } /* Skip heading in data file */ Vscanf("%*[^\n]"); Vscanf("%lf %lf %lf %lf %lf", &e, &a, &b, &c, &d); /* nag_quartic_roots (c02alc). * Zeros of a real quartic polynomial with real coefficients */ nag_quartic_roots (e, a, b, c, d, zeror, zeroi, errest, &fail); if (fail.code == NE_NOERROR) { Vprintf(" Roots of quartic equation Error estimates\n"); Vprintf(" (machine-dependent)\n\n"); for (i = 0; i <= 3; ++i) { Vprintf("%s %10.5f %10.5f%s %g\n", " z =", zeror[i], zeroi[i], "*i", errest[i]); } } else { Vprintf("Error from nag_quartic_roots (c02alc).\n%s\n", fail.message); exit_status = 1; goto END; } END: if (errest) NAG_FREE (errest); if (zeroi) NAG_FREE (zeroi); if (zeror) NAG_FREE (zeror); return exit_status; }