Previous: utde3 Up: ../plot79_u.html Next: utdpcr

UTDPCP

       SUBROUTINE  UTDPCP (C, N, X, Y, WORKA, WORKB)
 C$    (Polynomial - Coefficients from Points)
 C$    Evaluate the coefficients of a polynomial of degree N which
 C$    passes exactly through N+1  given points (X(K),Y(K)).   The
 C$    arguments (except N) are  all DOUBLE PRECISION arrays;  use
 C$    of maximal precision is desirable for numerical  stability.
 C$
 C$    The output argument is:
 C$
 C$    C(*)......Coefficients such that the polynomial is
 C$              P(x) = C(1) + C(2)*x + C(3)*x**2 + ... +
 C$              C(N+1)*x**N
 C$
 C$    The input arguments are:
 C$
 C$    N.........Order of the polynomial.
 C$    X(*)......N+1 x values.
 C$    Y(*)......N+1 y values, such that y(K) = P(x(K)).
 C$
 C$    The scratch arguments are:
 C$
 C$    WORKA(*).....Working array of at least N+2 locations.
 C$    WORKB(*).....Working array of at least N+1 locations.
 C$
 C$    Adapted from Harwell Subroutine Library routine PC02AD.
 C$    (27-JAN-83)