matlab 脚本转化成vb

function [n,x]=sor22(A,b,X,nm,w,ww)
%用超松弛迭代法求解方程组Ax=b
%输入：A为方程组的系数矩阵，b为方程组右端的列向量，X为迭代初值构成的列向量，nm为最大迭代次数，w为误差精度，ww为松弛因子
%输出：x为求得的方程组的解构成的列向量，n为迭代次数
n=1;
m=length(A);
D=diag(diag(A));        %令A=D-L-U,计算矩阵D
L=tril(-A)+D;           %令A=D-L-U,计算矩阵L
U=triu(-A)+D;          %令A=D-L-U,计算矩阵U
M=inv(D-ww*L)*((1-ww)*D+ww*U);    %计算迭代矩阵
g=ww*inv(D-ww*L)*b;                %计算迭代格式中的常数项
%下面是迭代过程
while n<=nm
    x=M*X+g;         %用迭代格式进行迭代
    if norm(x-X,'inf')<w
      disp('迭代次数为');n
      disp('方程组的解为');x
      return;
      %上面：达到精度要求就结束程序，输出迭代次数和方程组的解
    end
    X=x;n=n+1;
end
%下面：如果达到最大迭代次数仍不收敛，输出警告语句及迭代的最终结果（并不是方程组的解）
disp('在最大迭代次数内不收敛!');
disp('最大迭代次数后的结果为');x
上面是完整的超松弛迭代法
 
 
example：
a=[5 2 1;-1 4 2;2 -3 10];
b=[-12; 20; 3;];易错，勿遗漏中间分号，否组矩阵维数错误
c=1000;
d=5e-6;
f=1.03;
k=[0 ;0; 0;];
g=sor22(a,b,k,c,d,f)