I'm trying to integrate the following function from -inf to +inf with repsect to x.
fun = @(x) -3.*v.^2 -4.*v.*a./((x.^2+a.^2)*pi) - a.^2./((x.^2+a.^2)^2*pi^2)fun = function_handle with value: @(x)-3.*v.^2-4.*v.*a./((x.^2+a.^2)*pi)-a.^2./((x.^2+a.^2)^2*pi^2)>> integral(fun,-inf,inf)Error using symengineNot a square matrix.Error in sym/privBinaryOp (line 1030) Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:});Error in ^ (line 330) B = privBinaryOp(A, p, 'symobj::mpower');Error in @(x)-3.*v.^2-4.*v.*a./((x.^2+a.^2)*pi)-a.^2./((x.^2+a.^2)^2*pi^2)Error in integralCalc/iterateScalarValued (line 314) fx = FUN(t);Error in integralCalc/vadapt (line 132) [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);Error in integralCalc (line 103) [q,errbnd] = vadapt(@minusInfToInfInvTransform,interval);Error in integral (line 88)Q = integralCalc(fun,a,b,opstruct);
I don't know why i'm getting this error?
by the way v and a are both constants
Best Answer