x取范圍1到10000，y取范圍10001到20000，計算z=x^2+y2的值。分別使用for循環(huán)，meshgrid后的矩陣進行運算，利用tic，toc記錄每個方法所需要的時間。

x = 1:5000;

y = 5001:10000;

z = zeros(5000);

ans0 = zeros(5000,5000);

tic

for i = 1:5000

    for j = 1:5000

    ans0(i,j) = x(i)^2 + y(j)^2;

    end

end

toc                                            %0.491108s

tic

[X,Y] = meshgrid(x,y);

ans1 = X.^2 + Y.^2;

toc                                            %0.187623s

運算時間如標(biāo)注所示，利用矩陣運算能節(jié)省一半以上的時間。