r=-30:30;
z=(25/118)*exp(-59/10000*r.^2)/pi^(1/2)*59^(1/2);
t=0:1:32;
w=pi;

for i=1:33
V=@(t)(5+0.2*t);
if cos(w*t(i))==1
   R=r+quad(V,0,t(i));
else 
   R=r;
end
Z=z;
plot(R,Z);
hold on
end
hold off
可以运行,但总是
提示:Warning: Minimum step size reached; singularity possible.
还有就是看看怎么能把图中的各个图形叠加起来,为的是看看最终的轮廓是什么样子?。。请路过的大侠帮帮忙。。。在这里先谢过了。。。