simplerandomwalk

n=50
p=@(m) factorial(n)/((2^n)*factorial((n+m)/2)*factorial((n-m)/2))
b=@(m) sqrt(2/(pi*n))*exp(-(m^2)/(2*n))
x=[10:2:40];
y=[];
z=[]
for i=10:2:40;
    a=p(i);
    y=[y a]
    z=[z b(i)]
end

plot(x,y)
hold on

plot(x,z)
hold off