Numerical integration: Difference between revisions
Content added Content deleted
(CL numerical integration) |
m (→{{header|Fortran}}: Move a comment to avoid a long line) |
||
Line 747: | Line 747: | ||
real, parameter :: h = (b - a) / n |
real, parameter :: h = (b - a) / n |
||
real, parameter, dimension(0:2*n) :: xpoints = (/ (a + h*i/(n*2), i = 0, 2*n) /) |
real, parameter, dimension(0:2*n) :: xpoints = (/ (a + h*i/(n*2), i = 0, 2*n) /) |
||
! gather up all the f(x) values needed for all methods once and distribute via pointers |
|||
real, dimension(0:2*n), target :: fpoints |
|||
real, dimension(:), pointer :: fleft, fmid, fright |
real, dimension(:), pointer :: fleft, fmid, fright |
||
real :: leftrect, midrect, rightrect, trapezoid, simpson |
real :: leftrect, midrect, rightrect, trapezoid, simpson |