clc;
clear;
close all;
xaxis=[-50:1:50];yaxis=[-50:1:50];
[X,Y]=meshgrid(xaxis,yaxis);
A=X.^2+Y.^2;
A(find(A>2500))=0
%2500为半径的平方
运行结果:
0 0 0 0 0 25 0 0 0 0 0
0 0 25 20 17 16 17 20 25 0 0
0 25 18 13 10 9 10 13 18 25 0
0 20 13 8 5 4 5 8 13 20 0
0 17 10 5 2 1 2 5 10 17 0
25 16 9 4 1 0 1 4 9 16 25
0 17 10 5 2 1 2 5 10 17 0
0 20 13 8 5 4 5 8 13 20 0
0 25 18 13 10 9 10 13 18 25 0
0 0 25 20 17 16 17 20 25 0 0
0 0 0 0 0 25 0 0 0 0 0