diff --git a/test_mex_ifsh.m b/test_mex_ifsh.m
new file mode 100644
index 00000000..bb784aaa
--- /dev/null
+++ b/test_mex_ifsh.m
@@ -0,0 +1,9 @@
+G = computeGeometry(cartGrid([30,30,1]));
+src = [];
+src = addSource(src, 1, 1);
+src = addSource(src, G.cells.num, -1);
+bc = [];
+rock.perm = ones(G.cells.num, 1);
+rock.poro = ones(G.cells.num, 1);
+x = initResSol(G, 0, 0);
+x = mex_ifsh(x, G, rock, bc, src)