Can you be more specific what you mean by "the current gaussian curve is very wrong"? Do you mean with the lines connecting all over the place? If you just plot points rather than lines, then you see the nice smooth curve:
options = ['A', 'B'];
%Empty vectors for value storage
RunSumA1=[];
RunSumB1=[];
ListOneProbability=0.013;
ListTwoProbability=0.987;
totalRuns = 100;
%Storage loop from function
for k=1:totalRuns
[RunSumA1,RunSumB1]= my_sorter(RunSumA1,RunSumB1,options,ListOneProbability,ListTwoProbability);
end
%Converting to Column Vector
RunSumA1 = RunSumA1';
%Average, STD, Normal PDF
mu = mean(RunSumA1);
sigma = std(RunSumA1);
y = normpdf(RunSumA1, mu, sigma);
%Plot
plot(RunSumA1, y, '.') % Changed this line to plot only the points
%Random choice function
function [RunSumA1,RunSumB1]= my_sorter(RunSumA1,RunSumB1,options,ListOneProbability,ListTwoProbability)
tempA1 = 0;
tempB1 = 0;
for j=1:5000 % I made the smaller, since it would not time out on the Answers forum
newChoice = randsample(options, 1,true, [ListOneProbability,ListTwoProbability]);
if newChoice == 'A'
tempA1=tempA1+1;
else
tempB1=tempB1+1;
end
end
RunSumA1=[RunSumA1,tempA1];
RunSumB1=[RunSumB1,tempB1];
end
Is that it? Or is there some other problem?
