diff --git a/toolbox/+otp/+allencahn/f_diffusion.m b/toolbox/+otp/+allencahn/f_diffusion.m
new file mode 100644
index 00000000..bbbb7778
--- /dev/null
+++ b/toolbox/+otp/+allencahn/f_diffusion.m
@@ -0,0 +1,9 @@
+% Process splitting: 
+%   f = f_diffusion + f_reaction
+%   Jacobian = Jacobian_diffusion + Jacobian_reaction
+% 
+function du = f_diffusion(~, u, L, alpha, ~, ~)
+
+du = alpha*L*u;
+
+end
diff --git a/toolbox/+otp/+allencahn/f_linear.m b/toolbox/+otp/+allencahn/f_linear.m
new file mode 100644
index 00000000..52498144
--- /dev/null
+++ b/toolbox/+otp/+allencahn/f_linear.m
@@ -0,0 +1,8 @@
+% Linear-nonlinear splitting: f = f_linear + f_nonlinear
+% f_linear(u) = Jac*u,  f_nonlinear(u) = f(u)-Jac*u
+
+function du = f_linear(~, u, L, alpha, beta, ~)
+
+du = alpha*L*u + beta * (u - 3*u.^3);
+
+end
diff --git a/toolbox/+otp/+allencahn/f_nonlinear.m b/toolbox/+otp/+allencahn/f_nonlinear.m
new file mode 100644
index 00000000..0df35223
--- /dev/null
+++ b/toolbox/+otp/+allencahn/f_nonlinear.m
@@ -0,0 +1,8 @@
+% Linear-nonlinear splitting: f = f_linear + f_nonlinear
+% f_linear(u) = Jac*u,  f_nonlinear(u) = f(u)-Jac*u
+
+function du = f_nonlinear(t, u, ~, ~, beta, forcing)
+
+du = beta * (u - u.^3) - beta * (u - 3*u.^3) + forcing(t);
+
+end
diff --git a/toolbox/+otp/+allencahn/f_reaction.m b/toolbox/+otp/+allencahn/f_reaction.m
new file mode 100644
index 00000000..774d5bbd
--- /dev/null
+++ b/toolbox/+otp/+allencahn/f_reaction.m
@@ -0,0 +1,9 @@
+% Process splitting: 
+%   f = f_diffusion + f_reaction
+%   Jacobian = Jacobian_diffusion + Jacobian_reaction
+% 
+function du = f_reaction(t, u, ~, ~, beta, forcing)
+
+du = beta * (u - u.^3) + forcing(t);
+
+end
diff --git a/toolbox/+otp/+allencahn/jacobian_diffusion.m b/toolbox/+otp/+allencahn/jacobian_diffusion.m
new file mode 100644
index 00000000..c6e84cb7
--- /dev/null
+++ b/toolbox/+otp/+allencahn/jacobian_diffusion.m
@@ -0,0 +1,9 @@
+% Process splitting: 
+%   f = f_diffusion + f_reaction
+%   Jacobian = Jacobian_diffusion + Jacobian_reaction
+%
+function j = jacobian_diffusion(~, ~, L, alpha, ~, ~)
+
+j = alpha*L;
+
+end
diff --git a/toolbox/+otp/+allencahn/jacobian_reaction.m b/toolbox/+otp/+allencahn/jacobian_reaction.m
new file mode 100644
index 00000000..206d3b96
--- /dev/null
+++ b/toolbox/+otp/+allencahn/jacobian_reaction.m
@@ -0,0 +1,9 @@
+% Process splitting: 
+%   f = f_diffusion + f_reaction
+%   Jacobian = Jacobian_diffusion + Jacobian_reaction
+% 
+function j = jacobian_reaction(~, u, L, ~, beta, ~)
+
+j = spdiags(beta * (1 - 3*u.^2), 0, size(L, 1), size(L, 2));
+
+end
diff --git a/toolboxPackaging.prj b/toolboxPackaging.prj
index c9ddefa4..b4dd6978 100644
--- a/toolboxPackaging.prj
+++ b/toolboxPackaging.prj
@@ -1,5 +1,5 @@
 <deployment-project plugin="plugin.toolbox" plugin-version="1.0">
-  <configuration name="ODE Test Problems" target="target.toolbox" target-name="Package Toolbox">
+  <configuration file="/Users/sandu/GIT/GitHub/ODE-Test-Problems/toolboxPackaging.prj" location="/Users/sandu/GIT/GitHub/ODE-Test-Problems" name="toolboxPackaging" target="target.toolbox" target-name="Package Toolbox">
     <param.appname>ODE Test Problems</param.appname>
     <param.authnamewatermark>Steven Roberts, Andrey A. Popov, Arash Sarshar, Adrian Sandu</param.authnamewatermark>
     <param.email />
@@ -14,7 +14,7 @@
     <param.products.version />
     <param.platforms />
     <param.guid>0970ed36-a788-484f-beca-105eecd2562a</param.guid>
-    <param.exclude.filters></param.exclude.filters>
+    <param.exclude.filters />
     <param.exclude.pcodedmfiles>true</param.exclude.pcodedmfiles>
     <param.examples />
     <param.demosxml />
@@ -85,5 +85,59 @@
     <fileset.depfun.included />
     <fileset.depfun.excluded />
     <fileset.package />
+    <build-deliverables>
+      <file location="${PROJECT_ROOT}" name="ODE Test Problems.mltbx" optional="false">/Users/sandu/GIT/GitHub/ODE-Test-Problems/ODE Test Problems.mltbx</file>
+    </build-deliverables>
+    <workflow />
+    <matlab>
+      <root>/Applications/MATLAB_R2024b.app</root>
+      <toolboxes>
+        <toolbox name="matlabcoder" />
+        <toolbox name="embeddedcoder" />
+        <toolbox name="fixedpoint" />
+        <toolbox name="matlabhdlcoder" />
+        <toolbox name="neuralnetwork" />
+      </toolboxes>
+      <toolbox>
+        <matlabcoder>
+          <enabled>true</enabled>
+        </matlabcoder>
+      </toolbox>
+      <toolbox>
+        <embeddedcoder>
+          <enabled>true</enabled>
+        </embeddedcoder>
+      </toolbox>
+      <toolbox>
+        <fixedpoint>
+          <enabled>true</enabled>
+        </fixedpoint>
+      </toolbox>
+      <toolbox>
+        <matlabhdlcoder>
+          <enabled>true</enabled>
+        </matlabhdlcoder>
+      </toolbox>
+      <toolbox>
+        <neuralnetwork>
+          <enabled>true</enabled>
+        </neuralnetwork>
+      </toolbox>
+    </matlab>
+    <platform>
+      <unix>true</unix>
+      <mac>true</mac>
+      <windows>false</windows>
+      <win2k>false</win2k>
+      <winxp>false</winxp>
+      <vista>false</vista>
+      <linux>false</linux>
+      <solaris>false</solaris>
+      <osver>15.5</osver>
+      <os32>true</os32>
+      <os64>false</os64>
+      <arch>maci64</arch>
+      <matlab>true</matlab>
+    </platform>
   </configuration>
 </deployment-project>
\ No newline at end of file