编写迭代函数factorial，输入一个数，返回该数的阶乘。
function k=factorial(n)
if nargin~=1,error;end
if nargout>1,error;end
if n error;
end
if n>1
k=n*factorial(n-1);
elseif n==0|n==1
k=1;
end
运行：
>> factorial(5)
ans =
120
>> factorial(10)
ans =
3628800
法二：利用内置函数prod( )
function k=my_fact(n)
prod(1:n)
end
运行：
my_fact(11)
ans =
39916800
实例5：单个轴窗口显示多个图形
function shili05
h0=figure('toolbar','none',...
'position',[200 150 450 250],...
'name','实例05');
t=0:pi/10:2*pi;
[x,y]=meshgrid(t);
subplot(2,2,1)
plot(sin(t),cos(t))
axis equal
subplot(2,2,2)
z=sin(x)-cos(y);
plot(t,z)
axis([0 2*pi -2 2])
subplot(2,2,3)
h=sin(x)+cos(y);
plot(t,h)
axis([0 2*pi -2 2])
subplot(2,2,4)
g=(sin(x).^2)-(cos(y).^2);
plot(t,g)
