function htsim(n, p)

% Simulate the game of heads or tails, tossing n times
% with probability p that heads comes up and plot Peter's cumulative
% wins/losses


flip=rand(1,n)<p; % create n-component vector where each component is 0 or 1,
% 0=tails; 1=heads
x=2*flip-1; % create a new vector for each 1 (head) in flip, 1, for each 0 (tail), -1 
fortune=cumsum(x); % the jth component is the sum of the first j components in x

plot(0:n,[0, cumsum(x)])