function y=factsec4(x)
  y=1;
  for i=1:x
    y=i*y;
  end
  
end
