matlab拟合曲线积分,求助大神，关于MATLAB带积分的拟合

Vc=3.33e-5;%clearance volume m3

B=86e-3;

l=175e-3;

a=43.035e-3;

nc=1.30;

x=[0.942187500000000        101.325000000000        389.827000000000        0.000533263830555081        3000        9.42477796076938        -0.0709441434350655;

0.942187500000000        101.325000000000        402.817000000000        0.000487559323265132        3000        8.72664625997165        -0.0700784601260763;

0.753750000000000        101.325000000000        377.300000000000        0.000533263830555081        3000        9.42477796076938        -0.0715951512460594;

0.753750000000000        101.325000000000        394.932000000000        0.000487559323265132        3000        8.72664625997165        -0.0702303037709998;

0.753750000000000        111.457500000000        388.831000000000        0.000533263830555081        3000        9.42477796076938        -0.0707556478758501;

0.753750000000000        111.457500000000        401.750000000000        0.000487559323265132        3000        8.72664625997165        -0.0697957167872532;

0.753750000000000        101.325000000000        389.970000000000        0.000533263830555081        3000        9.42477796076938        -0.0699301071396568;

0.753750000000000        101.325000000000        403.150000000000        0.000487559323265132        3000        8.72664625997165        -0.0685914396033772;

0.753750000000000        101.325000000000        381.380000000000        0.000487559323265132        2000        8.72664625997165        -0.0800233462039400];

b=[-0.0073 1.4829 110.3672 -0.1488];

y=[1 1 1 1 1 1 1];

f = @(b,x) quadl(x(:,1).^b(2)./((b(1).*exp((b(3).*x(:,2).*(x(:,4)./(Vc+((pi.*B.^2)./4).*(l+a-a.*cos(theta)-sqrt(l.^2-(a.^2).*sin(theta).^2)))).^nc).^b(4))./ ...

(x(:,3).*(x(:,4)./(Vc+((pi.*B.^2)./4).*(l+a-a.*cos(theta)-sqrt(l.^2-(a.^2).*sin(theta).^2)))).^(nc-1))).*x(:,5)),x(:,6),x(:,7));

b=lsqcurvefit(f,b,x,y);