求助，总出错

MC夜雨冷风

clear;
%close all;
clc;

for T=-10:1:10
T=10.^(T./10);
p1=[];
syms r;
f4=2.*pi.*0.25.*r.*exp(-pi.*0.25.*r.^2).*E_S_mw(T,r);
f5=matlabFunction(f4);
p=quadgk(f5,0,+Inf);
end

function [output] = E_S_mw(T,r)
        T=T;
        syms u x
        f= 1./(1+T.*u.^2.*x.^(-4)).*pi.*0.25.*exp(-0.25.*pi.*u);
        f1=matlabFunction(f);
        y1=quadgk(f1,0,+Inf);

        syms x
        r=r;
        f2=exp(-2.*pi.*0.25.*(1-y1).*x);
        f3=matlabFunction(f2);
        y2=quadgk(f3,r,+Inf);
        output=y2;
end

深度大脑