unit2.pas

来自「《Delphi常用数值算法集》的配书源码」· PAS 代码 · 共 34 行

unit Unit2;

interface
uses
  unit1;
Procedure SVDVAR(V:matrx2; MA, NP:integer; W:array of real;
                  var CVM:matrx2; NCVM:integer);
implementation
Procedure SVDVAR(V:matrx2; MA, NP:integer; W:array of real;
                  var CVM:matrx2; NCVM:integer);
var
    WTI:array[0..20] of real;
    I,K,J:integer;  SUM1:real;
begin
    For I:=1 To MA do
    begin
        WTI[I]:=0 ;
        If W[I] <> 0  Then WTI[I]:=1 / (W[I] * W[I]);
    end; 
    For I:=1 To MA do
    begin
        For J:=1 To I do
        begin
            SUM1:=0;
            For K:=1 To MA do
                SUM1:=SUM1 + V[I, K] * V[J, K] * WTI[K];
            CVM[I, J]:=SUM1;
            CVM[J, I]:=SUM1;
        end; 
    end; 
end;

end.