1 SUBROUTINE gety2(y, y2, h, nots, nb)
8 REAL(rprec),
DIMENSION(*) :: y, y2, h
13 REAL(rprec),
DIMENSION(nots) :: aspline, bspline, dspline
17 IF (nb .eq. ideriv) jspmin(1) = 1
23 IF (nb .eq. ideriv) y2(1) = 6.0*(y(2)-y(1))/h(1)
24 aspline(2:jmax) = h(2:jmax)
25 bspline(2:jmax) = h(:jmax-1)
26 dspline(2:jmax) = 2.0*(h(2:jmax)+h(:jmax-1))
27 y2(2:jmax) = 6.0*((y(3:jmax+1)-y(2:jmax))/h(2:jmax)
28 1 -(y(2:jmax)-y(:jmax-1))/h(:jmax-1))
30 CALL tridslv(aspline,dspline,bspline,y2,jspmin,jmax,0,nots,1)