想要画一个图，y=e.^2+sqrt(r)从（0.1，6）求积分，想要看积分值同e和r的关系
function f=fun(x)
global e;
global r;
f=e.^2.*x+r.*sqrt(x);

主程序：
global e;               
global r;
y=ones(5,10);                先赋予y一定的空间
i=1;
for e=1:10
   for r=1:5
     y(i)=quad('fun',1,6);
     i=i+1;
   end
end
e=1:10;
r=1:5;
[ee,rr]=meshgrid(e,r);
mesh(ee,rr,y);


因为quad不能求向量的积分，所以只好用了两个循环，请高手指导一下这样的程序是否有错，因为得出的结果不是想要的图形。