Complex Trilogarithm

版本 1.0.1 (2.3 KB) 作者: Didier Clamond
= Spence's tri-Logarithm.
0.0
(0)
1.0 次下载
更新时间 2023/9/7

查看许可证

function t = trilog(z)
trilog = Spence's tri-Logarithm.
t = trilog(z) = Li_3(z)
= Sum From n=1 To n=Infinity Of z^n/n^3 if |z|<=1.
INPUT z: real or complex, scalar, vector or matrix.
OUTPUT t: component-wise trilogarithm of z (accuracy around machine precision).
Depending on the size of the matrix z and of the magnitude of its
components, trilog(z) is 2 to 5 (and more) orders of magnitude
faster than the polylog(3,z) from the Matlab symbolic toolbox.
EXAMPLE: trilog(z) vs polylog(3,z).
N = 10000;
z = 200*( rand(N) + i*rand(N) ) - 100 - i*100;
tic; t=trilog(z); toc;
Elapsed time is 12.081650 seconds.
N = 160;
z = 200*( rand(N) + i*rand(N) ) - 100 - i*100;
tic; t=polylog(3,z); toc;
Elapsed time is 12.289606 seconds.

引用格式

Didier Clamond (2024). Complex Trilogarithm (https://www.mathworks.com/matlabcentral/fileexchange/134986-complex-trilogarithm), MATLAB Central File Exchange. 检索时间: .

MATLAB 版本兼容性
创建方式 R2023a
兼容任何版本
平台兼容性
Windows macOS Linux
标签 添加标签

Community Treasure Hunt

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

Start Hunting!

版本 已发布 发行说明
1.0.1

Update title
1.0.0