diff --git a/src/harmonicmaps.parset b/src/harmonicmaps.parset
index 52ffe12d0c79e5ca060ef5498a8f717d5fce88da..0735c6b494c53561dbf2d8b1480dd242ad41e401 100644
--- a/src/harmonicmaps.parset
+++ b/src/harmonicmaps.parset
@@ -52,3 +52,6 @@ structuredGrid = true
 lower = -10 -10
 upper =  10  10
 elements = 100 100
+
+# Inverse stereographic projection
+initialIterate = "[2*x[0] / (x[0]*x[0]+x[1]*x[1]+1), 2*x[1] / (x[0]*x[0]+x[1]*x[1]+1), (x[0]*x[0]+x[1]*x[1]-1)/ (x[0]*x[0]+x[1]*x[1]+1)]"