pwitrtexample.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 6 行

clear;clc;
A=[10,2,3;2,20,4;3,4,30];
N=10;
lmd=PwItrt(A,N);
lmd0=max(eig(A));
plot(1:N,lmd,1:N,ones(1,N)*lmd0);