ＭＡＴＬＡＢで下記のプログラムです。
これを、ＶＣ＋＋６に書き直すとどうなるのでしょうか？

Normは行列のノルム、rankはやはり行列の階数だと思うのですが、
このへんの演算も、ＶＣ＋＋６では
関数で書かないといけないのでしょうか？
nergin < とかisempty()とか
svd()とか、sDiag()とか（対角行列か？）
いきなり括弧でくくる[U,S,V] = とかB = []とか
このへんも、何をやっているのかわかりません。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function[Y,W,b,M] = abc(X, numIC,econv,max_itr)
% Initialization

if (nergin < 5), display = 0; end
if isempty(display), display = 0; end

if nargin < 4, max_itr = 1000, end
if isempty(max_itr), display = 1000; end

if nargin < 3, econv = 0.00001; end
if isempty(numIC), numIC = 0; end

if nargin < 2, numIC = 0; end
if isempty(numIC), numIC = 0; end
[rowX,colX] = size(X);
if unmIC == 0 | numIC>rank(X) numIC=rank(X); end

% Globe Transform
[U,S,V] = svd(X,0)
Sdiag = diag(S(1:numIC,1:numIC));
Sdiag0 = find(Sdiag == 0)
for i=;length(Sdiag0)
Sdiag(1,Sdiag0(i) = inf;
end

M = V(:,1:numIC)*diag(1./Sdiag)*sqrt(rowX); % shering
U = U(:,1:numIC)*sqrt(rowX); % normizing PCs (scores)
clear S V Sdiag Sdiag0

%------ ABC ---------
B = [];
irr_fail = 0;
i = 1;
while i <=numIC
w_old = rand(numIC, 1) - 0.5;
if i > 1, w_old = w_old-B*B'*w_old; end
w_old = w_old/norm(w_old);
for itr=1:max_itr
w = (U'*((U*w_old).^3))/roxX-3*w_old;
if i>1
w = w -B*B'*w;
end
w = w/norm(w);
if min(norm(w-w_old),norm(w+w_old))<econv
irr_fail = 0;
i = i + 1;
B = [B w];
break
end
w_old = w;
if isr == 1000
itr_fail = itr_fail + 1;
irr_fail = 5,error('Error !'), end
break
end
end
end

% Result
Y = U*B;
W = M*B;