[laosam280]

x=[-2 -1 0 1 2 3];
y=x.^2;
for i=1:1:6
a=x(i)+cos([-pi:0.01i]);
b=y(i)+sin([-pi:0.01i]);
plot(a,b);
hold on
end
plot(x,y+1,'.');
hold on
plot(x,y-1,'-');
基本的运算就是这样的，当然我这个程序没有写出包络。但是过程如下：
首先找到圆心连线的方程，然后包络是所有圆心的连线的方程进行平移得到的。