vendor/gonum.org/v1/gonum/lapack/internal/testdata/dsterftest/testdsterf.f90

   1 program testdsterf
   2 implicit none
   3 integer, parameter :: n = 4
   4 real(kind=8), dimension(n) :: d
   5 real(kind=8), dimension(n-1) :: e
   6 integer :: info,i
   7
   8 d(1:4) = (/1D+00, 3D+00, 4D+00, 6D+00/)
   9 e(1:3) = (/2D+00, 4D+00, 5D+00/)
  10
  11 call dsterf(n,d,e,info)
  12 DO i = 1, n
  13     print *, d(i)
  14 end do
  15 end