Hint for HW 2.6

 

C f' and f" from five-point formulas

 C

       DO      100  I = 3, N-2

         F1(I) = (FI(I-2) - 8.*FI(I-1)+8.*FI(I+1)-FI(I+2))/(12.*H)

         F2(I)=(-FI(I-2)+16.*FI(I-1)-30.*FI(I)+16.*FI(I+1)-FI(I+2))

         F2(I) = F2(I)/(12.*H*H)

   100 CONTINUE

C

 C Linear extrapolation for the boundary points

 C

       F1(2) = 2.*F1(3) - F1(4)

       F1(1) = 2.*F1(2) - F1(3)

       F1(N-1) = 2.*F1(N-2) - F1(N-3)

       F1(N) = 2.*F1(N-1) - F1(N-2)

       F2(2) = 2.*F2(3) - F2(4)

       F2(1) = 2.*F2(2) - F2(3)

       F2(N-1) = 2.*F2(N-2) - F2(N-3)

       F2(N) = 2.*F2(N-1) - F2(N-2)