#놈
norm <- function(x){
     return(max(abs(x)));
}

#오차허용한계에 맞게 자리수를 맞춰주는 함수
figure <- function(x,TOL){
    round(x,abs(log10(TOL)));
}

# 행렬 A; 벡터 x; 오차허용한계 TOL; 최대반복수 M
powerMethod <- function(A,x,TOL,M){
    n <- length(A[,1]); #행렬의 차원
    y <- c();
    mu <- 0; #지배고유치를 저장할 변수
    ERR <- c();     

    #step1
    k=1;

    #step2
    p <- 0;
    for(i in 1:n){
        if(abs(x[i])==norm(x)){
            p <- i;
            break;
        }
    }

    #step3
    x <- x/x[p];

    #step4
    while(k<=M){

        #step5
        y <- A%*%x;    

        #step6
        mu <- y[p];

        #step7    
        for(i in 1:n){
            if(abs(y[i])==norm(y)){
                p <- i;
                break;
            }
        }        

        #step8
        if(y[p]==0){
            cat("Eigenvector",x,"\\n");
            cat("A has the eigenvalue 0, select a new vector x and restart");
            break;
        }

        #step9
        ERR <- norm(x-y/y[p]);
        x <- y/y[p];
    
        #step10    
        if(ERR<TOL){    
            cat("지배고유치: ",figure(mu,TOL),"\\n");
            cat("고유벡터: ",figure(x,TOL),"\\n");
            break;
        }
    
        #step11    
        k <- k+1;
    }

    #step12
    if(k>M){
        cat("Maximum number of iterations exceed\\n");
    }
    
}


A <- matrix(c(-4,-5,-1,14,13,0,0,0,2), ncol=3); #행렬 A
x <- seq(1,1,length=ncol(A));
TOL <- 1.0e-4;
M <- 30;
powerMethod(A,x,TOL,M);



자바 해본 적이 없어서 어케해야할지 모르겠다;;;;
살려주세요