In [1]:
## Intro to stability

function hilb(n::Integer)
  A = zeros(n,n)
  for j=1:n
    for i=1:n
      A[i,j] = 1.0/(i+j-1)
    end
  end
  return A
end
display(hilb(5))
5×5 Array{Float64,2}:
 1.0       0.5       0.333333  0.25      0.2     
 0.5       0.333333  0.25      0.2       0.166667
 0.333333  0.25      0.2       0.166667  0.142857
 0.25      0.2       0.166667  0.142857  0.125   
 0.2       0.166667  0.142857  0.125     0.111111