fnspl1(dspl) = -dspl*dspl + 0.75d0
      fnspl2(dspl) = 0.5d0*( dspl-1.5d0 )*( dspl-1.5d0 )