Inconsistent results from "int" function

2022 11 27
1 个回答
回答已采纳
更新时间：2022 11 27
4 次查看（30 天）

0 个投票

Hi, everyone. For the same question, the "int" function provided me with different results.
syms x y
int(x*y/24+y^2/48,x)
int((y*(2*x + y))/48,x)
The second line gave me "(x*y*(x + y))/48" which was what I originally wanted,
but the third line gave me "'(y*(2*x + y)^2)/192" which is wrong.
Can anyone explain the reason?
Thanks for your time and help in advance!

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

 采纳的回答

Alan Stevens
Alan Stevens 2022-11-27

2 个投票

You can express '(y*(2*x + y)^2)/192 as (x*y*(x + y))/48+y^3/198
As far as the integral is concerned y is a constant, so the two results are the same to within a constant, which is all you can guarantee from an integral with no specific limits.

5 个评论
显示 3更早的评论 隐藏 3更早的评论

YoungJoo Lee
YoungJoo Lee 2022-11-27
Dear Alan,
Thanks for your kind answer. I am also aware that the problem is about an indefinite integral, and the result may get involved with a constant. However, I think there is no mathematically-valid reason that the constant can be specified as y^3/192. If this is how the "int" function works, it seems the function cannot be used for indefinite integral problems...
Torsten
Torsten 2022-11-27
Ran in:
Ran in:
However, I think there is no mathematically-valid reason that the constant can be specified as y^3/192.
There is a mathematically-valid reason that each function f that only depends on y can be specified as the constant of integration. The reason is d/dx(f(y)) = 0.
MATLAB does not perform an "evaluation of the most beautiful antiderivative".
In this case, you can use
syms x y
int(x*y/24+y^2/48,x,0,x)
ans = 
int((y*(2*x + y))/48,x,0,x)
ans = 
But in other examples, the lower limit might be different from 0 to get a "nice" antiderivative.
YoungJoo Lee
YoungJoo Lee 2022-11-27
Dear Torsten,
I appreciate you on your detailed explanation. I agree with you that "any" function f that only depends on y can be the constant of the integration in my problem. However, I am still doubting whether it is reasonable to provide only one specific function for y (e.g., y^3/192 in my case) as output, as if it were the only and unique answer...
I believe, for my problem, providing something like "(x*y*(x + y))/48+C (where C is a constant)" as output would be more appropriate.
Anyway, I got your point and also thank you for sharing your tip.
Torsten
Torsten 2022-11-27
I believe, for my problem, providing something like "(x*y*(x + y))/48+C (where C is a constant)" as output would be more appropriate.
"(x*y*(x + y))/48+C(y) (where C is a function solely depending on y)"
YoungJoo Lee
YoungJoo Lee 2022-11-27
Oh, sure. Thanks for the correction!

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Number Theory 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by