diff --git a/Examples/ff/SasCosineRipple.py b/Examples/ff/SasCosineRipple.py
deleted file mode 100755
index 27715f001d7def1b1b96afffaf12a32f9f8aafde..0000000000000000000000000000000000000000
--- a/Examples/ff/SasCosineRipple.py
+++ /dev/null
@@ -1,20 +0,0 @@
-#!/usr/bin/env python3
-"""
-Plot form factors.
-"""
-import bornagain as ba
-from bornagain import nm, deg
-import bornplot as bp
-
-det = ba.SphericalDetector(bp.simargs['n'], 5*deg, 2.5*deg, 2.5*deg)
-n = 4
-results = []
-for i in range(n):
-    omega = 90*i/(n - 1)
-    title = r'$\omega=%d^\circ$' % omega
-    ff = ba.CosineRipple(25*nm, 10*nm, 8*nm)
-    trafo = ba.RotationZ(omega*deg)
-    data = bp.run_simulation(det, ff, trafo)
-    results.append(bp.Result(data, title))
-
-bp.make_plot(results, det, "ff_CosineRipple")
diff --git a/Examples/ff/sim_demo_1quadrants.py b/Examples/ff/sim_demo_1quadrants.py
deleted file mode 100755
index e5432bf8d5e66f989ee1b2a396c6a09ab34796db..0000000000000000000000000000000000000000
--- a/Examples/ff/sim_demo_1quadrants.py
+++ /dev/null
@@ -1,21 +0,0 @@
-#!/usr/bin/env python3
-"""
-Plot form factor.
-"""
-import bornagain as ba
-from bornagain import nm, deg
-import bornplot as bp
-
-det = ba.SphericalDetector(bp.simargs['n'], 5*deg, 2.5*deg, 2.5*deg)
-pars = [0, 5, 10, 20]
-n = len(pars)
-results = []
-for i in range(n):
-    theta = pars[i]
-    title = r'$\vartheta=%d^\circ$' % theta
-    ff = ba.TruncatedSphere(4.2*nm, 6.1*nm, 0)
-    trafo = ba.RotationY(theta*deg)
-    data = bp.run_simulation(det, ff, trafo)
-    results.append(bp.Result(data, title))
-
-bp.make_plot(results, det, "ff_demo_1quadrants")
diff --git a/Examples/ff/sim_demo_4quadrants.py b/Examples/ff/sim_demo_4quadrants.py
deleted file mode 100755
index 508961656158c42f8022d0baa60a446a89b89447..0000000000000000000000000000000000000000
--- a/Examples/ff/sim_demo_4quadrants.py
+++ /dev/null
@@ -1,21 +0,0 @@
-#!/usr/bin/env python3
-"""
-Plot form factor.
-"""
-import bornagain as ba
-from bornagain import nm, deg
-import bornplot as bp
-
-det = ba.SphericalDetector(bp.simargs['n'], 10*deg, 0, 0)
-pars = [0, 5, 10, 20]
-n = len(pars)
-results = []
-for i in range(n):
-    theta = pars[i]
-    title = r'$\vartheta=%d^\circ$' % theta
-    ff = ba.TruncatedSphere(4.2*nm, 6.1*nm, 0)
-    trafo = ba.RotationY(theta*deg)
-    data = bp.run_simulation(det, ff, trafo)
-    results.append(bp.Result(data, title))
-
-bp.make_plot(results, det, "ff_demo_4quadrants")
diff --git a/Examples/ff/sim_det_box.py b/Examples/ff/sim_det_box.py
deleted file mode 100755
index fb82d28453469459fd15cd19bbb9adad9601d746..0000000000000000000000000000000000000000
--- a/Examples/ff/sim_det_box.py
+++ /dev/null
@@ -1,19 +0,0 @@
-#!/usr/bin/env python3
-"""
-Plot form factors.
-"""
-import bornagain as ba
-from bornagain import nm, deg
-import bornplot as bp
-
-det = ba.SphericalDetector(bp.simargs['n'], 5*deg, 2.5*deg, 2.5*deg)
-n = 3
-results = []
-for i in range(n):
-    siz = [10, 100, 1000][i]
-    title = r'$L_x = %d~nm$' % siz
-    ff = ba.Box(siz, 10, 10)
-    data = bp.run_simulation(det, ff)
-    results.append(bp.Result(data, title))
-
-bp.make_plot(results, det, "ff_det_box")