c I(t,p,q,r,u,v)
	 complex function I(t,p,q,r,u,v)
  
         implicit double precision(a-h,o-z)
	 double precision t,p,q,r,u,v
	 complex A1,A2,A3,A4,A5,A6,A7
	 complex F0,F02,F1,F12,F2,F22,F3,F32,y1,y2
	
	 A1=cdsqrt(q*q-4*r); A2=-2*p+q
	 A3=cdsqrt(t*t+t*u+v); A4=-2*q*t-q*u+2*t*u+4*v
	 A5=-q*q*(2*t+u)+8*r*t+4*r*u
	 A6=(q*q-2*r+2*v-q*u)/2
	 A7=(q-u)*A1/2
	 F0=A1+A2; F02=-A2+A1
	 F1=cdsqrt(A7+A6); F12=cdsqrt(-A7+A6)
	 F2=q+2*t+A1; F22=q+2*t-A1
	 F3=A4*A1+A5; F32=-A4*A1+A5
	 y1=cdlog((4*A1*A3+F3/F1)/F0/F2)
	 y2=cdlog((4*A1*A3-F32/F12)/F02/F22)
	 I=-(F0*y1/F1+F02*y2/F12)/2/A1

c	 write(*,*) 't',t,p,q,r,u,v
c	 pause
c	 write(*,*) 'A',A1,A2,A3,A4,A5,A6,A7
c	 pause
c	 write(*,*) 'F',F0,F02,F1,F12,F2,F22,F3,F32
c	 pause
c	 write(*,*) 'y1,y2',y1,y2
c	 write(*,*) 'Iinside',I
	 if (A1.eq.0)  then
	    I=0
	 endif
	 

c	 pause
	 return
	 end