subroutine soptions(default)
c SOPTIONS returns the default values for the optional parameters
c used by SolvOpt.
      double precision default(13) 
       default(1)=-1.d0
       default(2)=1.d-4
       default(3)=1.d-6
       default(4)=15.d3
       default(5)=0.d0
       default(6)=1.d-8
       default(7)=2.5d0
       default(8)=1.d-11
       default(9)=0.d0
       default(10)=0.d0
       default(11)=0.d0
       default(12)=0.d0
       default(13)=0.d0
      end