diff --git a/tests/.gitignore b/tests/.gitignore
new file mode 100644
index 0000000..1d97a8a
--- /dev/null
+++ b/tests/.gitignore
@@ -0,0 +1,2 @@
+/output.actual
+/test
diff --git a/tests/Makefile b/tests/Makefile
new file mode 100644
index 0000000..b33bb57
--- /dev/null
+++ b/tests/Makefile
@@ -0,0 +1,15 @@
+CFLAGS = -I../libastro
+
+run-test: test
+	./test > output.actual
+	diff -u output.expected output.actual
+
+test: test.o ../libastro/libastro.a
+	$(CC) $< ../libastro/libastro.a -lm -o $@
+
+test.o: test.c
+
+../libastro/libastro.a: .FORCE
+	make -C ../libastro
+
+.PHONY: .FORCE
diff --git a/tests/output.expected b/tests/output.expected
new file mode 100644
index 0000000..c459432
--- /dev/null
+++ b/tests/output.expected
@@ -0,0 +1,15 @@
+Precession J2000 -> J2015
+RA  +0.0000000000000000 -> +0.0033543649473764
+Dec -0.0000000000000000 -> +0.0014575063671504
+RA  +1.0000000000000000 -> +1.0014450289313686
+Dec -1.0000000000000000 -> -0.9992134064184917
+RA  +2.0000000000000000 -> +5.1478363700948195
+Dec -2.0000000000000000 -> -1.1409819639181964
+RA  +3.0000000000000000 -> +6.1449758360061262
+Dec -3.0000000000000000 -> -0.1401493826857086
+RA  +4.0000000000000000 -> +0.8630426367929015
+Dec -4.0000000000000000 -> +0.8593574646184613
+RA  +5.0000000000000000 -> +4.9986250865903390
+Dec -5.0000000000000000 -> +1.2835978054231891
+RA  +6.0000000000000000 -> +6.0032362075042789
+Dec -6.0000000000000000 -> +0.2845854258344847
diff --git a/tests/test.c b/tests/test.c
new file mode 100644
index 0000000..39c1871
--- /dev/null
+++ b/tests/test.c
@@ -0,0 +1,15 @@
+#include "astro.h"
+
+int main() {
+     printf("Precession J2000 -> J2015\n");
+     double ra, pra;
+     double dec, pdec;
+     for (double i=0.0; i<7.0; i+=1.0) {
+          ra = pra = i;
+          dec = pdec = -i;
+          precess(J2000, J2015, &pra, &pdec);
+          printf("RA  %+.16f -> %+.16f\n", ra, pra);
+          printf("Dec %+.16f -> %+.16f\n", dec, pdec);
+     }
+     return 0;
+}