Getting error in this code.

Question

AVINASH SAHU 2022-6-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1732525-getting-error-in-this-code

评论： Jan 2022-6-3

betasqr=0.2;
gamma=0.3;
psi = 0.001;
mu = 1.0;
sbar = 50;
A = @(x) (4+sbar)-(sbar*betasqr*gamma).*(x.*(1-x)).^0.5;
B = @(x) (1-(betasqr)*((x.*(1-x)).^(0.5)));
E = @(x) ((6*(2+sbar))-((2*betasqr*sbar*gamma)*((x.*(1-x)).^(0.5))));
G = @(x) (12*psi*(1+sbar)) + A/B;
C = @(x) E/G;
IntEbyG = integral(C,0,1)

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Jan 2022-6-2

编辑：Jan 2022-6-2

在 MATLAB Online 中打开

A simpler version - note that sqrt(x) is much cheaper than x^0.5:

betasqr = 0.2;
gamma   = 0.3;
psi     = 0.001;
mu      = 1.0;
sbar    = 50;
A = @(x) 4 + sbar - sbar*betasqr*gamma .* sqrt(x .* (1-x));
B = @(x) 1 - betasqr * sqrt(x .* (1-x));
E = @(x) 12 + 6 * sbar - 2*betasqr*sbar*gamma * sqrt(x .* (1-x));
G = @(x) 12 * psi * (1 + sbar) + A/B;  % ERROR
C = @(x) E/G;                          % ERROR
IntEbyG = integral(C,0,1)
Error using /
Arguments must be numeric, char, or logical.

Error in solution (line 9)
G = @(x) 12 * psi * (1 + sbar) + A/B;

Error in solution (line 10)
C = @(x) E(x)/G(x);

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 75)
        [q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 87)
Q = integralCalc(fun,a,b,opstruct);

I've displayed the error message. This is very useful, if you ask a question, because solving a problem is much easier than guessing, what the problem is. Whenever you mention an error in the forum, attach a copy of the complete message.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2022-6-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1732525-getting-error-in-this-code#answer_976875

在 MATLAB Online 中打开

A simpler version - note that sqrt(x) is much cheaper than x^0.5:

betasqr = 0.2;
gamma   = 0.3;
psi     = 0.001;
mu      = 1.0;
sbar    = 50;
A = @(x) 4 + sbar - sbar*betasqr*gamma .* sqrt(x .* (1-x));
B = @(x) 1 - betasqr * sqrt(x .* (1-x));
E = @(x) 12 + 6 * sbar - 2*betasqr*sbar*gamma * sqrt(x .* (1-x));
G = @(x) 12 * psi * (1 + sbar) + A(x) ./ B(x);  % <-- bug fix
C = @(x) E(x) ./ G(x); % <-- bug fix
IntEbyG = integral(C,0,1)
IntEbyG = 5.3442