求助matlab积分，其中积分上限是变化的

南北左右

我尝试编了程序，但是一直运行都没有结果，不知道是不是由于for循环太多了 还是程序有错。

最后希望得到的就是面积M,N关于x的一个矩阵



function f=y1(x)
A1=3.84333;
A2=84.19866 ;
LOGx0= -0.04921;
p= -0.05426;
f=A1+(A2-A1)./(1+10.^((LOGx0-x).*p));

%%%%%%%%%%%%%%%%%%%%
A1=3.84333;
A2=84.19866 ;
LOGx0= -0.04921;
p= -0.05426;

m=-45;%积分下限
n=45;%
x=m:0.625:n;
y1=A1+(A2-A1)./(1 + 10.^((LOGx0-x).*p));

%%%%%%%x<0
intg1=zeros((n-m)/0.625+1,1);
i=1;
for x=m:0.625:0;
    intg1(i,1)=quad('y1',m,x);
    i=i+1;
M=A2*(-m)-(x-y1).*y1-intg1;
end

%%%%%%%x>0
intg2=quad(‘y1’,m,0);    %求出可得3.343490936379536e+03
for x=0:0.625:n;
    intg3(i,1)=quad('y1',0,x);
    s=x*y1；
   N=A2*(-m)-intg2-(intg3-s);
    i=i+1;
end