# causal integral
#
subroutine causint( conj, add,        n,pp,  qq   )
integer i,       n, conj, add;   real pp(n), qq(n )
temporary real tt( n)
call conjnull(      conj, add,        pp,n,  qq,n )
if( conj == 0){				tt(1) = pp(1) / 2.
				do i= 2, n
					tt(i) = tt(i-1) + (pp(i) + pp(i-1)) / 2.
		do i= 1, n
		  qq(i) = qq(i) + tt(i)
		}

else {					tt(n) = qq(n) / 2.
				do i= n, 2, -1 
					tt(i-1) = tt(i) + (qq(i) + qq(i-1)) / 2.
		do i= 1, n
		  pp(i) = pp(i) + tt(i)
		}
return; end