C
C *** THE HASKELL'S SUB MATRIX ***
C
      SUBROUTINE ARALSB(SBR,SBL,CV2,DNS,BM2
     *                 ,CQ,STQ,SDQ,CP,STP,SDP)
C
      IMPLICIT REAL*8(A-H,O-Z)
      COMPLEX*16 SBR(4,4),SBL(2,2),GM,SM,BM2,CPX
     *          ,CQ,STQ,SDQ,CP,STP,SDP
     *          ,C11,C12,C13,C14,C15,C16
     *          ,D11,D12,D13,D14,D15
C
      CPX=DCMPLX(0.D0,1.D0)
      GM=2.D0*BM2/CV2
      SM=DNS*BM2
C
      C11=GM
      C12=GM-1.D0
      C13=GM*GM
      C14=C12*C12
      C16=DNS*CV2
      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
      SBR(1,1)= C11*CP-C12*CQ
      SBR(1,2)= (C12*SDP+C11*STQ)*CPX
      SBR(1,3)=-C15*D11
      SBR(1,4)= C15*D12*CPX
      SBR(2,1)=-(C11*STP+C12*SDQ)*CPX
      SBR(2,2)=-C12*CP+C11*CQ
      SBR(2,3)= C15*D13*CPX
      SBR(2,4)= SBR(1,3)
      SBR(3,1)= C16*C11*C12*D11
      SBR(3,2)= C16*D14*CPX
      SBR(3,3)= SBR(2,2)
      SBR(3,4)= SBR(1,2)
      SBR(4,1)= C16*D15*CPX
      SBR(4,2)= SBR(3,1)
      SBR(4,3)= SBR(2,1)
      SBR(4,4)= SBR(1,1)
C
      SBL(1,1)= CQ
      SBL(1,2)= CPX*SDQ/SM
      SBL(2,1)= CPX*STQ*SM
      SBL(2,2)= CQ
C
      RETURN
      END
