diff --git a/test/TestIEEE754.py b/test/TestIEEE754.py
new file mode 100755
index 0000000..44c5236
--- /dev/null
+++ b/test/TestIEEE754.py
@@ -0,0 +1,28 @@
+#!/usr/bin/python3
+
+import struct
+from fractions import Fraction
+
+def to_float_binary(value):
+    # Use 'struct.pack' to pack the value into bytes using IEEE 754 floating-point format
+    # '>f' specifies big-endian single-precision float. Change to '>d' for double precision.
+    packed = struct.pack('>f', value)
+    # Convert the bytes to an integer and then to a binary string
+    return ''.join(f'{byte:08b}' for byte in packed)
+
+def to_float_bits(number) :
+    float_bytes = struct.pack('>f', number)
+    return struct.unpack('>I', float_bytes)[0], float_bytes
+
+def fixedPoint(n, shift) :
+    if shift >= 0 :
+        return Fraction(n) * (1<<shift)
+    else :
+        return Fraction(n) / (1<<-shift)
+
+number = fixedPoint(0xffffff, 127-23) # Largest number
+number = fixedPoint(0x1, -128) # Smallest number
+
+float_int32, float_bytes = to_float_bits(number)
+print ("float hex:",  float_bytes.hex(), hex(int(number*(1<<128))))
+print ("float hex:",  hex(float_int32), hex(int(number*(1<<128))))