c pt2pt/2-add_seq real*8 a,b,c n=3 do i=1,n a=dble(i) b=dble(i*i) call add(a,b,c) write(*,*) 'a=',a,' b=',b,' c=',c enddo stop end c--------------------------------------------- subroutine add(a,b,c) real*8 a,b,c c=a+b return end