function doit(x)
           doubleprecision x
           doubleprecision doit
           doit = sqrt(x)
           return
       end

       subroutine itimes(x, y, n, w)
           integer x(*)
           integer y(*)
           integer w(*)
           integer n
           integer i
           do 100 i = 1, n
              w(i) = x(i) * y(i)
100        continue
           return
       end

       subroutine times(x, y, n, w)
           doubleprecision x(*)
           doubleprecision y(*)
           doubleprecision w(*)
           integer n
           integer i
           do 100 i = 1, n
              w(i) = x(i) * y(i)
100        continue
           return
       end

       subroutine ctimes(x, y, n, w)
           complex x(*)
           complex y(*)
           complex w(*)
           integer n
           integer i
           do 100 i = 1, n
              w(i) = x(i) * y(i)
100        continue
           return
       end

       function bang( r, ra, n, k, c)
            integer n, k
            real r, ra(n)
            real bang
            complex c(2)
            bang = ra(k) * r * aimag(c(2))
            return
       end

       function two (a, n, m, k, b, c)
         real two
         integer n, m, k
         integer i, j, l
         real a(n, m)
         real b(m, k)
         real c(n, m)
         two = 0.0
         do 50 i = 1, m
            two = two + b(i, 1)
 50      continue
         do 100 i = 1, n
            do 200 j = 1, m
               c(i, j) = a(i, j)
200         continue
100      continue
       end 

       subroutine copy2 (a, n, m,  c)
         integer n, m
         integer i, j
         real a(n, m)
         real c(n, m)
         do 100 i = 1, n
               do 300 j = 1, m
                  c(i, j) = a(i, j) 
300           continue
100      continue
       end 

       subroutine copy1 (a, n,  c) 
         integer n
         integer i
         real a(n)
         real c(n)
         do 100 i = 1, n
              c(i) = a(i) 
100      continue
       end 

       subroutine no_return()
       return
       end

       function sarg (s, x)
       character*(*) s
       integer sarg, x
       if(s.eq."yes") then
          sarg = 1
       else
          sarg = 0
       endif
       x = len(s)
       return
       end

       function ctest (n, c, ca, ct, cout, caout)
          complex ctest
          integer n
          complex c, ca(n), cout, caout(n), ct(3)
          integer i
          cout = c * c
          ct(1) = c
          ct(2) = ca(1)
          ct(3) = cout
          ctest = ct(1) * ct(2) * ct(3)
          do 100 i = 1, n
             caout(i) = ca(i)
100       continue
          return
        end

        subroutine inoutme (n, m)
        integer n, m(n)
        integer i
        do 100 i = 1, n
           m(i) = 2 * m(i)
100        continue
        return
        end 

        function nocheck (n, m, z)
        integer n, m(1)
        integer nocheck
        integer z(3, *)
        nocheck = m(n) + z(3, 3)
        return
        end

        function work (n, x, ly, y)
        real work
        integer n, ly
        real x(n), y(ly)
        integer i
        do 100 i = 1, n
           y(i) = x(i)
           y(n+i) = -x(i)
 100       continue
        work = 0.0
        do 150 i = 1, ly
           work = work + y(i)
 150       continue
        return
        end


       function mirror (a, k, m, n, b, c)
         real mirror
         integer n, m, k
         integer i, j, l
         real a(m, n)
         real b(k, m)
         real c(m, n)
         mirror = 0.0
         do 50 i = 1, m
            mirror = mirror + b(1, i)
 50      continue
         do 100 i = 1, m
            do 200 j = 1, n
               c(i, j) = a(i, j)
200         continue
100      continue
       end 

        function chartest (c)
        character c, chartest
        if (c .eq. "x") then
            chartest = "y"
        else
            chartest = "n"
        endif
        return
        end