1 c
   2 c Demo code of interface with freePascal mtx09.pas unit
   3 c
   4       subroutine sub1(n,x)
   5 c
   6       integer n
   7       real*8 x(n)
   8 c
   9       integer i
  10 c
  11       do 10 i = 1, n
  12       x(i) = x(i)*2.0
  13    10 continue
  14 c
  15       return
  16       end
  17 c
  18 c compile with: gfortran -c sub1.f
  19 c to get sub1.o
  20 c