注意，julia中没有diag函数，如果想要定义4*4的对角矩阵，可以使用下面方法：

1. 使用Matrix函数

 using LinearAlgebra
  Matrix{Float64}(I,4,4)

注意，Matrix{Float64}(I,4,4)中是字母I，而不是数字1

julia> using LinearAlgebra

julia> Matrix{Float64}(I,4,4)
4×4 Array{Float64,2}:
 1.0  0.0  0.0  0.0
 0.0  1.0  0.0  0.0
 0.0  0.0  1.0  0.0
 0.0  0.0  0.0  1.0

2. 创建zeros矩阵

b2 = zeros(4,4)
b2

通过for循环实现对角线为1

搞定。