From 1c71badbefb55bf08bd31c4462b6a19480115fe2 Mon Sep 17 00:00:00 2001 From: Joakim Skogholt Date: Sat, 20 May 2023 08:36:15 +0200 Subject: [PATCH] Fixed another error in GL derivative --- src/TR.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/TR.jl b/src/TR.jl index f005cba..0abfa2c 100644 --- a/src/TR.jl +++ b/src/TR.jl @@ -76,7 +76,7 @@ elseif regType == "GL" # Grünwald-Letnikov fractional derivative regulariztion # regParam1 is alpha (order of fractional derivative) C = ones(p)*1.0; for k in 2:p - C[k] = (-1)^(k-1) * (1-(regParam1+1)/k) * C[k-1]; + C[k] = (1-(regParam1+1)/(k-1)) * C[k-1]; end regMat = zeros(p,p);