лаба 1 по инф безопасности -движемся дальше

disp('=============');
 
U = [1,2,3,4,1];
V = reshape(U, 1, []);
p=seqperiod(V);
disp(p);
disp(U);

 
U = []; S = 0.5; n=1; 
for k = 1:1:10
    disp 'ля-ля-ля';
    disp (S); 
    S = (1/pi)*acos(cos(10^n)*S); U(k)=S;  
end; 
plot(U); 
 
V = reshape(U, 1, []);
p=seqperiod(V);
disp(p);
disp(U);
 
hold on ;
U = []; S = 0.00001; n=4; 
for k = 1:1:10   
    S = (1/pi)*acos(cos(10^n)*S); U(k)=S;   
end; plot(U);
 
V = reshape(U, 1, []);
p=seqperiod(V);
disp(p);
disp(U);

a = ppscq(0.1, 4);
disp('Предпериод 111 =');
disp(a);

a = ppscq(0.00001, 32);
disp('Предпериод 222 =');
disp(a);

a = ppscq(0.1, 2);
disp('Предпериод 333 =');
disp(a);

clear();
  figure
  [S, n] = meshgrid(1:5, 1:5);
  Z = ppscq(S,n);
  surf(S, n, Z);



%Z  = ppscq(S, n);
 % surf(S, n, f11(S,n));
   
   
  % fill3(S, n, ppscq(S, n),'b')
%Z = ppscq(S, n);
% disp ('Исследуем выход');
% R = ppscq([S, n]);
% disp (R);
%Z = R)./R;
%imagesc()
%mesh(Z);

  % fill3(S, n);
  

% for s = 0:0000.1:10   
%     for k = 1:1:10   
%         for k = 1:1:100   
%             S = (1/pi)*acos(cos(10^n)*S); U(k)=S;  
%         end;
%     end; 
% end; 
% x{end + 1} = a


%fill3(x, y, z, 'b')
 % scatter3(S(:), n(:), n+S+S ,'filled'), view(-60,100)
%view(40,35)
% ppscq(S,n)
% v = [1 2 3 4 5 12 2 1 6]
% preperiod(v);
% 
% clear();
%  figure
%  [S, n]=meshgrid(0:0.1:1, 1:1:3);
%  
%  surf(S, n, ppscq(S,n));


% clear();
% figure
% [x,y]=meshgrid(-3:1:4, -3:1:4);
% 
% surf(x,y,x^2+2*x+y^2-y);

обёртка для препериода:

function R = ppscq(S, n)
  
    U = []; 
    for k = 1:1:10
           disp 'ля-ля-ля';
    disp (S); 
        S = (1/pi)*acos(cos(10^n)*S); U(k)=S;  
    end; 

    R2 = S + n -S - n + preperiod(U);
    R = R2;
end