diff --git a/src/math/modf.c b/src/math/modf.c
new file mode 100644
index 0000000..43b8657
--- /dev/null
+++ b/src/math/modf.c
@@ -0,0 +1,90 @@
+#include <stdio.h>
+#include <math.h>
+#include "test.h"
+
+static struct {
+	double x;
+	double yf;
+	double yi;
+} t[] = {
+	0.7, 0x1.6666666666666p-1, 0.0,
+	-0.7, -0x1.6666666666666p-1, -0.0,
+	1.7, 0x1.6666666666666p-1, 1.0,
+	-1.7, -0x1.6666666666666p-1, -1.0,
+	2.0, 0.0, 2.0,
+	-2.0, -0.0, -2.0,
+	-0x1p99, -0.0, -0x1p99,
+	-0.0, -0.0, -0.0,
+	INFINITY, 0.0, INFINITY,
+	-INFINITY, -0.0, -INFINITY,
+	NAN, NAN, NAN,
+};
+
+static struct {
+	float x;
+	float yf;
+	float yi;
+} tf[] = {
+	0.7, 0x1.666666p-1, 0.0,
+	-0.7, -0x1.666666p-1, -0.0,
+	1.7, 0x1.666668p-1, 1.0,
+	-1.7, -0x1.666668p-1, -1.0,
+	2.0, 0.0, 2.0,
+	-2.0, -0.0, -2.0,
+	-0x1p99, -0.0, -0x1p99,
+	-0.0, -0.0, -0.0,
+	INFINITY, 0.0, INFINITY,
+	-INFINITY, -0.0, -INFINITY,
+	NAN, NAN, NAN,
+};
+
+static struct {
+	long double x;
+	long double yf;
+	long double yi;
+} tl[] = {
+	0.7L, 0x1.6666666666666666p-1L, 0.0,
+	-0.7L, -0x1.6666666666666666p-1L, -0.0,
+	1.7L, 0x1.6666666666666668p-1L, 1.0,
+	-1.7L, -0x1.6666666666666668p-1L, -1.0,
+	2.0, 0, 2.0,
+	-2.0, -0.0, -2.0,
+	-0x1p99, -0.0, -0x1p99,
+	-0.0, -0.0, -0.0,
+	INFINITY, 0.0, INFINITY,
+	-INFINITY, -0.0, -INFINITY,
+	NAN, NAN, NAN,
+};
+
+#define eq(a, b) (isnan(a) && isnan(b) || ((a) == (b) && signbit(a) == signbit(b)))
+
+void test_modf()
+{
+	double yf, yi;
+	int i;
+	for (i = 0; i < sizeof t/sizeof *t; i++) {
+		yf = modf(t[i].x, &yi);
+		if (!eq(yf,t[i].yf) || !eq(yi,t[i].yi))
+			error("modf(%a) want %a %a got %a %a\n", t[i].x, t[i].yf, t[i].yi, yf, yi);
+	}
+}
+void test_modff()
+{
+	float yf, yi;
+	int i;
+	for (i = 0; i < sizeof tf/sizeof *tf; i++) {
+		yf = modff(tf[i].x, &yi);
+		if (!eq(yf,tf[i].yf) || !eq(yi,tf[i].yi))
+			error("modff(%a) want %a %a got %a %a\n", tf[i].x, tf[i].yf, tf[i].yi, yf, yi);
+	}
+}
+void test_modfl()
+{
+	long double yf, yi;
+	int i;
+	for (i = 0; i < sizeof tl/sizeof *tl; i++) {
+		yf = modfl(tl[i].x, &yi);
+		if (!eq(yf,tl[i].yf) || !eq(yi,tl[i].yi))
+			error("modfl(%La) want %La %La got %La %La\n", tl[i].x, tl[i].yf, tl[i].yi, yf, yi);
+	}
+}