Info

此问题已关闭。 请重新打开它进行编辑或回答。

Help on double integral

1 次查看（过去 30 天）
Natnael Hamda
Natnael Hamda 2012-3-5
关闭： MATLAB Answer Bot 2021-8-20
I wrote a simple code for evaluating a double integral of a function of the following but it does n't work;
R = 0.0067; D = 1.1111;
f = inline('exp(-R*(y-x))*((((y+D)/(x+D)))^(R*(K+D)+1))*(1/((y+D)*y)))','x','y');
dblquad(f,0,100,100,inf)
Any with a solution? Thanks for your time
  1 个评论
Walter Roberson
Walter Roberson 2012-3-5
You have an extra ) at the end of the integral expression.

此问题已关闭。

回答（1 个）

Walter Roberson
Walter Roberson 2012-3-5
The above cannot be integrated, as it involves the undefined variable K .
  4 个评论
显示 2更早的评论
Natnael Hamda
Natnael Hamda 2012-3-6
Thanks a lot; The anonymous functions creates the function in better way. But still there is an error when i try to evaluate the double integral. I included the error message:
??? Error using ==> mldivide
Matrix dimensions must agree.
Error in ==> @(x,y)exp(-R*(y-x))*((((y+D)/(x+D)))^(R*(K+D)+1))*(1/((y+D)*y))
Error in ==> quad at 77
y = f(x, varargin{:});
Error in ==> dblquad>innerintegral at 78
Q(i) = quadf(intfcn, xmin, xmax, tol, trace, y(i), varargin{:});
Error in ==> quad at 77
y = f(x, varargin{:});
Error in ==> dblquad at 54
Q = quadf(@innerintegral, ymin, ymax, tol, trace, intfcn, ...
Mike Hosea
Mike Hosea 2012-3-8
Your integral function must be able to accept arrays and return arrays. Don't use *, as this is matrix multiply. Instead use .*, elementwise multiplication. Similarly, don't use ^, rather .^. Don't use /, rather ./ . Basically, your integrand isn't working because it is being passed arrays, and it's trying to do matrix multiplications, matrix powers, and linear system solves when what you really just want is simple arithmetic, elementy-by-element.

此问题已关闭。

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by