* testing ob fortran
#+name: ftest
#+begin_src fortran :var m=20 :var n=30 :results output
  real, dimension(m,n) :: a
  real, dimension(n,m) :: c
  real, dimension(m) :: b
  integer :: i
  print *, m, n
  a = reshape( (/ (real(i), i=1,m*n) /), (/ m, n /) )
  b = (/ (2*i, i=1,m) /)
  c = reshape( a, (/ n, m /) )
  print *, b
#+end_src

#+results: ftest
: 20          30
: 2.00000000       4.00000000       6.00000000       8.00000000       10.0000000       12.0000000       14.0000000       16.0000000       18.0000000       20.0000000       22.0000000       24.0000000       26.0000000       28.0000000       30.0000000       32.0000000       34.0000000       36.0000000       38.0000000       40.0000000

#+results:
: 2.00000000       4.00000000       6.00000000       8.00000000       10.0000000       12.0000000       14.0000000       16.0000000       18.0000000       20.0000000       22.0000000       24.0000000       26.0000000       28.0000000       30.0000000       32.0000000       34.0000000       36.0000000       38.0000000       40.0000000

#+call: ftest(m=10,n=40)

#+results:
: 10          40
: 2.00000000       4.00000000       6.00000000       8.00000000       10.0000000       12.0000000       14.0000000       16.0000000       18.0000000       20.0000000

#+call: ftest(m=20,n=10)

#+results:
: 20          10
: 2.00000000       4.00000000       6.00000000       8.00000000       10.0000000       12.0000000       14.0000000       16.0000000       18.0000000       20.0000000       22.0000000       24.0000000       26.0000000       28.0000000       30.0000000       32.0000000       34.0000000       36.0000000       38.0000000       40.0000000