| Revision: | 1.1 |
| Committed: | Sat Aug 21 02:17:16 2010 UTC (13 years, 1 month ago) by rick |
| Branch: | MAIN |
| CVS Tags: | Ver_6-0, Ver_6-1, Ver_6-2, Ver_6-3, Ver_6-4, Ver_9-1, Ver_5-14, Ver_5-13, Ver_5-12, Ver_5-11, Ver_5-10, Ver_LATEST, Ver_9-3, Ver_9-41, Ver_9-2, Ver_8-8, Ver_8-2, Ver_8-3, Ver_8-0, Ver_8-1, Ver_8-6, Ver_8-7, Ver_8-4, Ver_8-5, Ver_7-1, Ver_7-0, Ver_9-5, Ver_9-4, Ver_8-10, Ver_8-11, Ver_8-12, Ver_9-0, HEAD |
| Log Message: | added for 5.9 |
| # | Content |
|---|---|
| 1 | module polyval_f |
| 2 | integer, parameter :: d=8 |
| 3 | contains |
| 4 | REAL(kind=d) FUNCTION POLYVAL(n,c,x) |
| 5 | ! |
| 6 | ! Evaluate the polynomial |
| 7 | ! |
| 8 | ! c(n+1)*x^n + c(n)*x^(n-1) ... c(2)*x + c(1) |
| 9 | ! |
| 10 | ! using Horner's scheme. |
| 11 | ! |
| 12 | implicit none |
| 13 | integer, parameter :: double=8 |
| 14 | integer :: n,i |
| 15 | real(kind=double), dimension(n+1) :: c |
| 16 | real(kind=double) :: x |
| 17 | |
| 18 | polyval = c(n+1) |
| 19 | do i = n,1,-1 |
| 20 | polyval = (x*polyval + c(i)) |
| 21 | end do |
| 22 | return |
| 23 | end function |
| 24 | end module |