program demo
c       This program evaluates polynomials. The user inputs the degree
c       and the coefficients of the polynomial, and a point where the
c       polynomial is to be evaluated. A function subprogram evaluates
c       the polynomial at the point.
        real a(0:20), x, poly
        integer n
        character ans*1

        print *, "What is the degree of the polynomial?"
        print *, "(no more than 20, please)"
        read *, n
        do i = 0, n
        print *, "Enter the coefficient of the term of degree", i
        read *, a(i)
        end do
 5      print *, "What is the point x of interest?"
        read *, x
        print *, "The value of the polynomial there is", poly(a,x,n)
        print *
        print *, "Do you want to evaluate at another point (y or n)?"
        read*, ans
        if (ans .eq. "y" .or. ans .eq. "Y") go to 5
        print *, "OK, then - bye"
        end

        function poly(a,x,n)
        real poly, a(0:20), x
        integer n
        poly = a(0)
        xpower = 1
        do i=1, n
                xpower = xpower * x
                poly = poly + a(i) * xpower
        end do
        end