function g=fct(x)
g=1;
if x>0
   for i=1:x
     g=g*i;
   end
end
