幂法求矩阵的主特征值和特征向量

#include<iostream>
using namespace std;
int main()
{
    double A[10]={1.0,1.0,0.5,1.0,1.0,0.25,0.5,0.25,2.0};//矩阵
    double v[3];
    double vk=-1;
    double uk[3];
    v[0]=1;
    v[1]=1;
    v[2]=1;
    uk[0]=uk[1]=uk[2]=1;
    int NUM=0;  //迭代次数
    while(NUM<10)
    {

    for(int i=0;i<3;i++)
    {
        double r=0;
        for(int j=i*3;j<i*3+3;j++)
        {
            r+=A[j]*uk[j%3];
        }
        v[i]=r;
    }
    vk=0;
    for(int i=0;i<3;i++)
    {
        if(vk<v[i])
        {
            vk=v[i];
        }
    }
    for(int i=0;i<3;i++)
    {
        uk[i]=v[i]/vk;
        cout<<uk[i]<<" ";
    }
    cout<<endl;
    NUM++;

        }

    return 0;


 }