Hi mikrom,
Thanks a lot. you are absolutely right.
can you look into my integral program. is it like
do i=1,m
sum=0. !phi(x)
do j=1,n
dy=0.01 !intervals of width, as the interval are very small
sum=sum+exp(-F(i,j)*dy)
enddo
write(*,*)x(i),log(sum)
enddo