program asymptote real sumnew, sumold, term ,error ,n sumold=0 sumnew=0 n=1 error=1 do while (error .gt. 0.01) term=1/n sumold=sumnew sumnew=sumnew+term error=(sumnew-sumold)/sumnew*100 n=n+1 end do print*,'area of asymptote is ',sumnew, ' after ',(n-1),' terms' end