diff --git a/src/state/calculation.wxm b/src/state/calculation.wxm
new file mode 100644
index 0000000000000000000000000000000000000000..d962f03dcb2b0f60ff3c53a118d998280e3455f0
--- /dev/null
+++ b/src/state/calculation.wxm
@@ -0,0 +1,24 @@
+/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
+/* [ Created with wxMaxima version 12.04.0 ] */
+
+/* [wxMaxima: input   start ] */
+d_eq    : 'diff(theta,t) = 1 - theta*(V/L);
+  ode2(d_eq, theta, t)$
+  subst(L/V[ref]*exp(alpha), theta, %)$
+  % / L * V[ref]$
+  ic1(%, t = 0, alpha = alpha[0])$
+d_alpha : log(expand(%));
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+r_eq : 'diff(theta,t) = -(V/L)*theta*log(theta*(V/L));
+  ode2(r_eq,theta,t)$
+  subst(L/V[ref] * exp(alpha),theta,%)$
+  ic1(%,t=0,alpha=alpha[0])$
+  map(exp,-(V/L)*%)$
+  block([logexpand : all], expand(% - log(V/V[ref])))$
+r_alpha : expand(logcontract(%));
+/* [wxMaxima: input   end   ] */
+
+/* Maxima can't load/batch files which end with a comment! */
+"Created with wxMaxima"$