#!/usr/bin/env julia

# Jacek Poziemski 272389

"""
    kahanmacheps(T::Type{<: AbstractFloat})::T

Calculate the Kahan machine epsilon
using the following expression: `3(4/3 - 1) - 1`.

# Arguments
- `T`: the floating point type
"""
kahanmacheps(T::Type{<: AbstractFloat})::T = T(3) * ((T(4) / T(3)) - one(T)) - one(T)

for T::Type{<: AbstractFloat} in [Float16, Float32, Float64]
    println("kahanmacheps($T): $(kahanmacheps(T))")
    println("eps($T): $(eps(T))")

    println()
end