diff --git a/esbo_etc/__init__.py b/esbo_etc/__init__.py
index 6aaf515..e3624c0 100644
--- a/esbo_etc/__init__.py
+++ b/esbo_etc/__init__.py
@@ -1,4 +1,4 @@
 from esbo_etc.classes import *
-from esbo_etc.modules import *
+from esbo_etc.lib import *
 
 # __root__ = "."
diff --git a/esbo_etc/lib/__init__.py b/esbo_etc/lib/__init__.py
index 42fc95b..7550af5 100644
--- a/esbo_etc/lib/__init__.py
+++ b/esbo_etc/lib/__init__.py
@@ -1 +1,2 @@
-from esbo_etc.lib.helpers import *
+from .helpers import *
+from .output import *
diff --git a/esbo_etc/lib/output.py b/esbo_etc/lib/output.py
new file mode 100644
index 0000000..3b1dbf7
--- /dev/null
+++ b/esbo_etc/lib/output.py
@@ -0,0 +1,38 @@
+from rich.console import Console
+from rich.table import Table
+import astropy.units as u
+
+
+def printSNR(exp_time: u.Quantity, snr: u.Quantity):
+    table = Table(show_header=True, header_style="bold magenta")
+    table.add_column("#", style="dim", width=4, justify="center")
+    table.add_column("Exposure Time", justify="right")
+    table.add_column("SNR", justify="right")
+    for i, exp_time_, snr_ in zip(range(len(exp_time)), exp_time.value, snr.value):
+        table.add_row(str(i), ("%1.4e " + exp_time.unit.to_string()) % exp_time_, "%1.4e" % snr_)
+    console = Console()
+    console.print(table)
+
+
+def printExposureTime(exp_time: u.Quantity, snr: u.Quantity):
+    table = Table(show_header=True, header_style="bold magenta")
+    table.add_column("#", style="dim", width=4, justify="center")
+    table.add_column("SNR", justify="right")
+    table.add_column("Exposure Time", justify="right")
+    for i, exp_time_, snr_ in zip(range(len(exp_time)), exp_time.value, snr.value):
+        table.add_row(str(i), "%1.4e" % snr_, ("%1.4e " + exp_time.unit.to_string()) % exp_time_)
+    console = Console()
+    console.print(table)
+
+
+def printSensitivity(exp_time: u.Quantity, snr: u.Quantity, sensitivity: u.Quantity):
+    table = Table(show_header=True, header_style="bold magenta")
+    table.add_column("#", style="dim", width=4, justify="center")
+    table.add_column("Exposure Time", justify="right")
+    table.add_column("SNR", justify="right")
+    table.add_column("Sensitivity", justify="right")
+    for i, exp_time_, snr_, sensitivity_ in zip(range(len(exp_time)), exp_time.value, snr.value, sensitivity.value):
+        table.add_row(str(i), ("%1.4e " + exp_time.unit.to_string()) % exp_time_, "%1.4e" % snr_,
+                      ("%1.4e " + sensitivity.unit.to_string()) % sensitivity_)
+    console = Console()
+    console.print(table)