simplerandomwalk2

n=input('n: ');
b=@(a,c,m) 2*n+m-4*c-3*a;
d=@(a,c,m) 2*a+3*c-m-n;
p=@(a,c,m) factorial(n)/((4^n)*factorial(a)*factorial(c)*factorial(b(a,c,m))*factorial(d(a,c,m)));
x=[0:2*n];
y=[];
for m=0:2*n
    matriz=posibles(n,m);
    ptotal=0;
    for i=1:size(matriz,1)
        a=matriz(i,1);
        c=matriz(i,2);
        ptotal=ptotal+p(a,c,m);
    end
    y=[y ptotal];
end

plot(x,y);
hold on
syms f(x);
f(x)=sqrt(2/(10*pi*n))*exp(-(x^2)/(5*n));

fplot(f,[0,100]);
hold off