C
C    ** HASKELL'S SUB MATRIX FOR R WAVE ***
C
      SUBROUTINE FRSB(DNS,CB2,CV,CQ,STQ,SDQ,CP,STP,SDP,SB)
C
      IMPLICIT REAL*8(A-H,O-Z)
      COMPLEX*16 SB(4,4)
C
      GM=2.D0/CB2
C
      C11=GM
      C12=GM-1.D0
      C13=GM*GM
      C14=C12*C12
      C16=DNS*CV*CV
      C15=1.D0/C16
C
      D11=CP-CQ
      D12=SDP+STQ
      D13=STP+SDQ
      D14=C14*SDP+C13*STQ
      D15=C13*STP+C14*SDQ
C
      SB(1,1)= DCMPLX(C11*CP-C12*CQ,0.0D0)
      SB(1,2)= DCMPLX(0.0D0,C12*SDP+C11*STQ)
      SB(1,3)=-DCMPLX(C15*D11,0.D0)
      SB(1,4)= DCMPLX(0.D0,C15*D12)
      SB(2,1)=-DCMPLX(0.0D0,C11*STP+C12*SDQ)
      SB(2,2)= DCMPLX(-C12*CP+C11*CQ,0.D0)
      SB(2,3)= DCMPLX(0.D0,C15*D13)
      SB(2,4)= SB(1,3)
      SB(3,1)= DCMPLX(C16*C11*C12*D11,0.0D0)
      SB(3,2)= DCMPLX(0.D0,C16*D14)
      SB(3,3)= SB(2,2)
      SB(3,4)= SB(1,2)
      SB(4,1)= DCMPLX(0.D0,C16*D15)
      SB(4,2)= SB(3,1)
      SB(4,3)= SB(2,1)
      SB(4,4)= SB(1,1)
C
      RETURN
      END
