double precision dp(-2:0),sp(-2:0),fp(-2:0)
      data dp/1d0,0d0,0d0/
      data sp/0d0,1d0,0d0/
      data fp/0d0,0d0,1d0/
      save dp,sp,fp