File: mysubs.f

package info (click to toggle)
fudgit 2.42-6
  • links: PTS
  • area: non-free
  • in suites: potato, woody
  • size: 2,468 kB
  • ctags: 2,375
  • sloc: ansic: 27,729; makefile: 793; yacc: 724; lex: 102; asm: 29; fortran: 15
file content (24 lines) | stat: -rw-r--r-- 687 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c    Everything has to be double precision
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
       double precision function myfunc(x, y)
       double precision x, y

       myfunc = x + y
       return
       end
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     Try with a vector
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
       subroutine mysub(vec, dn)
       double precision dn
       double precision vec(*)
       integer i, top

       top = int(dn)
       do 12 i = 1, top
            vec(i) = dreal(i) + sin(dreal(i))
12     continue
       return
       end