doublement.m

来自「3种伪随机数列发生器的matlab编码以及其各种特性的分析」· M 代码 · 共 10 行

function res = doublement(u0,n)
u = zeros(n,1); %initialisation du tableau avec des zeros
u(1)=u0; %initialisation de la recurrence
 
for k=1:n-1; %recurrence
    u(k+1)=2*u(k)^2-1;
end

res=u;