clear all
%y(1)=6999995; %large positive number leading to positive infinitive
%y(1)=-59999994; %large negative number leading to negative infinitive
y(1)=65;

%When using while loop, comment out for loop statements
%k=1;
%while (y(k)~=1);
%    y(k+1)=((y(k)*3+1)*(sin((y(k).*pi)/2).^2))+((cos((y(k).*pi)/2).^2)*(y(k)/2));
%    k=k+1;
%end

%When using for loop, comment out while loop statements
for k=1:40
    y(k+1)=((y(k)*3+1)*(sin((y(k).*pi)/2).^2))+((cos((y(k).*pi)/2).^2)*(y(k)/2));
end

x=1:length(y);

%for large scales use loglog
loglog(x,y,'r*')

%for smales scales use plot
%plot(x,y,'r*')

axis square
grid
clc