1) f90 -O dr_dmv10.f
2) a.out
 XEND= 10.0  H= 0.01  NSTEP= 1000
 --- INITIAL CONDITION
   AM     1.8 0.4 -0.9
   QQ     1.0 0.0E+0 0.0E+0 0.0E+0
 --- SOLUTION AT ENDPOINT
   AM     1.7659625085182097 -0.6954641970330873 0.779554981481872
   QQ     -0.7360613641607925 -0.5707869512403533 -0.29938857183078976
 -0.20683908610467177
 ERR HAM= 1.66273651113908E-15  TIME= 1.199677E-3