subroutine Start_Values()
use work
! computes x(-1), x(-2), x(-3), x(-4), ...
! Via BDF of order 1, 2, 3, 4


xa(0,:) =x0(:)



!  0.000000000000000E+000  0.000000000000000E+000
!  1.000000000000000E-002 -1.000008351517986E-002
!  2.000000000000000E-002 -2.000266696365973E-002
!  3.000000000000000E-002 -3.001575475715295E-002

!  4.000000000000000E-002 -4.005337915863078E-002

! x(-1):
xa(-1,2) = -3.001575475715295E-002
xa(-2,2) = -2.000266696365973E-002
xa(-3,2) = -1.000008351517986E-002
xa(-4,2) = 0.

!  0.000000000000000E+000   1.00000000000000     
!  1.000000000000000E-002  0.999949999763629     
!  2.000000000000000E-002  0.999799990502306     
!  3.000000000000000E-002  0.999549911986695     

!  4.000000000000000E-002  0.999199594037768     

xa(-1,1) =0.999549911986695
xa(-2,1) =0.999799990502306 
xa(-3,1) =0.999949999763629
xa(-4,1) =1.00000000000000


END subroutine Start_Values
