lambertw, all branches
|
01-23-2024, 03:54 PM
(This post was last modified: 01-23-2024 06:42 PM by Albert Chan.)
Post: #39
|
|||
|
|||
RE: lambertw, all branches
(01-23-2024 02:35 PM)Gil Wrote: For W (k=-1, a=+3*PI/2), p2> a = 3*pi/2 p2> a mpf('4.7123889803846897') p2> W(a, -1) mpc(real='-1.4351953215276557e-16', imag='-4.7123889803846897') Note that a is not *exactly* \( 3\pi/2 \), but only an approximation. Even with more precisions, W(a, -1) will have tiny real part. We accept this numerical limits, and confirm result instead. x = -a*I e^x = cis(-a) = cis(-a+2*pi) = cis(pi/2) = I x * e^x = -a*I * I = a |
|||
« Next Oldest | Next Newest »
|
User(s) browsing this thread: