Commit d2070230 authored by ulrich_y's avatar ulrich_y

Follow-up to a7004f12

parent a7004f12
......@@ -255,8 +255,8 @@ CONTAINS
FUNCTION BERNOULLI_POLYNOMIAL(n, x) result(res)
integer, parameter :: maxn = 15
integer n
restype :: x, res
restype :: xpow(maxn+1)
complex(kind=prec) :: x, res
complex(kind=prec) :: xpow(maxn+1)
integer, parameter :: coeffN(maxn+1, maxn) = reshape((/ &
- 1, + 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, &
+ 1, - 1, + 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, &
......@@ -292,7 +292,7 @@ CONTAINS
+ 1, + 2, + 1, + 6, + 1, + 2, + 1, + 2, + 1, + 6, + 1, + 2, + 1, + 2, + 2, + 1 /), &
(/maxn+1, maxn/))
real(kind=prec), parameter :: coeff(maxn+1,maxn) = coeffN/real(coeffD)
real(kind=prec), parameter :: coeff(maxn+1,maxn) = coeffN/real(coeffD,kind=prec)
integer i
if (n>maxn) then
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment