SOR方法用矩阵运算如下
A = [4 3 0;3 4 -1;-1 0 4]
b = [24, 30, -24]
tol = 1e-3
max_iter = 100
x = [1.0, 1.0, 1.0]
err = tol + 1
for k in 1:max_iter
global x #声明为全局变量
x_new = A \ (b - A*x)
err = maximum(abs.(x_new - x)) / maximum(abs.(x_new))
if err <= tol
break
end
x = x_new
end
println("Gauss-Seidel Method Solution: ", x)
x = [1.0, 1.0, 1.0]
w = 1.25
for k in 1:max_iter
global x
x_new = w * A * x + (1 - w) * A \ (b - A * x)
err = maximum(abs.(x_new - x)) / maximum(abs.(x_new))
if err <= tol
break
end
x = x_new
end
println("SOR Method Solution: ", x)
Julia代码运行出来报错：ERROR: DimensionMismatch: dimensions must match: a has dims (Base.OneTo(4), Base.OneTo(1)), b has dims (Base.OneTo(4), Base.OneTo(4)), mismatch at 2。如何修改代码？