n=0:16;                  
y=1./((n-2).^2+1)+1./((n-9^.2+4))+5;    #.运算      进行矩阵运算
plot(n,y,'mh','markersize',15)   #  m是品红色,h是 六角形,15是 六角形的大小
axis([0 17 5 6.2])     # 设置坐标轴的范围
grid on