Start with this:
p = input('Enter p '); % E.g. p=0.5
Y = 0; % Y <-- 0
C = 0; % C <-- 0
for t = 1 : 7
if rand < p
X(t) = 1;
else
X(t)=0;
if X(t) == 1
C = C + 1; %C <-- C+1
else
C = 0; % C <-- 0
if C >= 3
Y = 1; % Y <--1
end
end
end
end
% Output Y
Y
