pfqz.m

来自「Matlab语言常用算法程序集.rar 里面有很多有用的Matlab数值运」· M 代码 · 共 10 行

function r = PFQZ(k,x0,n)
r = zeros(n,1);
x = zeros(n,1);
x(1) = x0;
r(1) = x(1)/(100^k);

for i=2:n
    x(i) = mod(x(i-1)^2 /(10^k),100^k);
    r(i) = x(i)/(100^k);
end