I use this to test: <div class="highlight highlight-source-python notranslate posi

Hello <a class="user-mention notranslate" data-hovercard-type="user" data-hovercard-ur

Show singularities at endpoints about differint HOT 1 CLOSED

differint commented on June 19, 2024

Show singularities at endpoints

from differint.

Comments (1)

cooperhatfield commented on June 19, 2024 1

Hello @CheapMeow, the asymptotic behavior at the left endpoint is a side effect of how the Grunwald-Letnikov differintegral is implemented. Comparing to the result gained from the Riemann-Liouville differintegral and the expected results, we see that it much more accurate.

These results were obtained using the following code:

import numpy as np
import differint.differint as df
import matplotlib.pyplot as plt

def f(t):
    return np.power(t, 3)
def neg_f(t):
    return f(-1 * t)


alpha = 0.6

t_2 = np.linspace(0, 1, 100)
GL_2 = df.GL(alpha, f, 0, 1, 100)

GL_2_rev = df.RL(alpha, f, 0, 1, 100)

expected = lambda t: 6/(df.Gamma(3-alpha+1)) * t ** (3 - alpha)


fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.plot(t_2, GL_2, label='GL')
ax.plot(t_2, GL_2_rev, marker='o', markevery=5, label='RL')
ax.plot(t_1, expected(t_1), label='expected')
ax.legend()

plt.show()

Another thing to point out, is that to conserve many of the properties one might expect from the fractional derivative of polynomials, the left endpoint should be kept to 0; I only included that case.

from differint.

Recommend Projects