Вот код MATLAB для выполнения Грамма Шмидта на странице 1http://web.mit.edu/18.06/www/Essays/gramschmidtmat.pdf

Я часами пытаюсь это сделать с помощью R, так как у меня нет MATLAB. Вот мой R

f=function(x){
m=nrow(x);
n=ncol(x);
Q=matrix(0,m,n);
R=matrix(0,n,n);

for(j in 1:n){
v=x[,j,drop=FALSE];

for(i in 1:j-1){
R[i,j]=t(Q[,i,drop=FALSE])%*%x[,j,drop=FALSE];
v=v-R[i,j]%*%Q[,i,drop=FALSE]
}

R[j,j]=max(svd(v)$d);
Q[,j,,drop=FALSE]=v/R[j,j]}

return(list(Q,R))}

Он продолжает говорить, что есть ошибки в любом:

v=v-R[i,j]%*%Q[,i,drop=FALSE] 

или же

R[j,j]=max(svd(v)$d);

Что я делаю неправильно, переводя код MATLAB в R ???