Wilkinson2

subroutine wilkinson(W, n, MDIM)
    integer MDIM
    real W(MDIM,n)
c numero di prove
    do k1 = 1, 10
        s = SECNDS(0.0)
c ciclo perditempo
        do k2 = 1,1000
c svolgimento
            do i = 1, n-1
                W(i,i) = 1
                W(i+1,i) = 1
                do j = i+2, n
                    W(j,i) = -W(j-1,i)
                    W(i,j) = 0
                end do
                W(i,n) = -W(i+1,1)
            end do
            W(n,n) = 1
c fine svolgimento
        end do
c fine ciclo perditempo
        s = SECNDS(s)
        write(*,*) s
c fine ciclo numero di prove
    end do
    end





    subroutine matrix_write(A,rig, col, MDIM)
    integer MDIM, rig, col
    real A(MDIM,col)
    do i = 1, rig
        write(*,*)(A(i,j), j= 1, col)
        write(*,*)
    end do
      end



    program main
    parameter(MDIM = 1000)
    real W(MDIM,MDIM)
    write(*,*)'Inserire l''ordine della matrice:'
    read(*,*) n
    call wilkinson(W, n, MDIM)
    write(*,*)'Matrice di Wilkinson:'
    write(*,*)
c   call matrix_write(W, n, n, MDIM)
    end