if (verze == 3) then
    obsah = 0.0
    do i=0,10
        sp = 1
        DO o = 1, i
            sp = sp *o
        END DO

        delta = (2.0**(i+0.5))/(4.0*sp)
        obsah = obsah +  delta * gamma((2.0*i-1.0)/4.0) * (eps/(kb*T))**((2.0*i+1.0)/4.0)

    end do
    S = - (2.0 * pi)/3.0 *Na* sigma**3 * obsah * 1e4 * 100
end if