matlab小程序

if u<=1是对u这个矩阵在进行判断，所以其判定结果是0，直接进行下一步else运算，及对u中的元素开方。

对一个矩阵中的具体元素做判断，不能用if u<=1来进行

这是因为用ezplot画图，默认的范围是[-2*pi,2*pi],当然画不出大于2*pi以外的图形，我改一下：

function r=yuan(a,b,c)

%yuan.m,给定参数绘制圆

%输入（a,b,c）分别为圆心横纵坐标，半径

%输出（a,b）为圆心，c为半径的圆

syms x y eq

eq=(x-a)^2+(y-b)^2==c^2

%r=ezplot(eq)

r=ezplot(eq,[(a-2*c) (a+2*c) (b-2*c) (b+2*c)])

set(r,'linestyle','--','color','r','linewidth',2,'fill','on')

axis auto

axis equal

axis fill

%axis([a-c a+c+1 b-c b+c+1])

end

x2,y2长度不同

y1，y2长度相同

建议改下 不知道你要求x2与y2关系如何

x1=-17:1:3

y1=1./((x1+3).^2+1)+1./((x1+9).^2+4)+5

x2=-17:0.02:3

y2=1./((x2+3).^2+1)+1./((x2+9).^2+4)+5

或者再加上y2=1./((y2+3).^2+1)+1./((y2+9).^2+4)+5

---