Mercurial > hg > octave-thorsten

diff examples/fortsub.f @ 6572:8e7148b84b59

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-04-25 04:13:44 by jwe]
author jwe
date Wed, 25 Apr 2007 04:14:49 +0000
parents
children 4270ded9ddc6
line wrap: on
line diff
new file mode 100644
--- /dev/null
+++ b/examples/fortsub.f
@@ -0,0 +1,21 @@
+      subroutine fortsub (n, a, s)
+      implicit none
+      character*(*) s
+      real*8 a(*)
+      integer*4 i, n, ioerr
+      do i = 1, n
+        if (a(i) .eq. 0d0) then
+          call xstopx ('fortsub: divide by zero')
+        else
+          a(i) = 1d0 / a(i)
+        endif
+      enddo
+      write (unit = s, fmt = '(a,i3,a,a)', iostat = ioerr)
+     $       'There are ', n, ' values in the input vector',
+     $       char(0)
+      if (ioerr .ne. 0) then
+        call xstopx ('fortsub: error writing string')
+      endif
+      return
+      end
+