double precision s1, s2, udummy 
      s1(udummy) = 0.75 - udummy**2
      s2(udummy) = 0.5*(udummy-1.5)**2