diff --git a/.gitignore b/.gitignore
index ccc6ca6c..3622b3e5 100644
--- a/.gitignore
+++ b/.gitignore
@@ -175,4 +175,7 @@ src/_version.py
 *.stp
 *.sql
 *.html
-*.png
\ No newline at end of file
+
+*.png
+*.gif
+*.svg
\ No newline at end of file
diff --git a/examples/spherical_tokamak_from_plasma_with_color_animation.py b/examples/spherical_tokamak_from_plasma_with_color_animation.py
new file mode 100644
index 00000000..0785ed10
--- /dev/null
+++ b/examples/spherical_tokamak_from_plasma_with_color_animation.py
@@ -0,0 +1,80 @@
+import paramak
+import imageio
+import cadquery_png_plugin.plugin
+
+# Original radial build values
+original_radial_build = [
+    (paramak.LayerType.GAP, 10),
+    (paramak.LayerType.SOLID, 50),
+    (paramak.LayerType.SOLID, 15),
+    (paramak.LayerType.GAP, 50),
+    (paramak.LayerType.PLASMA, 300),
+    (paramak.LayerType.GAP, 60),
+    (paramak.LayerType.SOLID, 15),
+    (paramak.LayerType.SOLID, 60),
+    (paramak.LayerType.SOLID, 10),
+]
+original_elongation = 2
+original_triangularity = 0.55
+
+# Function to create a reactor with modified radial build
+def create_reactor(radial_build, elongation, triangularity):
+    return paramak.spherical_tokamak_from_plasma(
+        radial_build=radial_build,
+        elongation=elongation,
+        triangularity=triangularity,
+        rotation_angle=180,
+        colors={
+            "layer_1": (0.4, 0.9, 0.4),
+            "layer_2": (0.6, 0.8, 0.6),
+            "plasma": (1., 0.7, 0.8, 0.6),
+            "layer_3": (0.1, 0.1, 0.9),
+            "layer_4": (0.4, 0.4, 0.8),
+            "layer_5": (0.5, 0.5, 0.8),
+        },
+    )
+
+# Function to export reactor to PNG
+def export_reactor_to_png(reactor, file_path):
+    reactor.exportPNG(
+        options={
+            "width": 1280,
+            "height": 1024,
+            "zoom": 1.4,
+        },
+        file_path=file_path
+    )
+
+# Generate reactors with varying radial build values
+frame = 0
+factors = [1.0, 1.25, 1.5, 1.75, 2.0, 1.75, 1.5, 1.25, 1.0]
+for i in range(len(original_radial_build)):
+    layer_type, original_value = original_radial_build[i]
+    for factor in factors:
+        modified_radial_build = original_radial_build.copy()
+        modified_radial_build[i] = (layer_type, original_value * factor)
+        reactor = create_reactor(modified_radial_build, original_elongation, original_triangularity)
+        export_reactor_to_png(reactor, f'spherical_tokamak_frame_{frame:03d}.png')
+        frame += 1
+
+for factor in factors:
+    modified_elongation = original_elongation * factor
+    reactor = create_reactor(original_radial_build, modified_elongation, original_triangularity)
+    export_reactor_to_png(reactor, f'spherical_tokamak_frame_{frame:03d}.png')
+    frame += 1
+
+for modified_triangularity in [0.55, 0.3667, 0.1833, 0.0, -0.1833, -0.3667, -0.55, -0.3667, -0.1833, 0.0, 0.1833, 0.3667, 0.55]:
+    reactor = create_reactor(original_radial_build, original_elongation, modified_triangularity)
+    export_reactor_to_png(reactor, f'spherical_tokamak_frame_{frame:03d}.png')
+    frame += 1
+
+# Convert all the frames PNG files into a GIF
+def create_gif_from_frames(frame_count, output_file):
+    with imageio.get_writer(output_file, mode='I', duration=0.1) as writer:
+        for frame in range(frame_count):
+            filename = f'spherical_tokamak_frame_{frame:03d}.png'
+            image = imageio.imread(filename)
+            writer.append_data(image)
+
+# Create the GIF
+create_gif_from_frames(frame, 'spherical_tokamak_animation.gif')
\ No newline at end of file