抱歉，不晓得

A=[0.5,0.2,0.3;0.3,0.8,0.3;0.2,0,0.4];
x0=[p1,p2,p3];                %自己填数值
E=x0'*x0;
k=1;

while E>=0.00001
        temp=A*x0;
        printf('x%d =[%f %f %f]\n',temp(1),temp(2),temp(3));
        k=k+1;
        E=(temp-x0)'*(temp-x0);
        x0=temp;
end