Matlab- Factorial

ch='y';
while ch~='n'
n=input('Pick a number:\n');
fact=zeros(1,n+1);
fact(1)=1;
fact(2)=1;
k=2;
while k<=n
    fact(k+1)=k*fact(k);
    k=k+1;
end
fprintf('factorial of %d is\n',n);
fprintf('%g\n',fact(n+1));
ch=input('Do you want fo continue: ');
end