diff --git a/.idea/workspace.xml b/.idea/workspace.xml
index 890c01c..996d9f0 100644
--- a/.idea/workspace.xml
+++ b/.idea/workspace.xml
@@ -49,6 +49,7 @@
       <updated>1754359057118</updated>
       <workItem from="1754359058201" duration="85000" />
       <workItem from="1754359152820" duration="18000" />
+      <workItem from="1756439414549" duration="601000" />
     </task>
     <servers />
   </component>
diff --git a/mathfunctionize/mathfunctionize.py b/mathfunctionize/mathfunctionize.py
index f3a5af1..4c72883 100644
--- a/mathfunctionize/mathfunctionize.py
+++ b/mathfunctionize/mathfunctionize.py
@@ -58,6 +58,63 @@ def derangements(n):
 def bayes_theorem(priorA, priorB, likelihoodA, likelihoodB):
     return (priorA * likelihoodA) / ((priorA * likelihoodA) + (priorB * likelihoodB))
 
+#complex numbers
+def complex_addition(a, b):
+    if a.find("+") != -1:
+        a1 = float(a[0:a.find("+")])
+        a2 = float(a[a.find("+")+1:a.find("i")])
+    elif a.find("-", 1) != -1:
+        a1 = float(a[0:a.find("-", 1)])
+        a2 = float(a[a.find("-", 1):a.find("i")])
+    elif a.find("i") != -1:
+        a1 = 0
+        a2 = float(a[0:a.find("i")])
+    elif a.find("i") == -1:
+        a1 = float(a)
+        a2 = 0
+    if b.find("+") != -1:
+        b1 = float(b[0:b.find("+")])
+        b2 = float(b[b.find("+")+1:b.find("i")])
+    elif b.find("-", 1) != -1:
+        b1 = float(b[0:b.find("-", 1)])
+        b2 = float(b[b.find("-", 1):b.find("i")])
+    elif b.find("i") != -1:
+        b1 = 0
+        b2 = float(b[0:b.find("i")])
+    elif b.find("i") == -1:
+        b1 = float(b)
+        b2 = 0
+    if a2 + b2 >= 0:
+        return str(a1 + b1) + "+" + str(a2 + b2)+"i"
+    return str(a1 + b1) + str(a2 + b2)+"i"
+def complex_subtraction(a, b):
+    if a.find("+") != -1:
+        a1 = float(a[0:a.find("+")])
+        a2 = float(a[a.find("+")+1:a.find("i")])
+    elif a.find("-", 1) != -1:
+        a1 = float(a[0:a.find("-", 1)])
+        a2 = float(a[a.find("-", 1):a.find("i")])
+    elif a.find("i") != -1:
+        a1 = 0
+        a2 = float(a[0:a.find("i")])
+    elif a.find("i") == -1:
+        a1 = float(a)
+        a2 = 0
+    if b.find("+") != -1:
+        b1 = float(b[0:b.find("+")])
+        b2 = float(b[b.find("+")+1:b.find("i")])
+    elif b.find("-", 1) != -1:
+        b1 = float(b[0:b.find("-", 1)])
+        b2 = float(b[b.find("-", 1):b.find("i")])
+    elif b.find("i") != -1:
+        b1 = 0
+        b2 = float(b[0:b.find("i")])
+    elif b.find("i") == -1:
+        b1 = float(b)
+        b2 = 0
+    if a2 - b2 >= 0:
+        return str(a1 - b1) + "+" + str(a2 - b2)+"i"
+    return str(a1 - b1) + str(a2 - b2)+"i"
 # quantitative analysis
 def localMinimum(arr):
     num = 0