function elem (n,k)
integer :: n,k
if (n == 0) then
elem = 1
else
elem = (-1.0)**n/((2*n+1)**(2*k+1))
endif
end function elem

program series
real :: sum, curelem, oldelem, eps, d1, d2, d3, d4
integer :: i

i = 1
oldelem = elem (0,1)
sum = oldelem
eps = 0.0001

curelem = 9990.0 !LOLWUT????
d1 = 100500.0
d2 = 100500.0
d3 = 100500.0
d1 = 100500.0

print*, "Step, summ, element, d1, d2, d3, d4"
print("(I6, F8.3, F8.3)"), 0, oldelem, sum

do while (d1>eps)

curelem = elem (i,1)
sum = sum + curelem
d1 = abs(curelem)
d2 = abs(curelem-oldelem)
d3 = abs(d1/sum)
d4 = abs(d2/sum)
print("(I6, F8.3, F8.3, F8.3, F8.3, F8.3, F8.3)"), i, sum, curelem, d1, d2, d3, d4
oldelem = curelem
i = i + 1
end do

print*, sum
end program