Wrong solution of differential equation using symbolic lambda

2 ビュー (過去 30 日間)
Kevin Oyen
Kevin Oyen 2021 年 11 月 23 日
編集済み: Srijith Kasaragod 2021 年 12 月 2 日
As can be seen in the screenshot I have a problem with the symbol lambda in matlab R2021B. If lambda is used instead of the variable L, a wrong solution for the differential equation is obtained. What can be the cause for this problem? Many thanks in advance!
clear all
syms R(r) L
assume(L, "real")
assume(L > 0)
vgl = r*diff(R(r),2,r) + diff(R(r),1,r) == L*r*R(r)
dsolve(vgl)
clear all
syms R(r) lambda
assume(lambda, "real")
assume(lambda > 0)
vgl = r*diff(R(r),2,r) + diff(R(r),1,r) == lambda*r*R(r)
dsolve(vgl)
  3 件のコメント
Paul
Paul 2021 年 11 月 23 日
I don't have an answer, but the result seems (strangely) to depend on the case (upper or lower) of the first character of the variable. Is there way you can check if the solutions are equivalent? I wasn't sure how to choose the Ci.
syms R(r) L
assume(L, "real")
assume(L > 0)
vgl = r*diff(R(r),2,r) + diff(R(r),1,r) == L*r*R(r);
dsolve(vgl)
ans = 
syms R(r) lambda
assume(lambda, "real")
assume(lambda > 0)
vgl = r*diff(R(r),2,r) + diff(R(r),1,r) == lambda*r*R(r);
dsolve(vgl)
ans = 
syms R(r) Lambda
assume(Lambda, "real")
assume(Lambda > 0)
vgl = r*diff(R(r),2,r) + diff(R(r),1,r) == Lambda*r*R(r);
dsolve(vgl)
ans = 
syms R(r) AAA
assume(AAA, "real")
assume(AAA > 0)
vgl = r*diff(R(r),2,r) + diff(R(r),1,r) == AAA*r*R(r);
dsolve(vgl)
ans = 
syms R(r) aAA
assume(aAA, "real")
assume(aAA > 0)
vgl = r*diff(R(r),2,r) + diff(R(r),1,r) == aAA*r*R(r);
dsolve(vgl)
ans = 

サインインしてコメントする。

採用された回答

Srijith Kasaragod
Srijith Kasaragod 2021 年 11 月 30 日
編集済み: Srijith Kasaragod 2021 年 12 月 2 日
Hi Kevin,
This is a bug and has been brought to the notice of our developers. It may be fixed in future releases. One possibility to resolve this problem would be to avoid solutions using the imaginary unit in representation, which in this case would prefer the representation using besseli and besselk functions.
Regards,
Srijith.
  1 件のコメント
Paul
Paul 2021 年 11 月 30 日
Can you better describe what the bug actually is? What is the title and description of the bug? Is there a link to a bug report?

サインインしてコメントする。

その他の回答 (0 件)

製品


リリース

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by