************************************************************************** 
*     simple program for testing the subroutine SUNRISE 
*     g77 examplesunem.for sunem.for -o sunem; ./sunem
************************************************************************** 
      program examplesunem
      implicit none
      real*8 S,SR2,SI2,S1R2,S1I2,SR4,SI4,S1R4,S1I4

      write(6,*)'Input squared momentum transfer S:'
      read(5,*) S

      call sunrise(S,SR2,SI2,S1R2,S1I2,SR4,SI4,S1R4,S1I4) 

      write(6,*)'First Master integral in d=2 dimensions'
      write(6,200) S,SR2
      write(6,201) S,SI2
      write(6,*)'Second Master integral in d=2 dimensions'
      write(6,210) S,S1R2
      write(6,211) S,S1I2

      write(6,*)'First Master integral in d=4 dimensions: finite parts'
      write(6,400) S,SR4
      write(6,401) S,SI4
      write(6,*)'Second Master integral in d=4 dimensions: finite parts'
      write(6,410) S,S1R4
      write(6,411) S,S1I4

 200  format('  Re S(2,',1PE22.15,')     = ',1PE22.15)
 201  format('  Im S(2,',1PE22.15,')     = ',1PE22.15)
 210  format('  Re S1(2,',1PE22.15,')    = ',1PE22.15)
 211  format('  Im S1(2,',1PE22.15,')    = ',1PE22.15)

 400  format('  Re S^(0)(4,',1PE22.15,') = ',1PE22.15)
 401  format('  Im S^(0)(4,',1PE22.15,') = ',1PE22.15)
 410  format('  Re S1^(0)(4,',1PE22.15,')= ',1PE22.15)
 411  format('  Im S1^(0)(4,',1PE22.15,')= ',1PE22.15)
 
      end
