pimath 0.1.23 → 0.1.25

package/dist/pimath.js

@@ -1,22 +1,22 @@
-var Si = Object.defineProperty;
-var yi = (o) => {
+var Li = Object.defineProperty;
+var gi = (o) => {
   throw TypeError(o);
 };
-var Ri = (o, e, t) => e in o ? Si(o, e, { enumerable: !0, configurable: !0, writable: !0, value: t }) : o[e] = t;
-var a = (o, e, t) => Ri(o, typeof e != "symbol" ? e + "" : e, t), Jt = (o, e, t) => e.has(o) || yi("Cannot " + t);
-var s = (o, e, t) => (Jt(o, e, "read from private field"), t ? t.call(o) : e.get(o)), p = (o, e, t) => e.has(o) ? yi("Cannot add the same private member more than once") : e instanceof WeakSet ? e.add(o) : e.set(o, t), h = (o, e, t, i) => (Jt(o, e, "write to private field"), i ? i.call(o, t) : e.set(o, t), t), q = (o, e, t) => (Jt(o, e, "access private method"), t);
-function zi(o) {
-  const e = vi(o), t = [];
+var Vi = (o, e, t) => e in o ? Li(o, e, { enumerable: !0, configurable: !0, writable: !0, value: t }) : o[e] = t;
+var a = (o, e, t) => Vi(o, typeof e != "symbol" ? e + "" : e, t), _t = (o, e, t) => e.has(o) || gi("Cannot " + t);
+var s = (o, e, t) => (_t(o, e, "read from private field"), t ? t.call(o) : e.get(o)), m = (o, e, t) => e.has(o) ? gi("Cannot add the same private member more than once") : e instanceof WeakSet ? e.add(o) : e.set(o, t), h = (o, e, t, i) => (_t(o, e, "write to private field"), i ? i.call(o, t) : e.set(o, t), t), A = (o, e, t) => (_t(o, e, "access private method"), t);
+function Di(o) {
+  const e = wi(o), t = [];
   let i, n;
   for (; e.length > 0; )
     i = e.shift() ?? 1, n = (e.length > 0 ? e.pop() : +i) ?? 1, t.push([i, n]);
   return t;
 }
-function Li(...o) {
-  const e = fi(...o);
+function Zi(...o) {
+  const e = ui(...o);
   return o.map((t) => t / e);
 }
-function vi(o) {
+function wi(o) {
   const e = Math.abs(o), t = Math.sqrt(e), i = [];
   for (let n = 1; n <= t; n++)
     o % n === 0 && (i.push(n), i.push(e / n));
@@ -24,7 +24,7 @@ function vi(o) {
     return n - r;
   }), [...new Set(i)];
 }
-function fi(...o) {
+function ui(...o) {
   const e = function(n, r) {
     return r === 0 ? n : e(r, n % r);
   };
@@ -39,52 +39,52 @@ function fi(...o) {
     ;
   return Math.abs(t);
 }
-function Di(...o) {
+function Fi(...o) {
   return o.reduce(function(e, t) {
-    return Math.abs(e * t / fi(e, t));
+    return Math.abs(e * t / ui(e, t));
   });
 }
-function Vi(o, e = 3) {
+function ji(o, e = 3) {
   return +o.toFixed(e);
 }
-function Zi(o) {
+function Ui(o) {
   if (Number.isSafeInteger(o) || o.toString().split(".")[0].length < 10)
     return 0;
   throw new Error("Periodic value: Not implemented yet");
 }
-function Fi(o) {
+function Wi(o) {
   const e = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 293, 307, 311, 313, 317, 331, 337, 347, 349, 353, 359, 367, 373, 379, 383, 389, 397, 401, 409, 419, 421, 431, 433, 439, 443, 449, 457, 461, 463, 467, 479, 487, 491, 499, 503, 509, 521, 523, 541, 547, 557, 563, 569, 571, 577, 587, 593, 599, 601, 607, 613, 617, 619, 631, 641, 643, 647, 653, 659, 661, 673, 677, 683, 691, 701, 709, 719, 727, 733, 739, 743, 751, 757, 761, 769, 773, 787, 797, 809, 811, 821, 823, 827, 829, 839, 853, 857, 859, 863, 877, 881, 883, 887, 907, 911, 919, 929, 937, 941, 947, 953, 967, 971, 977, 983, 991, 997, 1009, 1013, 1019, 1021, 1031, 1033, 1039, 1049, 1051, 1061, 1063, 1069, 1087, 1091, 1093, 1097, 1103, 1109, 1117, 1123, 1129, 1151, 1153, 1163, 1171, 1181, 1187, 1193, 1201, 1213, 1217, 1223, 1229, 1231, 1237, 1249, 1259, 1277, 1279, 1283, 1289, 1291, 1297, 1301, 1303, 1307, 1319, 1321, 1327, 1361, 1367, 1373, 1381, 1399, 1409, 1423, 1427, 1429, 1433, 1439, 1447, 1451, 1453, 1459, 1471, 1481, 1483, 1487, 1489, 1493, 1499, 1511, 1523, 1531, 1543, 1549, 1553, 1559, 1567, 1571, 1579, 1583, 1597, 1601, 1607, 1609, 1613, 1619, 1621, 1627, 1637, 1657, 1663, 1667, 1669, 1693, 1697, 1699, 1709, 1721, 1723, 1733, 1741, 1747, 1753, 1759, 1777, 1783, 1787, 1789, 1801, 1811, 1823, 1831, 1847, 1861, 1867, 1871, 1873, 1877, 1879, 1889, 1901, 1907, 1913, 1931, 1933, 1949, 1951, 1973, 1979, 1987, 1993, 1997, 1999, 2003, 2011, 2017, 2027, 2029, 2039, 2053, 2063, 2069, 2081, 2083, 2087, 2089, 2099, 2111, 2113, 2129, 2131, 2137, 2141, 2143, 2153, 2161, 2179, 2203, 2207, 2213, 2221, 2237, 2239, 2243, 2251, 2267, 2269, 2273, 2281, 2287, 2293, 2297, 2309, 2311, 2333, 2339, 2341, 2347, 2351, 2357, 2371, 2377, 2381, 2383, 2389, 2393, 2399, 2411, 2417, 2423, 2437, 2441, 2447, 2459, 2467, 2473, 2477, 2503, 2521, 2531, 2539, 2543, 2549, 2551, 2557, 2579, 2591, 2593, 2609, 2617, 2621, 2633, 2647, 2657, 2659, 2663, 2671, 2677, 2683, 2687, 2689, 2693, 2699, 2707, 2711, 2713, 2719, 2729, 2731, 2741, 2749, 2753, 2767, 2777, 2789, 2791, 2797, 2801, 2803, 2819, 2833, 2837, 2843, 2851, 2857, 2861, 2879, 2887, 2897, 2903, 2909, 2917, 2927, 2939, 2953, 2957, 2963, 2969, 2971, 2999, 3001, 3011, 3019, 3023, 3037, 3041, 3049, 3061, 3067, 3079, 3083, 3089, 3109, 3119, 3121, 3137, 3163, 3167, 3169, 3181, 3187, 3191, 3203, 3209, 3217, 3221, 3229, 3251, 3253, 3257, 3259, 3271, 3299, 3301, 3307, 3313, 3319, 3323, 3329, 3331, 3343, 3347, 3359, 3361, 3371, 3373, 3389, 3391, 3407, 3413, 3433, 3449, 3457, 3461, 3463, 3467, 3469, 3491, 3499, 3511, 3517, 3527, 3529, 3533, 3539, 3541, 3547, 3557, 3559, 3571, 3581, 3583, 3593, 3607, 3613, 3617, 3623, 3631, 3637, 3643, 3659, 3671, 3673, 3677, 3691, 3697, 3701, 3709, 3719, 3727, 3733, 3739, 3761, 3767, 3769, 3779, 3793, 3797, 3803, 3821, 3823, 3833, 3847, 3851, 3853, 3863, 3877, 3881, 3889, 3907, 3911, 3917, 3919, 3923, 3929, 3931, 3943, 3947, 3967, 3989, 4001, 4003, 4007, 4013, 4019, 4021, 4027, 4049, 4051, 4057, 4073, 4079, 4091, 4093, 4099, 4111, 4127, 4129, 4133, 4139, 4153, 4157, 4159, 4177, 4201, 4211, 4217, 4219, 4229, 4231, 4241, 4243, 4253, 4259, 4261, 4271, 4273, 4283, 4289, 4297, 4327, 4337, 4339, 4349, 4357, 4363, 4373, 4391, 4397, 4409, 4421, 4423, 4441, 4447, 4451, 4457, 4463, 4481, 4483, 4493, 4507, 4513, 4517, 4519, 4523, 4547, 4549, 4561, 4567, 4583, 4591, 4597, 4603, 4621, 4637, 4639, 4643, 4649, 4651, 4657, 4663, 4673, 4679, 4691, 4703, 4721, 4723, 4729, 4733, 4751, 4759, 4783, 4787, 4789, 4793, 4799, 4801, 4813, 4817, 4831, 4861, 4871, 4877, 4889, 4903, 4909, 4919, 4931, 4933, 4937, 4943, 4951, 4957, 4967, 4969, 4973, 4987, 4993, 4999, 5003, 5009, 5011, 5021, 5023, 5039, 5051, 5059, 5077, 5081, 5087, 5099, 5101, 5107, 5113, 5119, 5147, 5153, 5167, 5171, 5179, 5189, 5197, 5209, 5227, 5231, 5233, 5237, 5261, 5273, 5279, 5281, 5297, 5303, 5309, 5323, 5333, 5347, 5351, 5381, 5387, 5393, 5399, 5407, 5413, 5417, 5419, 5431, 5437, 5441, 5443, 5449, 5471, 5477, 5479, 5483, 5501, 5503, 5507, 5519, 5521, 5527, 5531, 5557, 5563, 5569, 5573, 5581, 5591, 5623, 5639, 5641, 5647, 5651, 5653, 5657, 5659, 5669, 5683, 5689, 5693, 5701, 5711, 5717, 5737, 5741, 5743, 5749, 5779, 5783, 5791, 5801, 5807, 5813, 5821, 5827, 5839, 5843, 5849, 5851, 5857, 5861, 5867, 5869, 5879, 5881, 5897, 5903, 5923, 5927, 5939, 5953, 5981, 5987, 6007, 6011, 6029, 6037, 6043, 6047, 6053, 6067, 6073, 6079, 6089, 6091, 6101, 6113, 6121, 6131, 6133, 6143, 6151, 6163, 6173, 6197, 6199, 6203, 6211, 6217, 6221, 6229, 6247, 6257, 6263, 6269, 6271, 6277, 6287, 6299, 6301, 6311, 6317, 6323, 6329, 6337, 6343, 6353, 6359, 6361, 6367, 6373, 6379, 6389, 6397, 6421, 6427, 6449, 6451, 6469, 6473, 6481, 6491, 6521, 6529, 6547, 6551, 6553, 6563, 6569, 6571, 6577, 6581, 6599, 6607, 6619, 6637, 6653, 6659, 6661, 6673, 6679, 6689, 6691, 6701, 6703, 6709, 6719, 6733, 6737, 6761, 6763, 6779, 6781, 6791, 6793, 6803, 6823, 6827, 6829, 6833, 6841, 6857, 6863, 6869, 6871, 6883, 6899, 6907, 6911, 6917, 6947, 6949, 6959, 6961, 6967, 6971, 6977, 6983, 6991, 6997, 7001, 7013, 7019, 7027, 7039, 7043, 7057, 7069, 7079, 7103, 7109, 7121, 7127, 7129, 7151, 7159, 7177, 7187, 7193, 7207, 7211, 7213, 7219, 7229, 7237, 7243, 7247, 7253, 7283, 7297, 7307, 7309, 7321, 7331, 7333, 7349, 7351, 7369, 7393, 7411, 7417, 7433, 7451, 7457, 7459, 7477, 7481, 7487, 7489, 7499, 7507, 7517, 7523, 7529, 7537, 7541, 7547, 7549, 7559, 7561, 7573, 7577, 7583, 7589, 7591, 7603, 7607, 7621, 7639, 7643, 7649, 7669, 7673, 7681, 7687, 7691, 7699, 7703, 7717, 7723, 7727, 7741, 7753, 7757, 7759, 7789, 7793, 7817, 7823, 7829, 7841, 7853, 7867, 7873, 7877, 7879, 7883, 7901, 7907, 7919, 7927, 7933, 7937, 7949, 7951, 7963, 7993, 8009, 8011, 8017, 8039, 8053, 8059, 8069, 8081, 8087, 8089, 8093, 8101, 8111, 8117, 8123, 8147, 8161, 8167, 8171, 8179, 8191, 8209, 8219, 8221, 8231, 8233, 8237, 8243, 8263, 8269, 8273, 8287, 8291, 8293, 8297, 8311, 8317, 8329, 8353, 8363, 8369, 8377, 8387, 8389, 8419, 8423, 8429, 8431, 8443, 8447, 8461, 8467, 8501, 8513, 8521, 8527, 8537, 8539, 8543, 8563, 8573, 8581, 8597, 8599, 8609, 8623, 8627, 8629, 8641, 8647, 8663, 8669, 8677, 8681, 8689, 8693, 8699, 8707, 8713, 8719, 8731, 8737, 8741, 8747, 8753, 8761, 8779, 8783, 8803, 8807, 8819, 8821, 8831, 8837, 8839, 8849, 8861, 8863, 8867, 8887, 8893, 8923, 8929, 8933, 8941, 8951, 8963, 8969, 8971, 8999, 9001, 9007, 9011, 9013, 9029, 9041, 9043, 9049, 9059, 9067, 9091, 9103, 9109, 9127, 9133, 9137, 9151, 9157, 9161, 9173, 9181, 9187, 9199, 9203, 9209, 9221, 9227, 9239, 9241, 9257, 9277, 9281, 9283, 9293, 9311, 9319, 9323, 9337, 9341, 9343, 9349, 9371, 9377, 9391, 9397, 9403, 9413, 9419, 9421, 9431, 9433, 9437, 9439, 9461, 9463, 9467, 9473, 9479, 9491, 9497, 9511, 9521, 9533, 9539, 9547, 9551, 9587, 9601, 9613, 9619, 9623, 9629, 9631, 9643, 9649, 9661, 9677, 9679, 9689, 9697, 9719, 9721, 9733, 9739, 9743, 9749, 9767, 9769, 9781, 9787, 9791, 9803, 9811, 9817, 9829, 9833, 9839, 9851, 9857, 9859, 9871, 9883, 9887, 9901, 9907, 9923, 9929, 9931, 9941, 9949, 9967, 9973];
   return o === void 0 ? e : e.slice(0, Math.min(e.length, o));
 }
-function ji(o, e) {
+function Gi(o, e) {
   const t = [], i = e === !0 ? +o : o ** 2;
   for (let n = 0; n <= o; n++)
     for (let r = 0; r <= o; r++)
       n ** 2 + r ** 2 === i && t.push([n, r, o]);
   return t;
 }
-function Ui(o, e = 2) {
+function Hi(o, e = 2) {
   return +`${Math.round(+`${o}e${e}`)}e-${e}`;
 }
-const G = {
-  decompose: zi,
-  dividers: vi,
-  divideNumbersByGCD: Li,
-  gcd: fi,
-  lcm: Di,
-  numberCorrection: Vi,
-  periodic: Zi,
-  primes: Fi,
-  pythagoreanTripletsWithTarget: ji,
-  round: Ui
+const U = {
+  decompose: Di,
+  dividers: wi,
+  divideNumbersByGCD: Zi,
+  gcd: ui,
+  lcm: Fi,
+  numberCorrection: ji,
+  periodic: Ui,
+  primes: Wi,
+  pythagoreanTripletsWithTarget: Gi,
+  round: Hi
 };
-var ct, b, w, ze;
-const $ = class $ {
+var ut, x, v, Re;
+const P = class P {
   constructor(e, t) {
-    p(this, ct, !1);
-    p(this, b, 1);
-    p(this, w, 1);
-    p(this, ze, "frac");
+    m(this, ut, !1);
+    m(this, x, 1);
+    m(this, v, 1);
+    m(this, Re, "frac");
     // ------------------------------------------
     /**
      * Parse the value to get the numerator and denominator
@@ -94,7 +94,7 @@ const $ = class $ {
     a(this, "parse", (e, t) => {
       let i;
       if (e === "")
-        return h(this, w, 0), h(this, b, 1), this;
+        return h(this, v, 0), h(this, x, 1), this;
       switch (typeof e) {
         case "string":
           if (i = e.split("/"), i.length > 2)
@@ -103,37 +103,36 @@ const $ = class $ {
             throw new Error(`The given value is not a valid fraction (${e})`);
           if (i.length === 1)
             return this.parse(+i[0]);
-          i.length === 2 ? i[1] === "0" ? (h(this, w, NaN), h(this, b, 1)) : (h(this, w, +i[0]), h(this, b, +i[1])) : (h(this, w, NaN), h(this, b, 1));
+          i.length === 2 ? i[1] === "0" ? (h(this, v, NaN), h(this, x, 1)) : (h(this, v, +i[0]), h(this, x, +i[1])) : (h(this, v, NaN), h(this, x, 1));
           break;
         case "number":
           if (Number.isSafeInteger(e))
-            h(this, w, +e), t === void 0 || !Number.isSafeInteger(t) ? h(this, b, 1) : h(this, b, +t);
+            h(this, v, +e), t === void 0 || !Number.isSafeInteger(t) ? h(this, x, 1) : h(this, x, +t);
           else {
             const [, n] = e.toString().split("."), r = n ? n.length : 0;
-            t === void 0 ? (h(this, w, e * Math.pow(10, r)), h(this, b, Math.pow(10, r))) : Number.isSafeInteger(t) && (h(this, w, e * Math.pow(10, r) - Math.floor(e * Math.pow(10, r - t))), this.denominator = Math.pow(10, r) - Math.pow(10, r - t)), this.reduce();
+            t === void 0 ? (h(this, v, e * Math.pow(10, r)), h(this, x, Math.pow(10, r))) : Number.isSafeInteger(t) && (h(this, v, e * Math.pow(10, r) - Math.floor(e * Math.pow(10, r - t))), this.denominator = Math.pow(10, r) - Math.pow(10, r - t)), this.reduce();
           }
           break;
         case "object":
-          e instanceof $ && (h(this, w, +e.numerator), h(this, b, +e.denominator));
+          e instanceof P && (h(this, v, +e.numerator), h(this, x, +e.denominator));
           break;
       }
       return this;
     });
     a(this, "clone", () => {
-      const e = new $();
-      return e.numerator = +s(this, w), e.denominator = +s(this, b), e;
+      const e = new P();
+      return e.numerator = +s(this, v), e.denominator = +s(this, x), e;
     });
-    a(this, "abs", () => (h(this, w, Math.abs(s(this, w))), h(this, b, Math.abs(s(this, b))), this));
+    a(this, "abs", () => (h(this, v, Math.abs(s(this, v))), h(this, x, Math.abs(s(this, x))), this));
     a(this, "add", (e) => {
-      if (e instanceof $) {
-        const t = s(this, w), i = s(this, b);
-        h(this, w, t * e.denominator + e.numerator * i), h(this, b, i * e.denominator);
+      if (e instanceof P) {
+        const t = s(this, v), i = s(this, x);
+        h(this, v, t * e.denominator + e.numerator * i), h(this, x, i * e.denominator);
       } else
-        return this.add(new $(e));
+        return this.add(new P(e));
       return this.reduce();
     });
-    a(this, "amplify", (e) => (Number.isSafeInteger(e) && (h(this, w, s(this, w) * e), h(this, b, s(this, b) * e)), this));
-    // TODO: The rest of the functions are not used or unnecessary ?
+    a(this, "amplify", (e) => (Number.isSafeInteger(e) && (h(this, v, s(this, v) * e), h(this, x, s(this, x) * e)), this));
     /**
      * Simple function to determine if it's a fraction
      */
@@ -147,7 +146,7 @@ const $ = class $ {
     a(this, "compare", (e, t) => {
       t === void 0 && (t = "=");
       let i;
-      switch (e instanceof $ ? i = e.clone() : i = new $(e), t) {
+      switch (e instanceof P ? i = e.clone() : i = new P(e), t) {
         case ">":
           return this.value > i.value;
         case ">=":
@@ -169,66 +168,66 @@ const $ = class $ {
       }
     });
     a(this, "divide", (e) => {
-      const t = new $(e);
+      const t = new P(e);
       if (t.numerator === 0)
-        return new $().infinite();
-      const i = +s(this, w), n = +s(this, b);
-      return h(this, w, i * t.denominator), h(this, b, n * t.numerator), this.reduce();
+        return new P().infinite();
+      const i = +s(this, v), n = +s(this, x);
+      return h(this, v, i * t.denominator), h(this, x, n * t.numerator), this.reduce();
     });
-    a(this, "infinite", () => (h(this, w, 1 / 0), h(this, b, 1), this));
-    a(this, "invalid", () => (h(this, w, NaN), h(this, b, 1), this));
+    a(this, "infinite", () => (h(this, v, 1 / 0), h(this, x, 1), this));
+    a(this, "invalid", () => (h(this, v, NaN), h(this, x, 1), this));
     a(this, "inverse", () => {
-      const e = +s(this, w);
-      return h(this, w, +s(this, b)), h(this, b, e), this;
+      const e = +s(this, v);
+      return h(this, v, +s(this, x)), h(this, x, e), this;
     });
-    a(this, "isApproximative", () => s(this, ct) || s(this, w).toString().length >= 15 && s(this, b).toString().length >= 15);
+    a(this, "isApproximative", () => s(this, ut) || s(this, v).toString().length >= 15 && s(this, x).toString().length >= 15);
     a(this, "isEqual", (e) => this.compare(e, "="));
     a(this, "isEven", () => this.isRelative() && this.value % 2 === 0);
     a(this, "isExact", () => !this.isApproximative());
     a(this, "isFinite", () => !this.isInfinity() && !this.isNaN());
     a(this, "isGeq", (e) => this.compare(e, ">="));
     a(this, "isGreater", (e) => this.compare(e, ">"));
-    a(this, "isInfinity", () => Math.abs(s(this, w)) === 1 / 0);
-    a(this, "isInverted", (e) => this.isEqual(new $().one().divide(e.clone())));
+    a(this, "isInfinity", () => Math.abs(s(this, v)) === 1 / 0);
+    a(this, "isInverted", (e) => this.isEqual(new P().one().divide(e.clone())));
     a(this, "isLeq", (e) => this.compare(e, "<="));
     /* Compare shortcuts */
     a(this, "isLesser", (e) => this.compare(e, "<"));
-    a(this, "isNaN", () => isNaN(s(this, w)));
+    a(this, "isNaN", () => isNaN(s(this, v)));
     a(this, "isNatural", () => this.isRelative() && this.isPositive());
     a(this, "isNegative", () => this.sign() === -1);
-    a(this, "isNegativeOne", () => s(this, w) === -1 && s(this, b) === 1);
+    a(this, "isNegativeOne", () => s(this, v) === -1 && s(this, x) === 1);
     a(this, "isNotEqual", (e) => this.compare(e, "<>"));
     // ------------------------------------------
-    a(this, "isNotZero", () => s(this, w) !== 0);
+    a(this, "isNotZero", () => s(this, v) !== 0);
     a(this, "isOdd", () => this.isRelative() && this.value % 2 === 1);
-    a(this, "isOne", () => s(this, w) === 1 && s(this, b) === 1);
+    a(this, "isOne", () => s(this, v) === 1 && s(this, x) === 1);
     a(this, "isOpposite", (e) => this.isEqual(e.clone().opposite()));
     a(this, "isPositive", () => this.sign() === 1);
     a(this, "isRational", () => !this.isRelative());
-    a(this, "isReduced", () => Math.abs(G.gcd(s(this, w), s(this, b))) === 1);
+    a(this, "isReduced", () => Math.abs(U.gcd(s(this, v), s(this, x))) === 1);
     a(this, "isRelative", () => this.clone().reduce().denominator === 1);
-    a(this, "isSquare", () => Math.sqrt(s(this, w)) % 1 === 0 && Math.sqrt(s(this, b)) % 1 === 0);
+    a(this, "isSquare", () => Math.sqrt(s(this, v)) % 1 === 0 && Math.sqrt(s(this, x)) % 1 === 0);
     a(this, "isStrictlyNegative", () => this.value < 0);
     a(this, "isStrictlyPositive", () => this.value > 0);
     // Mathematical operations specific to fractions
-    a(this, "isZero", () => s(this, w) === 0);
+    a(this, "isZero", () => s(this, v) === 0);
     a(this, "multiply", (e) => {
-      const t = new $(e);
-      return h(this, w, s(this, w) * t.numerator), h(this, b, s(this, b) * t.denominator), this.reduce();
+      const t = new P(e);
+      return h(this, v, s(this, v) * t.numerator), h(this, x, s(this, x) * t.denominator), this.reduce();
     });
-    a(this, "one", () => (h(this, w, 1), h(this, b, 1), this));
-    a(this, "opposite", () => (h(this, w, -s(this, w)), this));
+    a(this, "one", () => (h(this, v, 1), h(this, x, 1), this));
+    a(this, "opposite", () => (h(this, v, -s(this, v)), this));
     a(this, "pow", (e) => {
-      if (e instanceof $)
+      if (e instanceof P)
         return this.pow(e.value);
       this.reduce(), e < 0 && this.inverse();
-      const t = Math.floor(Math.pow(s(this, w), Math.abs(e))), i = Math.floor(Math.pow(s(this, b), Math.abs(e)));
-      return t ** Math.abs(e) === s(this, w) && i ** Math.abs(e) === s(this, b) ? (h(this, w, s(this, w) ** Math.abs(e)), h(this, b, s(this, b) ** Math.abs(e))) : (h(this, w, s(this, w) ** Math.abs(e)), h(this, b, s(this, b) ** Math.abs(e))), this;
+      const t = Math.floor(Math.pow(s(this, v), Math.abs(e))), i = Math.floor(Math.pow(s(this, x), Math.abs(e)));
+      return t ** Math.abs(e) === s(this, v) && i ** Math.abs(e) === s(this, x) ? (h(this, v, s(this, v) ** Math.abs(e)), h(this, x, s(this, x) ** Math.abs(e))) : (h(this, v, s(this, v) ** Math.abs(e)), h(this, x, s(this, x) ** Math.abs(e))), this;
     });
     // ------------------------------------------
     a(this, "reduce", () => {
-      const e = G.gcd(s(this, w), s(this, b));
-      return h(this, w, s(this, w) / e), h(this, b, s(this, b) / e), s(this, b) < 0 && (h(this, b, -s(this, b)), h(this, w, -s(this, w))), this;
+      const e = U.gcd(s(this, v), s(this, x));
+      return h(this, v, s(this, v) / e), h(this, x, s(this, x) / e), s(this, x) < 0 && (h(this, x, -s(this, x)), h(this, v, -s(this, v))), this;
     });
     a(this, "root", (e) => {
       if (e === 0)
@@ -239,111 +238,111 @@ const $ = class $ {
         throw new Error("The root of a negative number must be odd.");
       const t = this.sign();
       this.abs(), this.reduce();
-      const i = Math.floor(Math.pow(s(this, w), Math.abs(1 / e))), n = Math.floor(Math.pow(s(this, b), Math.abs(1 / e)));
-      return h(this, w, Math.pow(s(this, w), Math.abs(1 / e))), h(this, b, Math.pow(s(this, b), Math.abs(1 / e))), (i !== s(this, w) || n !== s(this, b)) && (h(this, w, s(this, w) / s(this, b)), h(this, b, 1), h(this, ct, !0)), this.multiply(t), this;
+      const i = Math.floor(Math.pow(s(this, v), Math.abs(1 / e))), n = Math.floor(Math.pow(s(this, x), Math.abs(1 / e)));
+      return h(this, v, Math.pow(s(this, v), Math.abs(1 / e))), h(this, x, Math.pow(s(this, x), Math.abs(1 / e))), (i !== s(this, v) || n !== s(this, x)) && (h(this, v, s(this, v) / s(this, x)), h(this, x, 1), h(this, ut, !0)), this.multiply(t), this;
     });
     // ------------------------------------------
     // Getter and setter
-    a(this, "sign", () => s(this, w) * s(this, b) >= 0 ? 1 : -1);
+    a(this, "sign", () => s(this, v) * s(this, x) >= 0 ? 1 : -1);
     a(this, "sqrt", () => this.root(2));
-    a(this, "subtract", (e) => e instanceof $ ? this.add(e.clone().opposite()) : this.add(-e));
-    a(this, "zero", () => (h(this, w, 0), h(this, b, 1), this));
+    a(this, "subtract", (e) => e instanceof P ? this.add(e.clone().opposite()) : this.add(-e));
+    a(this, "zero", () => (h(this, v, 0), h(this, x, 1), this));
     return e !== void 0 && this.parse(e, t), this;
   }
   get denominator() {
-    return s(this, b);
+    return s(this, x);
   }
   set denominator(e) {
-    h(this, b, e);
+    h(this, x, e);
   }
   get dfrac() {
-    return h(this, ze, "dfrac"), this;
+    return h(this, Re, "dfrac"), this;
   }
   get display() {
-    return this.isExact() ? s(this, b) === 1 ? `${s(this, w)}` : `${s(this, w)}/${s(this, b)}` : this.value.toFixed(3);
+    return this.isExact() ? s(this, x) === 1 ? `${s(this, v)}` : `${s(this, v)}/${s(this, x)}` : this.value.toFixed(3);
   }
   get frac() {
-    return h(this, ze, "frac"), this;
+    return h(this, Re, "frac"), this;
   }
   // ------------------------------------------
   get numerator() {
-    return s(this, w);
+    return s(this, v);
   }
   set numerator(e) {
-    h(this, w, e);
+    h(this, v, e);
   }
   // Display getter
   get tex() {
-    return this.isInfinity() ? `${this.sign() === 1 ? "+" : "-"}\\infty` : this.isExact() ? s(this, b) === 1 ? `${s(this, w)}` : s(this, w) < 0 ? `-\\${s(this, ze)}{ ${-s(this, w)} }{ ${s(this, b)} }` : `\\${s(this, ze)}{ ${s(this, w)} }{ ${s(this, b)} }` : this.value.toFixed(3);
+    return this.isInfinity() ? `${this.sign() === 1 ? "+" : "-"}\\infty` : this.isExact() ? s(this, x) === 1 ? `${s(this, v)}` : s(this, v) < 0 ? `-\\${s(this, Re)}{ ${-s(this, v)} }{ ${s(this, x)} }` : `\\${s(this, Re)}{ ${s(this, v)} }{ ${s(this, x)} }` : this.value.toFixed(3);
   }
   get texWithSign() {
     return this.isPositive() ? `+${this.tex}` : this.tex;
   }
   get tfrac() {
-    return h(this, ze, "tfrac"), this;
+    return h(this, Re, "tfrac"), this;
   }
   get value() {
-    const e = s(this, w) / s(this, b);
+    const e = s(this, v) / s(this, x);
     return e === 0 ? 0 : e;
   }
 };
-ct = new WeakMap(), b = new WeakMap(), w = new WeakMap(), ze = new WeakMap(), a($, "average", (...e) => {
-  const t = new $().zero();
+ut = new WeakMap(), x = new WeakMap(), v = new WeakMap(), Re = new WeakMap(), a(P, "average", (...e) => {
+  const t = new P().zero();
   for (const i of e)
     t.add(i);
   return t.divide(e.length), t;
-}), a($, "max", (...e) => {
-  let t = new $(e[0]);
+}), a(P, "max", (...e) => {
+  let t = new P(e[0]);
   for (const i of e) {
-    const n = new $(i);
+    const n = new P(i);
     n.isGreater(t) && (t = n.clone());
   }
   return t;
-}), a($, "min", (...e) => {
-  let t = new $(e[0]);
+}), a(P, "min", (...e) => {
+  let t = new P(e[0]);
   for (const i of e) {
-    const n = new $(i);
+    const n = new P(i);
     n.isLesser(t) && (t = n.clone());
   }
   return t;
-}), a($, "sort", (e, t) => {
-  const n = e.map((r) => r instanceof $ ? r : new $(r)).sort((r, l) => r.value - l.value);
+}), a(P, "sort", (e, t) => {
+  const n = e.map((r) => r instanceof P ? r : new P(r)).sort((r, l) => r.value - l.value);
   return t && n.reverse(), n;
 }), // ------------------------------------------
 // Compare functions
-a($, "unique", (e) => {
+a(P, "unique", (e) => {
   const t = {}, i = [];
   return e.forEach((n) => {
-    n instanceof $ || (n = new $(n)), t[n.clone().reduce().tex] || (i.push(n.clone()), t[n.tex] = !0);
+    n instanceof P || (n = new P(n)), t[n.clone().reduce().tex] || (i.push(n.clone()), t[n.tex] = !0);
   }), i;
-}), a($, "xMultiply", (...e) => {
-  const t = new $();
+}), a(P, "xMultiply", (...e) => {
+  const t = new P();
   for (const i of e) {
-    const n = new $(i);
+    const n = new P(i);
     t.numerator = t.numerator * n.numerator, t.denominator = t.denominator * n.denominator;
   }
   return t;
 });
-let u = $;
-var Z, ee, ie, Xe;
+let u = P;
+var Z, ie, ne, Ye;
 class vt {
   constructor(...e) {
-    p(this, Z);
-    p(this, ee);
-    p(this, ie);
-    p(this, Xe);
+    m(this, Z);
+    m(this, ie);
+    m(this, ne);
+    m(this, Ye);
     // ------------------------------------------
     // Creation / parsing functions
     // ------------------------------------------
-    a(this, "parse", (e, t, i) => (h(this, ie, i ?? 1), h(this, ee, t ?? 2), h(this, Z, e), s(this, ee) % 2 === 0 && s(this, Z) < 0 && h(this, Xe, !1), this));
+    a(this, "parse", (e, t, i) => (h(this, ne, i ?? 1), h(this, ie, t ?? 2), h(this, Z, e), s(this, ie) % 2 === 0 && s(this, Z) < 0 && h(this, Ye, !1), this));
     // ------------------------------------------
     // Mathematical operations
     // ------------------------------------------
     a(this, "reduce", () => {
-      let e = Math.floor(Math.pow(s(this, Z), 1 / s(this, ee)));
+      let e = Math.floor(Math.pow(s(this, Z), 1 / s(this, ie)));
       for (; e > 1; ) {
-        if (s(this, Z) % Math.pow(e, s(this, ee)) === 0) {
-          h(this, ie, s(this, ie) * e), h(this, Z, s(this, Z) / Math.pow(e, s(this, ee))), e = Math.floor(Math.pow(s(this, Z), 1 / s(this, ee)));
+        if (s(this, Z) % Math.pow(e, s(this, ie)) === 0) {
+          h(this, ne, s(this, ne) * e), h(this, Z, s(this, Z) / Math.pow(e, s(this, ie))), e = Math.floor(Math.pow(s(this, Z), 1 / s(this, ie)));
           continue;
         }
         e--;
@@ -354,8 +353,8 @@ class vt {
     // ------------------------------------------
     // Help functions
     // ------------------------------------------
-    a(this, "hasRadical", () => !(s(this, Z) === 1 || s(this, Z) === 0 || !s(this, Xe)));
-    h(this, Z, 1), h(this, ie, 1), h(this, ee, 2), h(this, Xe, !0), e.length > 0 && this.parse(e[0], e[1], e[2]);
+    a(this, "hasRadical", () => !(s(this, Z) === 1 || s(this, Z) === 0 || !s(this, Ye)));
+    h(this, Z, 1), h(this, ne, 1), h(this, ie, 2), h(this, Ye, !0), e.length > 0 && this.parse(e[0], e[1], e[2]);
   }
   // ------------------------------------------
   // Getter and setter
@@ -367,274 +366,315 @@ class vt {
     h(this, Z, e);
   }
   get nth() {
-    return s(this, ee);
+    return s(this, ie);
   }
   set nth(e) {
-    Number.isSafeInteger(e) && e >= 2 ? h(this, ee, e) : (console.log("Error setting the nth root"), h(this, ee, 2));
+    Number.isSafeInteger(e) && e >= 2 ? h(this, ie, e) : (console.log("Error setting the nth root"), h(this, ie, 2));
   }
   get coefficient() {
-    return s(this, ie);
+    return s(this, ne);
   }
   set coefficient(e) {
-    h(this, ie, e);
+    h(this, ne, e);
   }
   get tex() {
     let e;
-    return s(this, ie) === 1 ? e = "" : s(this, ie) === -1 ? e = "-" : e = s(this, ie).toString(), s(this, Z) === 1 ? `${s(this, ie)}` : s(this, ee) === 2 ? `${e}\\sqrt{${s(this, Z)}}` : `${e}\\sqrt[${s(this, ee)}]{${s(this, Z)}}`;
+    return s(this, ne) === 1 ? e = "" : s(this, ne) === -1 ? e = "-" : e = s(this, ne).toString(), s(this, Z) === 1 ? `${s(this, ne)}` : s(this, ie) === 2 ? `${e}\\sqrt{${s(this, Z)}}` : `${e}\\sqrt[${s(this, ie)}]{${s(this, Z)}}`;
   }
   get display() {
     let e;
-    return s(this, ie) === 1 ? e = "" : s(this, ie) === -1 ? e = "-" : e = s(this, ie).toString(), s(this, Z) === 1 ? `${s(this, ie)}` : s(this, ee) === 2 ? `${e}sqrt{${s(this, Z)}}` : `${e}root(${s(this, ee)}){${s(this, Z)}}`;
+    return s(this, ne) === 1 ? e = "" : s(this, ne) === -1 ? e = "-" : e = s(this, ne).toString(), s(this, Z) === 1 ? `${s(this, ne)}` : s(this, ie) === 2 ? `${e}sqrt{${s(this, Z)}}` : `${e}root(${s(this, ie)}){${s(this, Z)}}`;
   }
   get value() {
-    return s(this, ie) * Math.pow(s(this, Z), 1 / s(this, ee));
+    return s(this, ne) * Math.pow(s(this, Z), 1 / s(this, ie));
   }
 }
-Z = new WeakMap(), ee = new WeakMap(), ie = new WeakMap(), Xe = new WeakMap();
-var fe, je, P, ot, ye, bi, ti, xi, Ei, Ni;
-const mi = class mi {
+Z = new WeakMap(), ie = new WeakMap(), ne = new WeakMap(), Ye = new WeakMap();
+var K, Fe, C, Ge, ge, vi, bi, xi, Ei, Ni, Ai, Ti;
+const pi = class pi {
   constructor(e, t, i = "x") {
-    p(this, P);
-    p(this, fe);
-    p(this, je);
-    if (h(this, je, i), Object.hasOwn(e, "moveLeft")) {
+    m(this, C);
+    m(this, K);
+    m(this, Fe);
+    if (h(this, Fe, i), Object.hasOwn(e, "moveLeft")) {
       const n = e;
-      h(this, fe, n.left.clone().subtract(n.right));
+      h(this, K, n.left.clone().subtract(n.right));
     } else
-      h(this, fe, e.clone().subtract(t ?? 0));
+      h(this, K, e.clone().subtract(t ?? 0));
   }
   solve() {
-    if (s(this, fe).degree().isOne())
-      return q(this, P, xi).call(this);
-    if (s(this, fe).degree().value === 2)
-      return q(this, P, Ei).call(this);
-    const e = q(this, P, bi).call(this);
-    if (e.length > 0)
-      return e;
-    if (s(this, fe).degree().value === 3)
-      return q(this, P, ti).call(this);
-    throw new Error("The equation degree is too high.");
+    const e = s(this, K).degree().value;
+    if (e === 0)
+      return [];
+    if (e === 1)
+      return A(this, C, Ni).call(this);
+    if (e === 2)
+      return A(this, C, Ai).call(this);
+    const t = A(this, C, xi).call(this);
+    return t.length > 0 ? t : A(this, C, vi).call(this);
   }
   solveAsCardan() {
-    if (s(this, fe).degree().value !== 3)
+    if (s(this, K).degree().value !== 3)
       throw new Error("The equation is not cubic.");
-    return q(this, P, ti).call(this);
+    return A(this, C, Ei).call(this);
   }
 };
-fe = new WeakMap(), je = new WeakMap(), P = new WeakSet(), ot = function(e, t) {
+K = new WeakMap(), Fe = new WeakMap(), C = new WeakSet(), Ge = function(e, t) {
   return {
-    variable: s(this, je),
+    variable: s(this, Fe),
     exact: !1,
     value: +e.toFixed(10),
     tex: (t == null ? void 0 : t.tex) ?? "",
     display: (t == null ? void 0 : t.display) ?? ""
   };
-}, ye = function(e) {
+}, ge = function(e) {
   if (e instanceof u && e.isApproximative())
-    return q(this, P, ot).call(this, e.value);
+    return A(this, C, Ge).call(this, e.value);
   const t = new u(e);
   return {
-    variable: s(this, je),
+    variable: s(this, Fe),
     exact: t,
     value: t.value,
     tex: t.tex,
     display: t.display
   };
-}, bi = function() {
-  let e = s(this, fe).clone(), t = [];
+}, vi = function() {
+  const e = [];
+  s(this, K).degree().value;
+  const t = s(this, K).getCoefficients().map((p) => p.value), [i, ...n] = s(this, K).getCoefficients(), r = 1 + Math.max(...n.map((p) => p.value / i.value)), c = 2 * r / 100, f = [];
+  for (let p = -r; p <= r; p += c) {
+    const w = U.numberCorrection(p);
+    f.push(
+      {
+        x: w,
+        fx: s(this, K).evaluate(w, !0)
+      }
+    );
+  }
+  f.sort((p, w) => p.x - w.x);
+  const d = [];
+  return f.forEach((p, w) => {
+    w > 0 && (p.fx === 0 ? d.push([p.x, p.x]) : f[w - 1].fx * p.fx < 0 && d.push([
+      f[w - 1].x,
+      p.x
+    ]));
+  }), d.forEach((p) => {
+    const [w, V] = p;
+    if (w === V)
+      e.push(A(this, C, ge).call(this, w));
+    else {
+      const ce = A(this, C, bi).call(this, t, w, V);
+      ce !== null && e.push(A(this, C, Ge).call(this, ce));
+    }
+  }), e;
+}, bi = function(e, t, i, n = 1e-10) {
+  let r = s(this, K).evaluate(t, !0), l = s(this, K).evaluate(i, !0);
+  if (r * l > 0)
+    return console.log("Pas de racine dans l'intervalle donné"), null;
+  let c;
+  for (; (i - t) / 2 > n; ) {
+    c = (t + i) / 2;
+    const f = s(this, K).evaluate(c, !0);
+    if (f === 0)
+      return c;
+    r * f < 0 ? (i = c, l = f) : (t = c, r = f);
+  }
+  return (t + i) / 2;
+}, xi = function() {
+  let e = s(this, K).clone(), t = [];
   const i = e.lcmDenominator();
   i !== 1 && e.multiply(i);
   const n = e.monomByDegree().coefficient;
   let r = e.monomByDegree(0).coefficient;
   for (; r.isZero(); )
-    t.length === 0 && t.push(q(this, P, ye).call(this, 0)), e = e.divide("x"), r = e.monomByDegree(0).coefficient;
-  const l = G.dividers(n.value), c = G.dividers(r.value);
-  for (const f of l)
-    for (const g of c) {
-      const N = new u(g, f);
-      e.evaluate(N).isZero() && !t.find((V) => V.value === N.value) && t.push(q(this, P, ye).call(this, N)), N.opposite(), e.evaluate(N).isZero() && !t.find((V) => V.value === N.value) && t.push(q(this, P, ye).call(this, N));
+    t.length === 0 && t.push(A(this, C, ge).call(this, 0)), e = e.divide("x"), r = e.monomByDegree(0).coefficient;
+  const l = U.dividers(n.value), c = U.dividers(r.value);
+  for (const d of l)
+    for (const p of c) {
+      const w = new u(p, d);
+      e.evaluate(w).isZero() && !t.find((V) => V.value === w.value) && t.push(A(this, C, ge).call(this, w)), w.opposite(), e.evaluate(w).isZero() && !t.find((V) => V.value === w.value) && t.push(A(this, C, ge).call(this, w));
     }
-  for (const f of t) {
-    if (f.exact !== !1 && f.exact.isZero())
+  for (const d of t) {
+    if (d.exact !== !1 && d.exact.isZero())
       continue;
-    const g = s(this, fe).clone().parse("x", f.exact.denominator, -f.exact.numerator);
-    for (; e.isDividableBy(g); )
-      e = e.divide(g);
+    const p = s(this, K).clone().parse("x", d.exact.denominator, -d.exact.numerator);
+    for (; e.isDividableBy(p); )
+      e = e.divide(p);
   }
   if (e.degree().isZero())
-    return t.sort((f, g) => f.value - g.value);
+    return t.sort((d, p) => d.value - p.value);
   if (e.degree().value > 3)
     return [];
-  const d = new mi(e, e.clone().parse("0"), s(this, je));
-  return t = t.concat(d.solve()), t.sort((f, g) => f.value - g.value);
-}, ti = function() {
-  const e = s(this, fe), t = e.monomByDegree(3).coefficient, i = e.monomByDegree(2).coefficient, n = e.monomByDegree(1).coefficient, r = e.monomByDegree(0).coefficient, l = i.clone().divide(t), c = n.clone().divide(t), d = r.clone().divide(t), f = c.clone().subtract(l.clone().pow(2).divide(3)), g = d.clone().subtract(l.clone().multiply(c).divide(3)).add(l.clone().pow(3).multiply(2).divide(27)), N = g.clone().opposite(), V = f.clone().opposite().pow(3).divide(27), pe = N.clone().pow(2).subtract(V.clone().multiply(4)).opposite();
-  if (pe.isNegative()) {
-    const he = g.clone().opposite().add(pe.clone().opposite().sqrt()).divide(2).root(3), ae = g.clone().opposite().subtract(pe.clone().opposite().sqrt()).divide(2).root(3), ge = he.clone().add(ae).subtract(l.clone().divide(3));
-    return [q(this, P, ye).call(this, ge)];
-  }
-  if (pe.isZero()) {
-    const he = g.clone().opposite().divide(2).root(3), ae = he.clone().opposite().subtract(l.clone().divide(3)), ge = he.clone().multiply(2).subtract(l.clone().divide(3));
-    return ae.isEqual(ge) ? [q(this, P, ye).call(this, ae)] : [
-      q(this, P, ye).call(this, ge),
-      q(this, P, ye).call(this, ae)
-    ].sort((me, ce) => me.value - ce.value);
-  }
-  if (pe.isPositive()) {
-    const he = [], ae = f.value, ge = g.value, me = l.value;
-    for (let ce = 0; ce < 3; ce++)
-      he.push(2 * Math.sqrt(-ae / 3) * Math.cos(Math.acos(3 * ge / (2 * ae) * Math.sqrt(-3 / ae)) / 3 + 2 * Math.PI * ce / 3) - me / 3);
-    return he.map((ce) => q(this, P, ot).call(this, ce)).sort((ce, nt) => ce.value - nt.value);
+  const f = new pi(e, e.clone().parse("0"), s(this, Fe));
+  return t = t.concat(f.solve()), t.sort((d, p) => d.value - p.value);
+}, Ei = function() {
+  const e = s(this, K), t = e.monomByDegree(3).coefficient, i = e.monomByDegree(2).coefficient, n = e.monomByDegree(1).coefficient, r = e.monomByDegree(0).coefficient, l = i.clone().divide(t), c = n.clone().divide(t), f = r.clone().divide(t), d = c.clone().subtract(l.clone().pow(2).divide(3)), p = f.clone().subtract(l.clone().multiply(c).divide(3)).add(l.clone().pow(3).multiply(2).divide(27)), w = p.clone().opposite(), V = d.clone().opposite().pow(3).divide(27), ce = w.clone().pow(2).subtract(V.clone().multiply(4)).opposite();
+  if (ce.isNegative()) {
+    const ue = p.clone().opposite().add(ce.clone().opposite().sqrt()).divide(2).root(3), fe = p.clone().opposite().subtract(ce.clone().opposite().sqrt()).divide(2).root(3), ye = ue.clone().add(fe).subtract(l.clone().divide(3));
+    return [A(this, C, ge).call(this, ye)];
+  }
+  if (ce.isZero()) {
+    const ue = p.clone().opposite().divide(2).root(3), fe = ue.clone().opposite().subtract(l.clone().divide(3)), ye = ue.clone().multiply(2).subtract(l.clone().divide(3));
+    return fe.isEqual(ye) ? [A(this, C, ge).call(this, fe)] : [
+      A(this, C, ge).call(this, ye),
+      A(this, C, ge).call(this, fe)
+    ].sort((me, de) => me.value - de.value);
+  }
+  if (ce.isPositive()) {
+    const ue = [], fe = d.value, ye = p.value, me = l.value;
+    for (let de = 0; de < 3; de++)
+      ue.push(2 * Math.sqrt(-fe / 3) * Math.cos(Math.acos(3 * ye / (2 * fe) * Math.sqrt(-3 / fe)) / 3 + 2 * Math.PI * de / 3) - me / 3);
+    return ue.map((de) => A(this, C, Ge).call(this, de)).sort((de, rt) => de.value - rt.value);
   }
   return [];
-}, xi = function() {
-  const e = s(this, fe).monomByDegree(0).coefficient.clone().opposite().divide(s(this, fe).monomByDegree(1).coefficient);
+}, Ni = function() {
+  const [e, t] = s(this, K).getCoefficients(), i = t.opposite().divide(e);
   return [
-    q(this, P, ye).call(this, e)
+    A(this, C, ge).call(this, i)
   ];
-}, Ei = function() {
-  const e = s(this, fe), t = e.monomByDegree(2).coefficient, i = e.monomByDegree(1).coefficient, n = e.monomByDegree(0).coefficient, r = i.clone().pow(2).subtract(t.clone().multiply(n).multiply(4));
+}, Ai = function() {
+  const e = s(this, K), [t, i, n] = e.getCoefficients(), r = i.clone().pow(2).subtract(t.clone().multiply(n).multiply(4));
   if (r.isNegative())
     return [];
   if (r.isSquare()) {
-    const l = r.sqrt(), c = i.clone().opposite().add(l).divide(t.clone().multiply(2)), d = i.clone().opposite().subtract(l).divide(t.clone().multiply(2));
-    return l.isZero() ? [q(this, P, ye).call(this, c)] : [
-      q(this, P, ye).call(this, c),
-      q(this, P, ye).call(this, d)
-    ].sort((f, g) => f.value - g.value);
-  }
-  return q(this, P, Ni).call(this, t, i, r);
-}, Ni = function(e, t, i) {
-  const n = G.dividers(i.value).filter((me) => Math.sqrt(me) % 1 === 0).map((me) => Math.sqrt(me)).pop() ?? 1, r = G.gcd(2 * e.value, t.value, n) * (e.isNegative() ? -1 : 1), l = t.clone().divide(r).opposite(), c = e.clone().divide(r).multiply(2), d = i.clone().divide(n ** 2), f = Math.abs(n / r), g = n === 1 ? "-" : `-${f} `, N = n === 1 ? "+" : `+${f} `;
-  function V(me, ce, nt, Kt) {
-    return `\\frac{ ${ce} ${nt}\\sqrt{ ${Kt} } }{ ${me} }`;
-  }
-  function pe(me, ce, nt, Kt) {
-    return `(${ce}${nt}sqrt(${Kt}))/${me}`;
-  }
-  const he = i.value ** 0.5, ae = (-t.value - he) / (2 * e.value), ge = (-t.value + he) / (2 * e.value);
+    const l = r.sqrt(), c = i.clone().opposite().add(l).divide(t.clone().multiply(2)), f = i.clone().opposite().subtract(l).divide(t.clone().multiply(2));
+    return l.isZero() ? [A(this, C, ge).call(this, c)] : [
+      A(this, C, ge).call(this, c),
+      A(this, C, ge).call(this, f)
+    ].sort((d, p) => d.value - p.value);
+  }
+  return A(this, C, Ti).call(this, t, i, r);
+}, Ti = function(e, t, i) {
+  const n = U.dividers(i.value).filter((me) => Math.sqrt(me) % 1 === 0).map((me) => Math.sqrt(me)).pop() ?? 1, r = U.gcd(2 * e.value, t.value, n) * (e.isNegative() ? -1 : 1), l = t.clone().divide(r).opposite(), c = e.clone().divide(r).multiply(2), f = i.clone().divide(n ** 2), d = Math.abs(n / r), p = n === 1 ? "-" : `-${d} `, w = n === 1 ? "+" : `+${d} `;
+  function V(me, de, rt, Jt) {
+    return `\\frac{ ${de} ${rt}\\sqrt{ ${Jt} } }{ ${me} }`;
+  }
+  function ce(me, de, rt, Jt) {
+    return `(${de}${rt}sqrt(${Jt}))/${me}`;
+  }
+  const ue = i.value ** 0.5, fe = (-t.value - ue) / (2 * e.value), ye = (-t.value + ue) / (2 * e.value);
   return [
-    q(this, P, ot).call(this, ae, {
-      tex: V(c.tex, l.tex, g.toString(), d.tex),
-      display: pe(c.display, l.display, g.toString(), d.display)
+    A(this, C, Ge).call(this, fe, {
+      tex: V(c.tex, l.tex, p.toString(), f.tex),
+      display: ce(c.display, l.display, p.toString(), f.display)
     }),
-    q(this, P, ot).call(this, ge, {
-      tex: V(c.tex, l.tex, N.toString(), d.tex),
-      display: pe(c.display, l.display, N.toString(), d.display)
+    A(this, C, Ge).call(this, ye, {
+      tex: V(c.tex, l.tex, w.toString(), f.tex),
+      display: ce(c.display, l.display, w.toString(), f.display)
     })
-  ].sort((me, ce) => me.value - ce.value);
+  ].sort((me, de) => me.value - de.value);
 };
-let Nt = mi;
-var Gi = Object.defineProperty, Ti = (o) => {
+let At = pi;
+var Xi = Object.defineProperty, Oi = (o) => {
   throw TypeError(o);
-}, Wi = (o, e, t) => e in o ? Gi(o, e, { enumerable: !0, configurable: !0, writable: !0, value: t }) : o[e] = t, _t = (o, e, t) => Wi(o, typeof e != "symbol" ? e + "" : e, t), Ai = (o, e, t) => e.has(o) || Ti("Cannot " + t), _ = (o, e, t) => (Ai(o, e, "read from private field"), t ? t.call(o) : e.get(o)), rt = (o, e, t) => e.has(o) ? Ti("Cannot add the same private member more than once") : e instanceof WeakSet ? e.add(o) : e.set(o, t), xe = (o, e, t, i) => (Ai(o, e, "write to private field"), e.set(o, t), t);
-const di = {
+}, Yi = (o, e, t) => e in o ? Xi(o, e, { enumerable: !0, configurable: !0, writable: !0, value: t }) : o[e] = t, ei = (o, e, t) => Yi(o, typeof e != "symbol" ? e + "" : e, t), qi = (o, e, t) => e.has(o) || Oi("Cannot " + t), te = (o, e, t) => (qi(o, e, "read from private field"), t ? t.call(o) : e.get(o)), ot = (o, e, t) => e.has(o) ? Oi("Cannot add the same private member more than once") : e instanceof WeakSet ? e.add(o) : e.set(o, t), Ee = (o, e, t, i) => (qi(o, e, "write to private field"), e.set(o, t), t);
+const fi = {
   pi: Math.PI,
   e: Math.exp(1)
 };
-var m = /* @__PURE__ */ ((o) => (o.VARIABLE = "variable", o.COEFFICIENT = "coefficient", o.OPERATION = "operation", o.CONSTANT = "constant", o.FUNCTION = "function", o.FUNCTION_ARGUMENT = "function-argument", o.MONOM = "monom", o.LEFT_PARENTHESIS = "(", o.RIGHT_PARENTHESIS = ")", o))(m || {}), Fe = /* @__PURE__ */ ((o) => (o.EXPRESSION = "expression", o.POLYNOM = "polynom", o.SET = "set", o.NUMERIC = "numeric", o))(Fe || {});
-function Hi(o, e) {
+var g = /* @__PURE__ */ ((o) => (o.VARIABLE = "variable", o.COEFFICIENT = "coefficient", o.OPERATION = "operation", o.CONSTANT = "constant", o.FUNCTION = "function", o.FUNCTION_ARGUMENT = "function-argument", o.MONOM = "monom", o.LEFT_PARENTHESIS = "(", o.RIGHT_PARENTHESIS = ")", o))(g || {}), Ze = /* @__PURE__ */ ((o) => (o.EXPRESSION = "expression", o.POLYNOM = "polynom", o.SET = "set", o.NUMERIC = "numeric", o))(Ze || {});
+function Qi(o, e) {
   if (o.length <= 1)
     return o;
-  const t = Object.keys(e).filter((N) => e[N].type === m.FUNCTION).map((N) => N);
-  t.sort((N, V) => V.length - N.length);
-  const i = new RegExp(`^(${t.join("|")})\\(`), n = Object.keys(di);
-  n.sort((N, V) => V.length - N.length);
+  const t = Object.keys(e).filter((w) => e[w].type === g.FUNCTION).map((w) => w);
+  t.sort((w, V) => V.length - w.length);
+  const i = new RegExp(`^(${t.join("|")})\\(`), n = Object.keys(fi);
+  n.sort((w, V) => V.length - w.length);
   const r = new RegExp(`^(${n.join("|")})`), l = /^(\d+(\.\d+)?)/;
-  let c = "", d, f, g;
+  let c = "", f, d, p;
   for (; o.length > 0; ) {
-    if (d = f, g = void 0, t.length > 0 && i.exec(o)) {
-      const N = t.find((V) => o.startsWith(V));
-      N && (g = N + "(", o = o.slice(N.length + 1), f = m.FUNCTION);
+    if (f = d, p = void 0, t.length > 0 && i.exec(o)) {
+      const w = t.find((V) => o.startsWith(V));
+      w && (p = w + "(", o = o.slice(w.length + 1), d = g.FUNCTION);
     } else if (n.length > 0 && r.exec(o)) {
-      const N = n.find((V) => o.startsWith(V));
-      N && (g = N, o = o.slice(N.length), f = m.CONSTANT);
+      const w = n.find((V) => o.startsWith(V));
+      w && (p = w, o = o.slice(w.length), d = g.CONSTANT);
     } else if (l.exec(o)) {
-      const N = l.exec(o);
-      N && (g = N[0], o = o.slice(N[0].length), f = m.COEFFICIENT);
+      const w = l.exec(o);
+      w && (p = w[0], o = o.slice(w[0].length), d = g.COEFFICIENT);
     } else
-      switch (g = o[0], o = o.slice(1), g) {
+      switch (p = o[0], o = o.slice(1), p) {
         case "(":
-          f = m.LEFT_PARENTHESIS;
+          d = g.LEFT_PARENTHESIS;
           break;
         case ")":
-          f = m.RIGHT_PARENTHESIS;
+          d = g.RIGHT_PARENTHESIS;
           break;
         case ",":
-          f = m.FUNCTION_ARGUMENT;
+          d = g.FUNCTION_ARGUMENT;
           break;
         case "+":
         case "-":
         case "*":
         case "/":
         case "^":
-          f = m.OPERATION;
+          d = g.OPERATION;
           break;
         default:
-          f = m.VARIABLE;
+          d = g.VARIABLE;
       }
-    if (g === void 0 || f === void 0)
+    if (p === void 0 || d === void 0)
       throw new Error("The token is undefined");
-    c += Xi(d, f), c += g;
+    c += Ki(f, d), c += p;
   }
   return c;
 }
-function Xi(o, e) {
-  return o === void 0 || o === m.OPERATION || e === m.OPERATION || o === m.LEFT_PARENTHESIS || o === m.FUNCTION || o === m.FUNCTION_ARGUMENT || e === m.RIGHT_PARENTHESIS || e === m.FUNCTION_ARGUMENT ? "" : "*";
+function Ki(o, e) {
+  return o === void 0 || o === g.OPERATION || e === g.OPERATION || o === g.LEFT_PARENTHESIS || o === g.FUNCTION || o === g.FUNCTION_ARGUMENT || e === g.RIGHT_PARENTHESIS || e === g.FUNCTION_ARGUMENT ? "" : "*";
 }
-const Yi = {
-  "^": { precedence: 4, associative: "right", type: m.OPERATION },
-  "*": { precedence: 3, associative: "left", type: m.OPERATION },
-  "/": { precedence: 3, associative: "left", type: m.OPERATION },
-  "+": { precedence: 2, associative: "left", type: m.OPERATION },
-  "-": { precedence: 2, associative: "left", type: m.OPERATION }
-}, Qi = {
-  "^": { precedence: 4, associative: "right", type: m.OPERATION },
-  "*": { precedence: 3, associative: "left", type: m.OPERATION },
-  "/": { precedence: 3, associative: "left", type: m.OPERATION },
-  "+": { precedence: 2, associative: "left", type: m.OPERATION },
-  "-": { precedence: 2, associative: "left", type: m.OPERATION },
-  "%": { precedence: 3, associative: "right", type: m.OPERATION },
-  sin: { precedence: 4, associative: "right", type: m.FUNCTION },
-  cos: { precedence: 4, associative: "right", type: m.FUNCTION },
-  tan: { precedence: 4, associative: "right", type: m.FUNCTION },
-  sqrt: { precedence: 4, associative: "right", type: m.FUNCTION },
-  nthrt: { precedence: 4, associative: "right", type: m.FUNCTION },
-  ",": { precedence: 2, associative: "left", type: m.FUNCTION_ARGUMENT }
-}, Ki = {
-  "^": { precedence: 4, associative: "right", type: m.OPERATION },
-  "*": { precedence: 3, associative: "left", type: m.OPERATION },
-  "/": { precedence: 3, associative: "left", type: m.OPERATION },
-  "+": { precedence: 2, associative: "left", type: m.OPERATION },
-  "-": { precedence: 2, associative: "left", type: m.OPERATION },
-  "%": { precedence: 3, associative: "right", type: m.OPERATION },
-  sin: { precedence: 4, associative: "right", type: m.FUNCTION },
-  cos: { precedence: 4, associative: "right", type: m.FUNCTION },
-  tan: { precedence: 4, associative: "right", type: m.FUNCTION },
-  sqrt: { precedence: 4, associative: "right", type: m.FUNCTION },
-  nthrt: { precedence: 4, associative: "right", type: m.FUNCTION },
-  ln: { precedence: 4, associative: "right", type: m.FUNCTION },
-  log: { precedence: 4, associative: "right", type: m.FUNCTION }
-}, Ji = {
-  "&": { precedence: 3, associative: "left", type: m.OPERATION },
-  "|": { precedence: 3, associative: "left", type: m.OPERATION },
-  "!": { precedence: 4, associative: "right", type: m.OPERATION },
-  "-": { precedence: 2, associative: "left", type: m.OPERATION }
+const Ji = {
+  "^": { precedence: 4, associative: "right", type: g.OPERATION },
+  "*": { precedence: 3, associative: "left", type: g.OPERATION },
+  "/": { precedence: 3, associative: "left", type: g.OPERATION },
+  "+": { precedence: 2, associative: "left", type: g.OPERATION },
+  "-": { precedence: 2, associative: "left", type: g.OPERATION }
+}, _i = {
+  "^": { precedence: 4, associative: "right", type: g.OPERATION },
+  "*": { precedence: 3, associative: "left", type: g.OPERATION },
+  "/": { precedence: 3, associative: "left", type: g.OPERATION },
+  "+": { precedence: 2, associative: "left", type: g.OPERATION },
+  "-": { precedence: 2, associative: "left", type: g.OPERATION },
+  "%": { precedence: 3, associative: "right", type: g.OPERATION },
+  sin: { precedence: 4, associative: "right", type: g.FUNCTION },
+  cos: { precedence: 4, associative: "right", type: g.FUNCTION },
+  tan: { precedence: 4, associative: "right", type: g.FUNCTION },
+  sqrt: { precedence: 4, associative: "right", type: g.FUNCTION },
+  nthrt: { precedence: 4, associative: "right", type: g.FUNCTION },
+  ",": { precedence: 2, associative: "left", type: g.FUNCTION_ARGUMENT }
+}, es = {
+  "^": { precedence: 4, associative: "right", type: g.OPERATION },
+  "*": { precedence: 3, associative: "left", type: g.OPERATION },
+  "/": { precedence: 3, associative: "left", type: g.OPERATION },
+  "+": { precedence: 2, associative: "left", type: g.OPERATION },
+  "-": { precedence: 2, associative: "left", type: g.OPERATION },
+  "%": { precedence: 3, associative: "right", type: g.OPERATION },
+  sin: { precedence: 4, associative: "right", type: g.FUNCTION },
+  cos: { precedence: 4, associative: "right", type: g.FUNCTION },
+  tan: { precedence: 4, associative: "right", type: g.FUNCTION },
+  sqrt: { precedence: 4, associative: "right", type: g.FUNCTION },
+  nthrt: { precedence: 4, associative: "right", type: g.FUNCTION },
+  ln: { precedence: 4, associative: "right", type: g.FUNCTION },
+  log: { precedence: 4, associative: "right", type: g.FUNCTION }
+}, ts = {
+  "&": { precedence: 3, associative: "left", type: g.OPERATION },
+  "|": { precedence: 3, associative: "left", type: g.OPERATION },
+  "!": { precedence: 4, associative: "right", type: g.OPERATION },
+  "-": { precedence: 2, associative: "left", type: g.OPERATION }
 };
-var We, ht, te, bt, Ve;
-class Qt {
+var He, at, se, bt, Ve;
+class Kt {
   constructor(e) {
-    rt(this, We), rt(this, ht, []), rt(this, te, {}), rt(this, bt, []), rt(this, Ve), xe(this, We, typeof e > "u" ? Fe.POLYNOM : e), this.tokenConfigInitialization();
+    ot(this, He), ot(this, at, []), ot(this, se, {}), ot(this, bt, []), ot(this, Ve), Ee(this, He, typeof e > "u" ? Ze.POLYNOM : e), this.tokenConfigInitialization();
   }
   // Getter
   get rpn() {
-    return _(this, ht);
+    return te(this, at);
   }
   get rpnToken() {
-    return _(this, ht).map((e) => e.token);
+    return te(this, at).map((e) => e.token);
   }
   tokenConfigInitialization() {
-    return _(this, We) === Fe.SET ? (xe(this, te, Ji), xe(this, Ve, !1)) : _(this, We) === Fe.NUMERIC ? (xe(this, te, Ki), xe(this, Ve, !0)) : _(this, We) === Fe.EXPRESSION ? (xe(this, te, Qi), xe(this, Ve, !0)) : (xe(this, te, Yi), xe(this, Ve, !0)), xe(this, bt, Object.keys(_(this, te)).sort((e, t) => t.length - e.length)), _(this, te);
+    return te(this, He) === Ze.SET ? (Ee(this, se, ts), Ee(this, Ve, !1)) : te(this, He) === Ze.NUMERIC ? (Ee(this, se, es), Ee(this, Ve, !0)) : te(this, He) === Ze.EXPRESSION ? (Ee(this, se, _i), Ee(this, Ve, !0)) : (Ee(this, se, Ji), Ee(this, Ve, !0)), Ee(this, bt, Object.keys(te(this, se)).sort((e, t) => t.length - e.length)), te(this, se);
   }
   /**
    * Get the next token to analyse.
@@ -644,31 +684,31 @@ class Qt {
   NextToken(e, t) {
     let i, n;
     if (i = "", n = void 0, e[t] === "(")
-      i = "(", n = m.LEFT_PARENTHESIS;
+      i = "(", n = g.LEFT_PARENTHESIS;
     else if (e[t] === ")")
-      i = ")", n = m.RIGHT_PARENTHESIS;
+      i = ")", n = g.RIGHT_PARENTHESIS;
     else if (e[t] === ",")
-      i = ",", n = m.FUNCTION_ARGUMENT;
+      i = ",", n = g.FUNCTION_ARGUMENT;
     else {
-      for (const r of _(this, bt))
+      for (const r of te(this, bt))
         if (e.substring(t, t + r.length) === r) {
-          i += r, n = _(this, te)[r].type;
+          i += r, n = te(this, se)[r].type;
           break;
         }
-      for (const r in di)
+      for (const r in fi)
         if (e.substring(t, t + r.length) === r) {
-          i += r, n = m.CONSTANT;
+          i += r, n = g.CONSTANT;
           break;
         }
       if (i === "")
         if (/[0-9.]/.exec(e[t])) {
           const r = /^([0-9.]+)/.exec(e.substring(t));
-          i = r ? r[0] : "", n = m.COEFFICIENT;
+          i = r ? r[0] : "", n = g.COEFFICIENT;
         } else if (/[a-zA-Z]/.exec(e[t])) {
           const r = /^([a-zA-Z])/.exec(e.substring(t));
-          i = r ? r[0] : "", n = m.VARIABLE;
+          i = r ? r[0] : "", n = g.VARIABLE;
         } else
-          console.log("Unidentified token", e[t], e, t), i = e[t], n = m.MONOM;
+          console.log("Unidentified token", e[t], e, t), i = e[t], n = g.MONOM;
     }
     if (n === void 0)
       throw new Error(`Token type is undefined for token ${i}`);
@@ -683,56 +723,56 @@ class Qt {
   parse(e, t) {
     const i = [], n = [];
     let r = "", l = 0, c;
-    (t ?? _(this, Ve)) && (e = Hi(e, _(this, te)));
-    const d = 50;
-    let f = 50, g;
+    (t ?? te(this, Ve)) && (e = Qi(e, te(this, se)));
+    const f = 50;
+    let d = 50, p;
     for (; l < e.length; ) {
-      if (f--, f === 0) {
+      if (d--, d === 0) {
         console.log("SECURITY LEVEL 1 EXIT");
         break;
       }
       switch ([r, l, c] = this.NextToken(e, l), c) {
-        case m.MONOM:
-        case m.COEFFICIENT:
-        case m.VARIABLE:
-        case m.CONSTANT:
+        case g.MONOM:
+        case g.COEFFICIENT:
+        case g.VARIABLE:
+        case g.CONSTANT:
           i.push({
             token: r,
             tokenType: c
           });
           break;
-        case m.OPERATION:
+        case g.OPERATION:
           if (n.length > 0) {
-            let N = n[n.length - 1];
-            for (g = +d; N.token in _(this, te) && //either o1 is left-associative and its precedence is less than or equal to that of o2,
-            (_(this, te)[r].associative === "left" && _(this, te)[r].precedence <= _(this, te)[N.token].precedence || //or o1 is right associative, and has precedence less than that of o2,
-            _(this, te)[r].associative === "right" && _(this, te)[r].precedence < _(this, te)[N.token].precedence); ) {
-              if (g--, g === 0) {
+            let w = n[n.length - 1];
+            for (p = +f; w.token in te(this, se) && //either o1 is left-associative and its precedence is less than or equal to that of o2,
+            (te(this, se)[r].associative === "left" && te(this, se)[r].precedence <= te(this, se)[w.token].precedence || //or o1 is right associative, and has precedence less than that of o2,
+            te(this, se)[r].associative === "right" && te(this, se)[r].precedence < te(this, se)[w.token].precedence); ) {
+              if (p--, p === 0) {
                 console.log("SECURITY LEVEL 2 OPERATION EXIT");
                 break;
               }
-              if (i.push(n.pop() ?? { token: "", tokenType: m.OPERATION }), n.length === 0)
+              if (i.push(n.pop() ?? { token: "", tokenType: g.OPERATION }), n.length === 0)
                 break;
-              N = n[n.length - 1];
+              w = n[n.length - 1];
             }
           }
           n.push({ token: r, tokenType: c });
           break;
-        case m.FUNCTION_ARGUMENT:
-          for (g = +d; n[n.length - 1].token !== "(" && n.length > 0; ) {
-            if (g--, g === 0) {
+        case g.FUNCTION_ARGUMENT:
+          for (p = +f; n[n.length - 1].token !== "(" && n.length > 0; ) {
+            if (p--, p === 0) {
               console.log("SECURITY LEVEL 2 FUNCTION ARGUMENT EXIT");
               break;
             }
             i.push(n.pop() ?? { token: r, tokenType: c });
           }
           break;
-        case m.LEFT_PARENTHESIS:
-          n.push({ token: r, tokenType: c }), e[l] === "-" && i.push({ token: "0", tokenType: m.COEFFICIENT });
+        case g.LEFT_PARENTHESIS:
+          n.push({ token: r, tokenType: c }), e[l] === "-" && i.push({ token: "0", tokenType: g.COEFFICIENT });
           break;
-        case m.RIGHT_PARENTHESIS:
-          for (g = +d; n[n.length - 1].token !== "(" && n.length > 1; ) {
-            if (g--, g === 0) {
+        case g.RIGHT_PARENTHESIS:
+          for (p = +f; n[n.length - 1].token !== "(" && n.length > 1; ) {
+            if (p--, p === 0) {
               console.log("SECURITY LEVEL 2 CLOSING PARENTHESIS EXIT");
               break;
             }
@@ -740,22 +780,22 @@ class Qt {
           }
           n.pop();
           break;
-        case m.FUNCTION:
+        case g.FUNCTION:
           n.push({ token: r, tokenType: c });
           break;
         default:
           throw new Error(`Token type ${r} is not handled`);
       }
     }
-    return xe(this, ht, i.concat(n.reverse())), this;
+    return Ee(this, at, i.concat(n.reverse())), this;
   }
 }
-We = /* @__PURE__ */ new WeakMap(), ht = /* @__PURE__ */ new WeakMap(), te = /* @__PURE__ */ new WeakMap(), bt = /* @__PURE__ */ new WeakMap(), Ve = /* @__PURE__ */ new WeakMap();
-class _i {
+He = /* @__PURE__ */ new WeakMap(), at = /* @__PURE__ */ new WeakMap(), se = /* @__PURE__ */ new WeakMap(), bt = /* @__PURE__ */ new WeakMap(), Ve = /* @__PURE__ */ new WeakMap();
+class is {
   constructor(e, t) {
-    _t(this, "_rpn"), _t(this, "_expression"), _t(this, "_isValid"), this._expression = e;
+    ei(this, "_rpn"), ei(this, "_expression"), ei(this, "_isValid"), this._expression = e;
     try {
-      this._rpn = new Qt(Fe.NUMERIC).parse(e, t).rpn;
+      this._rpn = new Kt(Ze.NUMERIC).parse(e, t).rpn;
     } catch {
       throw this._rpn = null, this._isValid = !1, new Error(`There was a problem parsing: ${e}`);
     }
@@ -784,7 +824,7 @@ class _i {
       return this._isValid = !1, 0;
     this._isValid = !0;
     for (const i of this._rpn)
-      if (i.tokenType === m.COEFFICIENT)
+      if (i.tokenType === g.COEFFICIENT)
         if (!isNaN(+i.token))
           t.push(+i.token);
         else {
@@ -793,11 +833,11 @@ class _i {
             throw this._isValid = !1, new Error("This coefficient is not a fraction");
           t.push(+n[0] / +n[1]);
         }
-      else if (i.tokenType === m.VARIABLE && e !== void 0)
+      else if (i.tokenType === g.VARIABLE && e !== void 0)
         Object.hasOwn(e, i.token) && t.push(+e[i.token]);
-      else if (i.tokenType === m.CONSTANT)
-        t.push(di[i.token]);
-      else if (i.tokenType === m.OPERATION) {
+      else if (i.tokenType === g.CONSTANT)
+        t.push(fi[i.token]);
+      else if (i.tokenType === g.OPERATION) {
         if (i.token === "*") {
           const n = t.pop(), r = t.pop();
           if (r === void 0 || n === void 0)
@@ -824,7 +864,7 @@ class _i {
             throw this._isValid = !1, new Error(`The base value ${r ?? "a"} or exponent ${n ?? "b"} are not defined`);
           t.push(Math.pow(r, n));
         }
-      } else if (i.tokenType === m.FUNCTION) {
+      } else if (i.tokenType === g.FUNCTION) {
         const n = t.pop();
         if (n === void 0)
           throw this._isValid = !1, new Error(`The parameters for ${i.token} is not defined`);
@@ -851,20 +891,20 @@ class _i {
     return +e.toFixed(t);
   }
 }
-var T, v, Ye, xt, Le, qt, kt;
+var T, b, Qe, xt, ze, Mt, kt;
 const I = class I {
   constructor(e) {
-    p(this, Ye);
-    p(this, T);
-    p(this, v);
+    m(this, Qe);
+    m(this, T);
+    m(this, b);
     /**
      * Clone the current Monom.
      */
     a(this, "clone", () => {
       const e = new I();
       e.coefficient = s(this, T).clone();
-      for (const t in s(this, v))
-        e.setLetter(t, s(this, v)[t].clone());
+      for (const t in s(this, b))
+        e.setLetter(t, s(this, b)[t].clone());
       return e;
     });
     /**
@@ -874,24 +914,24 @@ const I = class I {
     a(this, "add", (...e) => {
       for (const t of e) {
         const i = t instanceof I ? t : new I(t);
-        this.isSameAs(i) ? (this.isZero() && q(this, Ye, xt).call(this, i), s(this, T).add(i.coefficient)) : console.log("Add monom: " + this.display + " is not similar with ", i.display);
+        this.isSameAs(i) ? (this.isZero() && A(this, Qe, xt).call(this, i), s(this, T).add(i.coefficient)) : console.log("Add monom: " + this.display + " is not similar with ", i.display);
       }
       return this;
     });
-    a(this, "containsRationalPower", () => Object.values(s(this, v)).some((e) => e.isRational()));
+    a(this, "containsRationalPower", () => Object.values(s(this, b)).some((e) => e.isRational()));
     /**
      * Get the degree of a monom. If no setLetter is given, the result will be the global degree.
      * @param letter (string) Letter to get to degree (power)
      */
-    a(this, "degree", (e) => this.variables.length === 0 ? new u().zero() : e === void 0 ? Object.values(s(this, v)).reduce((t, i) => t.clone().add(i)) : this.hasVariable(e) ? s(this, v)[e].clone() : new u().zero());
+    a(this, "degree", (e) => this.variables.length === 0 ? new u().zero() : e === void 0 ? Object.values(s(this, b)).reduce((t, i) => t.clone().add(i)) : this.hasVariable(e) ? s(this, b)[e].clone() : new u().zero());
     /**
      * Derivative the monom
      * @param letter
      */
     a(this, "derivative", (e) => {
       if (e === void 0 && (e = "x"), this.hasVariable(e)) {
-        const t = s(this, v)[e].clone(), i = this.clone();
-        return s(i, v)[e].subtract(1), s(i, T).multiply(new u(t.clone())), i;
+        const t = s(this, b)[e].clone(), i = this.clone();
+        return s(i, b)[e].subtract(1), s(i, T).multiply(new u(t.clone())), i;
       } else
         return new I().zero();
     });
@@ -904,7 +944,7 @@ const I = class I {
         const i = t instanceof I ? t : new I(t);
         s(this, T).divide(i.coefficient);
         for (const n in i.literal)
-          s(this, v)[n] = this.hasVariable(n) ? s(this, v)[n].subtract(i.literal[n]) : i.literal[n].clone().opposite(), s(this, v)[n].isZero() && this.removeVariable(n);
+          s(this, b)[n] = this.hasVariable(n) ? s(this, b)[n].subtract(i.literal[n]) : i.literal[n].clone().opposite(), s(this, b)[n].isZero() && this.removeVariable(n);
       }
       return this;
     });
@@ -916,16 +956,16 @@ const I = class I {
     a(this, "evaluate", (e, t) => {
       if (t === !0) {
         if (e instanceof u)
-          return s(this, Le).call(this, e.value);
+          return s(this, ze).call(this, e.value);
         if (e instanceof vt)
           return new u().invalid();
         if (typeof e == "number")
-          return s(this, Le).call(this, e);
+          return s(this, ze).call(this, e);
         if (typeof e == "object") {
           const n = {};
           for (const r in e)
             n[r] = new u(e[r]).value;
-          return s(this, Le).call(this, n);
+          return s(this, ze).call(this, n);
         }
       }
       const i = this.coefficient.clone();
@@ -938,9 +978,9 @@ const I = class I {
       if (typeof e == "object") {
         if (this.variables.length === 0)
           return this.coefficient;
-        for (const n in s(this, v)) {
+        for (const n in s(this, b)) {
           const r = new u(e[n]);
-          i.multiply(r.pow(s(this, v)[n]));
+          i.multiply(r.pow(s(this, b)[n]));
         }
       }
       return i;
@@ -950,11 +990,11 @@ const I = class I {
      * Determine if a monom contains a setLetter in it's literal part
      * @param letter
      */
-    a(this, "hasVariable", (e) => Object.hasOwn(s(this, v), e ?? "x"));
+    a(this, "hasVariable", (e) => Object.hasOwn(s(this, b), e ?? "x"));
     a(this, "inverse", () => {
       s(this, T).opposite();
-      for (const e in s(this, v))
-        s(this, v)[e].opposite();
+      for (const e in s(this, b))
+        s(this, b)[e].opposite();
       return this;
     });
     a(this, "isDivisible", (e) => {
@@ -992,7 +1032,7 @@ const I = class I {
         return !1;
       if (!this.isZero() && !e.isZero()) {
         for (const r of n)
-          if (!this.hasVariable(r) || !e.hasVariable(r) || !s(this, v)[r].isEqual(e.literal[r]))
+          if (!this.hasVariable(r) || !e.hasVariable(r) || !s(this, b)[r].isEqual(e.literal[r]))
             return !1;
       }
       return !0;
@@ -1011,14 +1051,14 @@ const I = class I {
         const i = t instanceof I ? t : new I(t);
         s(this, T).multiply(i.coefficient);
         for (const n in i.literal)
-          this.hasVariable(n) ? s(this, v)[n].add(i.literal[n]) : s(this, v)[n] = i.literal[n].clone();
+          this.hasVariable(n) ? s(this, b)[n].add(i.literal[n]) : s(this, b)[n] = i.literal[n].clone();
       }
       return this;
     });
     /**
      * Create a one value monom
      */
-    a(this, "one", () => (h(this, T, new u().one()), h(this, v, {}), this));
+    a(this, "one", () => (h(this, T, new u().one()), h(this, b, {}), this));
     /**
      * Get the opposite
      * Returns a monom.
@@ -1030,8 +1070,8 @@ const I = class I {
      */
     a(this, "pow", (e) => {
       s(this, T).pow(e);
-      for (const t in s(this, v))
-        s(this, v)[t].multiply(e);
+      for (const t in s(this, b))
+        s(this, b)[t].multiply(e);
       return this;
     });
     a(this, "primitive", (e) => {
@@ -1042,8 +1082,8 @@ const I = class I {
     });
     a(this, "reduce", () => {
       this.coefficient.reduce();
-      for (const e in s(this, v))
-        s(this, v)[e].isZero() && this.removeVariable(e);
+      for (const e in s(this, b))
+        s(this, b)[e].isZero() && this.removeVariable(e);
       return this;
     });
     /**
@@ -1057,15 +1097,15 @@ const I = class I {
      * @param letter (string) Letter to change
      * @param pow (number) Power of the setLetter (must be positive integer.
      */
-    a(this, "setLetter", (e, t) => t instanceof u ? (this.hasVariable(e) && t.isZero() && this.removeVariable(e), s(this, v)[e] = t.clone(), this) : this.setLetter(e, new u(t)));
+    a(this, "setLetter", (e, t) => t instanceof u ? (this.hasVariable(e) && t.isZero() && this.removeVariable(e), s(this, b)[e] = t.clone(), this) : this.setLetter(e, new u(t)));
     /**
      * Return the square root of a monom
      */
     a(this, "sqrt", () => {
       if (this.isSquare()) {
         s(this, T).sqrt();
-        for (const e in s(this, v))
-          s(this, v)[e].clone().divide(2);
+        for (const e in s(this, b))
+          s(this, b)[e].clone().divide(2);
       }
       return this;
     });
@@ -1076,57 +1116,57 @@ const I = class I {
     a(this, "subtract", (...e) => {
       for (const t of e) {
         const i = t instanceof I ? t : new I(t);
-        this.isSameAs(i) ? (this.isZero() && q(this, Ye, xt).call(this, i), s(this, T).add(i.clone().coefficient.opposite())) : console.log("Subtract: Is not similar: ", i.display);
+        this.isSameAs(i) ? (this.isZero() && A(this, Qe, xt).call(this, i), s(this, T).add(i.clone().coefficient.opposite())) : console.log("Subtract: Is not similar: ", i.display);
       }
       return this;
     });
     /**
      * Create a zero value monom
      */
-    a(this, "zero", () => (h(this, T, new u().zero()), h(this, v, {}), this));
-    p(this, Le, (e) => {
+    a(this, "zero", () => (h(this, T, new u().zero()), h(this, b, {}), this));
+    m(this, ze, (e) => {
       let t = this.coefficient.value;
       if (typeof e == "number") {
         const i = {}, n = this.variables[0];
-        return i[n] = e, s(this, Le).call(this, i);
+        return i[n] = e, s(this, ze).call(this, i);
       }
       if (e instanceof u) {
         const i = {};
-        return i[this.variables[0]] = new u(e).value, s(this, Le).call(this, i);
+        return i[this.variables[0]] = new u(e).value, s(this, ze).call(this, i);
       }
       if (e instanceof vt)
         return NaN;
       if (typeof e == "object") {
         if (this.variables.length === 0)
           return this.coefficient.value;
-        for (const i in s(this, v)) {
+        for (const i in s(this, b)) {
           const n = e[i];
-          n instanceof u ? t *= n.value ** s(this, v)[i].value : t *= n ** s(this, v)[i].value;
+          n instanceof u ? t *= n.value ** s(this, b)[i].value : t *= n ** s(this, b)[i].value;
         }
       }
       return t;
     });
-    p(this, qt, (e) => {
-      const i = new Qt().parse(e).rpn, n = [];
+    m(this, Mt, (e) => {
+      const i = new Kt().parse(e).rpn, n = [];
       if (i.length === 0)
         return this.zero(), this;
       if (i.length === 1) {
         const r = i[0];
-        return this.one(), r.tokenType === m.COEFFICIENT ? this.coefficient = new u(r.token) : r.tokenType === m.VARIABLE && this.setLetter(r.token, 1), this;
+        return this.one(), r.tokenType === g.COEFFICIENT ? this.coefficient = new u(r.token) : r.tokenType === g.VARIABLE && this.setLetter(r.token, 1), this;
       } else
         for (const r of i)
           s(this, kt).call(this, n, r);
       return this.one(), this.multiply(n[0]), this;
     });
-    p(this, kt, (e, t) => {
-      var d;
+    m(this, kt, (e, t) => {
+      var f;
       let i, n, r, l, c;
-      if (t.tokenType === m.COEFFICIENT)
+      if (t.tokenType === g.COEFFICIENT)
         e.push(new I(new u(t.token)));
-      else if (t.tokenType === m.VARIABLE) {
-        const f = new I().one();
-        f.setLetter(t.token, 1), e.push(f.clone());
-      } else if (t.tokenType === m.OPERATION)
+      else if (t.tokenType === g.VARIABLE) {
+        const d = new I().one();
+        d.setLetter(t.token, 1), e.push(d.clone());
+      } else if (t.tokenType === g.OPERATION)
         switch (t.token) {
           case "-":
             n = e.pop() ?? new I().zero(), i = e.pop() ?? new I().zero(), e.push(i.subtract(n));
@@ -1138,12 +1178,12 @@ const I = class I {
             n = e.pop() ?? new I().one(), i = e.pop() ?? new I().one(), e.push(i.divide(n));
             break;
           case "^": {
-            c = ((d = e.pop()) == null ? void 0 : d.coefficient) ?? new u().one(), r = e.pop() ?? new I().one(), l = r.variables[0], l && r.setLetter(l, c), e.push(r);
+            c = ((f = e.pop()) == null ? void 0 : f.coefficient) ?? new u().one(), r = e.pop() ?? new I().one(), l = r.variables[0], l && r.setLetter(l, c), e.push(r);
             break;
           }
         }
     });
-    return h(this, T, new u().zero()), h(this, v, {}), e !== void 0 && this.parse(e), this;
+    return h(this, T, new u().zero()), h(this, b, {}), e !== void 0 && this.parse(e), this;
   }
   // -----------------------------------------
   /**
@@ -1151,7 +1191,7 @@ const I = class I {
    * @param inputStr
    */
   parse(e) {
-    return h(this, T, new u()), h(this, v, {}), typeof e == "string" ? s(this, qt).call(this, e) : typeof e == "number" ? h(this, T, new u(e)) : e instanceof u ? h(this, T, e.clone()) : e instanceof I && (h(this, T, s(e, T).clone()), q(this, Ye, xt).call(this, e)), this;
+    return h(this, T, new u()), h(this, b, {}), typeof e == "string" ? s(this, Mt).call(this, e) : typeof e == "number" ? h(this, T, new u(e)) : e instanceof u ? h(this, T, e.clone()) : e instanceof I && (h(this, T, s(e, T).clone()), A(this, Qe, xt).call(this, e)), this;
   }
   /**
    * Get the coefficient \\(k\\) of the Monom \\(k\\cdot x^{n}\\)
@@ -1173,9 +1213,9 @@ const I = class I {
    */
   get display() {
     let e = "";
-    const t = Object.keys(s(this, v)).sort();
+    const t = Object.keys(s(this, b)).sort();
     for (const i of t)
-      s(this, v)[i].isNotZero() && (e += i, s(this, v)[i].isNotEqual(1) && (e += `^(${s(this, v)[i].display})`));
+      s(this, b)[i].isNotZero() && (e += i, s(this, b)[i].isNotEqual(1) && (e += `^(${s(this, b)[i].display})`));
     return e === "" ? s(this, T).value != 0 ? s(this, T).display : "" : s(this, T).value === 1 ? e : s(this, T).value === -1 ? `-${e}` : s(this, T).value === 0 ? "0" : `${s(this, T).display}${e}`;
   }
   get dividers() {
@@ -1185,7 +1225,7 @@ const I = class I {
       return [this.clone()];
     if (this.coefficient.numerator > 1e6)
       return [this.clone()];
-    const e = G.dividers(Math.abs(this.coefficient.numerator));
+    const e = U.dividers(Math.abs(this.coefficient.numerator));
     let t = [];
     for (const n in this.literal)
       t = this._getLiteralDividers(t, n);
@@ -1217,28 +1257,27 @@ const I = class I {
    * @returns {literalType}
    */
   get literal() {
-    return s(this, v);
+    return s(this, b);
   }
   /**
    * Set the literal part of the monom. Must be a dictionary {x: Fraction, y: Fraction, ...}
    * @param {literalType<Fraction>} L
    */
   set literal(e) {
-    h(this, v, e);
+    h(this, b, e);
   }
   /**
    * Get the literal square roots of the Monom.
-   * TODO: remove this getter ? Is it used and is it correct ?
    * @returns {literalType<Fraction>}
    */
   get literalSqrt() {
     if (this.isLiteralSquare()) {
       const e = {};
-      for (const t in s(this, v))
-        e[t] = s(this, v)[t].clone().sqrt();
+      for (const t in s(this, b))
+        e[t] = s(this, b)[t].clone().sqrt();
       return e;
     } else
-      return s(this, v);
+      return s(this, b);
   }
   /**
    * Set the literal part of the monom from a string
@@ -1246,28 +1285,28 @@ const I = class I {
    */
   set literalStr(e) {
     for (const t of [...e.matchAll(/([a-z])\^([+-]?[0-9]+)/g)])
-      t[1] in s(this, v) || (s(this, v)[t[1]] = new u().zero()), s(this, v)[t[1]].add(+t[2]);
+      t[1] in s(this, b) || (s(this, b)[t[1]] = new u().zero()), s(this, b)[t[1]].add(+t[2]);
     for (const t of [...e.matchAll(/([a-z](?!\^))/g)])
-      t[1] in s(this, v) || (s(this, v)[t[1]] = new u().zero()), s(this, v)[t[1]].add(1);
+      t[1] in s(this, b) || (s(this, b)[t[1]] = new u().zero()), s(this, b)[t[1]].add(1);
   }
   get plotFunction() {
     let e = "";
-    const t = Object.keys(s(this, v)).sort();
+    const t = Object.keys(s(this, b)).sort();
     for (const i of t)
-      s(this, v)[i].isNotZero() && (e += (e === "" ? "" : "*") + i, s(this, v)[i].isNotEqual(1) && (e += `^(${s(this, v)[i].display})`));
+      s(this, b)[i].isNotZero() && (e += (e === "" ? "" : "*") + i, s(this, b)[i].isNotEqual(1) && (e += `^(${s(this, b)[i].display})`));
     return e === "" ? s(this, T).value != 0 ? s(this, T).display : "" : s(this, T).value === 1 ? e : s(this, T).value === -1 ? `-${e}` : s(this, T).value === 0 ? "0" : `${s(this, T).display}*${e}`;
   }
   removeVariable(e) {
-    delete s(this, v)[e];
+    delete s(this, b)[e];
   }
   /**
    * Get the tex output of the monom
    */
   get tex() {
     let e = "";
-    const t = Object.keys(s(this, v)).sort();
+    const t = Object.keys(s(this, b)).sort();
     for (const i of t)
-      s(this, v)[i].isNotZero() && (e += i, s(this, v)[i].isNotEqual(1) && (e += `^{ ${s(this, v)[i].tfrac.tex} }`));
+      s(this, b)[i].isNotZero() && (e += i, s(this, b)[i].isNotEqual(1) && (e += `^{ ${s(this, b)[i].tfrac.tex} }`));
     return e === "" ? s(this, T).value != 0 ? s(this, T).frac.tex : "0" : s(this, T).value === 1 ? e : s(this, T).value === -1 ? `-${e}` : s(this, T).value === 0 ? "0" : `${s(this, T).frac.tex}${e}`;
   }
   // Getter helpers.
@@ -1298,14 +1337,14 @@ const I = class I {
     return i;
   }
 };
-T = new WeakMap(), v = new WeakMap(), Ye = new WeakSet(), xt = function(e) {
+T = new WeakMap(), b = new WeakMap(), Qe = new WeakSet(), xt = function(e) {
   for (const t in e.literal)
-    s(this, v)[t] = e.literal[t].clone();
-}, Le = new WeakMap(), qt = new WeakMap(), kt = new WeakMap(), a(I, "gcd", (...e) => {
+    s(this, b)[t] = e.literal[t].clone();
+}, ze = new WeakMap(), Mt = new WeakMap(), kt = new WeakMap(), a(I, "gcd", (...e) => {
   for (const r of e)
     if (r.containsRationalPower())
       return new I().zero();
-  const t = new I(), i = G.gcd(...e.map((r) => r.coefficient.numerator)), n = G.lcm(...e.map((r) => r.coefficient.denominator));
+  const t = new I(), i = U.gcd(...e.map((r) => r.coefficient.numerator)), n = U.lcm(...e.map((r) => r.coefficient.denominator));
   t.coefficient = new u(i, n).reduce();
   for (const r of e) {
     for (const l in t.literal)
@@ -1324,36 +1363,36 @@ a(I, "xMultiply", (...e) => {
     t.multiply(i);
   return t;
 });
-let M = I;
-function wi(o, e = !0) {
+let k = I;
+function yi(o, e = !0) {
   return e ? `\\left( ${o} \\right)` : `(${o})`;
 }
-function Re(o, e, t, i, n) {
-  return o.map((r, l) => i !== void 0 && l < i || n !== void 0 && l > n ? r : r === e ? t : r);
+function Et(o, e, t, i, n) {
+  return o.map((r, l) => r === e ? t : r);
 }
-var Ue, y, Qe, ut, Ke, ft, Mt, Ct, $t, Je, It, dt, Pt, Bt, St, Rt, Oi, zt, Lt;
-const C = class C {
+var je, y, Ke, ft, Je, dt, Ct, $t, Pt, _e, It, pt, Bt, St, Rt, zt, Mi, Lt, Vt;
+const $ = class $ {
   constructor(e, ...t) {
-    p(this, Rt);
-    p(this, Ue);
-    p(this, y);
-    p(this, Qe);
-    p(this, ut, !1);
+    m(this, zt);
+    m(this, je);
+    m(this, y);
+    m(this, Ke);
+    m(this, ft, !1);
     /**
      * Parse a string to a polynom.
      * @param inputStr
      * @param values
      */
     a(this, "parse", (e, ...t) => {
-      if (h(this, y, []), h(this, Ue, []), typeof e == "string")
-        return q(this, Rt, Oi).call(this, e, ...t);
-      if ((typeof e == "number" || e instanceof u || e instanceof M) && t.length === 0)
-        s(this, y).push(new M(e));
-      else if (e instanceof M && t.length > 0)
-        s(this, y).push(new M(e)), t.forEach((i) => {
-          s(this, y).push(new M(i));
+      if (h(this, y, []), h(this, je, []), typeof e == "string")
+        return A(this, zt, Mi).call(this, e, ...t);
+      if ((typeof e == "number" || e instanceof u || e instanceof k) && t.length === 0)
+        s(this, y).push(new k(e));
+      else if (e instanceof k && t.length > 0)
+        s(this, y).push(new k(e)), t.forEach((i) => {
+          s(this, y).push(new k(i));
         });
-      else if (e instanceof C)
+      else if (e instanceof $)
         for (const i of e.monoms)
           s(this, y).push(i.clone());
       return this;
@@ -1362,18 +1401,18 @@ const C = class C {
      * Clone the polynom
      */
     a(this, "clone", () => {
-      const e = new C(), t = [];
+      const e = new $(), t = [];
       for (const i of s(this, y))
         t.push(i.clone());
       return e.monoms = t, e;
     });
     a(this, "add", (...e) => {
       for (const t of e)
-        t instanceof C ? h(this, y, s(this, y).concat(t.monoms)) : t instanceof M ? s(this, y).push(t.clone()) : typeof t == "number" && Number.isSafeInteger(t) ? s(this, y).push(new M(t.toString())) : s(this, y).push(new M(t));
+        t instanceof $ ? h(this, y, s(this, y).concat(t.monoms)) : t instanceof k ? s(this, y).push(t.clone()) : typeof t == "number" && Number.isSafeInteger(t) ? s(this, y).push(new k(t.toString())) : s(this, y).push(new k(t));
       return this.reduce();
     });
     a(this, "commonMonom", () => {
-      const e = new M().one(), t = this.gcdNumerator(), i = this.gcdDenominator(), n = this.degree();
+      const e = new k().one(), t = this.gcdNumerator(), i = this.gcdDenominator(), n = this.degree();
       e.coefficient = new u(t, i);
       for (const r of this.variables) {
         e.setLetter(r, n);
@@ -1390,28 +1429,28 @@ const C = class C {
       return t;
     });
     a(this, "derivative", (e) => {
-      const t = new C();
+      const t = new $();
       for (const i of s(this, y))
         t.add(i.derivative(e));
-      return t;
+      return t.reduce();
     });
     a(this, "divide", (e) => {
       if (e instanceof u)
-        return s(this, ft).call(this, e);
+        return s(this, dt).call(this, e);
       if (typeof e == "number" && Number.isSafeInteger(e))
-        return s(this, Mt).call(this, e);
-      if (e instanceof M)
-        return this.divide(new C(e));
-      if (e instanceof C) {
+        return s(this, Ct).call(this, e);
+      if (e instanceof k)
+        return this.divide(new $(e));
+      if (e instanceof $) {
         if (e.monoms.length === 1 && e.variables.length === 0)
-          return s(this, ft).call(this, e.monoms[0].coefficient);
+          return s(this, dt).call(this, e.monoms[0].coefficient);
         {
           const { quotient: t, reminder: i } = this.euclidean(e);
           if (i.isZero())
             return h(this, y, t.monoms), this;
         }
       } else if (typeof e == "string")
-        return this.divide(new C(e));
+        return this.divide(new $(e));
       throw new Error(`Cannot divide by ${e}`);
     });
     a(this, "empty", () => (h(this, y, []), this));
@@ -1421,21 +1460,21 @@ const C = class C {
      * returns {quotient: Polynom, reminder: Polynom}
      */
     a(this, "euclidean", (e) => {
-      const t = e.variables[0], i = new C().zero(), n = this.clone().reorder(t);
+      const t = e.variables[0], i = new $().zero(), n = this.clone().reorder(t);
       if (e.variables.length === 0)
         return {
           quotient: this.clone().divide(e).reduce(),
-          reminder: new C().zero()
+          reminder: new $().zero()
         };
       const r = e.monomByDegree(void 0, t), l = e.degree(t);
-      let c, d = this.degree(t).value * 2;
-      for (; n.degree(t).isGeq(l) && d > 0 && (d--, c = n.monomByDegree(void 0, t).clone().divide(r), !(!c.isZero() && (i.add(c), n.subtract(e.clone().multiply(c)).reduce(), c.degree(t).isZero()))); )
+      let c, f = this.degree(t).value * 2;
+      for (; n.degree(t).isGeq(l) && f > 0 && (f--, c = n.monomByDegree(void 0, t).clone().divide(r), !(!c.isZero() && (i.add(c), n.subtract(e.clone().multiply(c)).reduce(), c.degree(t).isZero()))); )
         ;
       return i.reduce(), n.reduce(), { quotient: i, reminder: n };
     });
     a(this, "evaluate", (e, t) => {
       if (t)
-        return s(this, Ct).call(this, e);
+        return s(this, $t).call(this, e);
       const i = new u().zero();
       return s(this, y).forEach((n) => {
         i.add(n.evaluate(e, t));
@@ -1450,7 +1489,7 @@ const C = class C {
       let t = [], i = this.clone().reorder();
       const n = i.commonMonom();
       if (i.monomByDegree().coefficient.isStrictlyNegative() && n.coefficient.isStrictlyPositive() && !n.isOne() && n.opposite(), !n.isOne()) {
-        const c = new C(n);
+        const c = new $(n);
         t = [c.clone()], i = i.euclidean(c).quotient;
       }
       let r = i.degree().clone().multiply(2).value, l = 1;
@@ -1464,22 +1503,22 @@ const C = class C {
         } else {
           let c = s(this, It).call(this, i, l, e ?? "x");
           for (l = i.degree(e).value; c.length > 0; ) {
-            const d = c[0];
-            if (!i.isDividableBy(d))
+            const f = c[0];
+            if (!i.isDividableBy(f))
               c.shift();
             else {
-              const f = i.euclidean(d);
-              t.push(d), i = f.quotient.clone(), c = c.filter((g) => {
-                const N = i.monoms[0], V = i.monoms[i.monoms.length - 1], pe = g.monoms[0], he = g.monoms[g.monoms.length - 1];
-                return V.isDivisible(he) ? N.isDivisible(pe) : !1;
+              const d = i.euclidean(f);
+              t.push(f), i = d.quotient.clone(), c = c.filter((p) => {
+                const w = i.monoms[0], V = i.monoms[i.monoms.length - 1], ce = p.monoms[0], ue = p.monoms[p.monoms.length - 1];
+                return V.isDivisible(ue) ? w.isDivisible(ce) : !1;
               });
             }
           }
         }
-      return i.isOne() || t.push(i.clone()), h(this, Ue, t), s(this, Ue);
+      return i.isOne() || t.push(i.clone()), h(this, je, t), s(this, je);
     });
-    a(this, "gcdDenominator", () => G.gcd(...this.getDenominators()));
-    a(this, "gcdNumerator", () => G.gcd(...this.getNumerators()));
+    a(this, "gcdDenominator", () => U.gcd(...this.getDenominators()));
+    a(this, "gcdNumerator", () => U.gcd(...this.getNumerators()));
     // Next functions are used for for commonMonom, which is used in the factorize method.
     a(this, "getDenominators", () => {
       const e = [];
@@ -1493,7 +1532,7 @@ const C = class C {
         e.push(t.coefficient.numerator);
       return e;
     });
-    a(this, "getZeroes", () => this.degree().isZero() ? [] : (this.roots = new Nt(this.clone()).solve(), this.roots));
+    a(this, "getZeroes", () => this.degree().isZero() ? [] : (this.roots = new At(this.clone()).solve(), this.roots));
     a(this, "integrate", (e, t, i = "x") => {
       const n = this.primitive(i), r = {}, l = {};
       return r[i] = new u(e), l[i] = new u(t), n.evaluate(l).subtract(n.evaluate(r));
@@ -1504,7 +1543,7 @@ const C = class C {
       if (i.includes("(") || i.includes(")"))
         return !1;
       try {
-        t = new C(e);
+        t = new $(e);
       } catch {
         return !1;
       }
@@ -1519,12 +1558,12 @@ const C = class C {
         return t.isZero();
       }
     });
-    a(this, "isEqual", (e) => s(this, Ke).call(this, e, "="));
-    a(this, "isOppositeAt", (e) => s(this, Ke).call(this, e.clone().opposite(), "="));
+    a(this, "isEqual", (e) => s(this, Je).call(this, e, "="));
+    a(this, "isOppositeAt", (e) => s(this, Je).call(this, e.clone().opposite(), "="));
     a(this, "isReduced", (e) => {
       if (!this.isDeveloped(e))
         return !1;
-      const t = new C(e);
+      const t = new $(e);
       if (t.monoms.length > this.monoms.length)
         return !1;
       for (const i of t.monoms)
@@ -1532,9 +1571,9 @@ const C = class C {
           return !1;
       return !1;
     });
-    a(this, "isSameAs", (e) => s(this, Ke).call(this, e, "same"));
-    a(this, "lcmDenominator", () => G.lcm(...this.getDenominators()));
-    a(this, "lcmNumerator", () => G.lcm(...this.getNumerators()));
+    a(this, "isSameAs", (e) => s(this, Je).call(this, e, "same"));
+    a(this, "lcmDenominator", () => U.lcm(...this.getDenominators()));
+    a(this, "lcmNumerator", () => U.lcm(...this.getNumerators()));
     a(this, "letters", () => {
       let e = /* @__PURE__ */ new Set();
       for (const t of s(this, y))
@@ -1556,7 +1595,7 @@ const C = class C {
       for (const n of s(i, y))
         if (n.degree(t).isEqual(e))
           return n.clone();
-      return new M().zero();
+      return new k().zero();
     });
     // Used in LinearSystem.tex
     a(this, "monomByLetter", (e) => {
@@ -1564,7 +1603,7 @@ const C = class C {
       for (const i of s(t, y))
         if (i.hasVariable(e))
           return i.clone();
-      return new M().zero();
+      return new k().zero();
     });
     a(this, "monomsByDegree", (e, t) => {
       if (e === void 0)
@@ -1574,8 +1613,8 @@ const C = class C {
         r.degree(t).isEqual(e) && i.push(r.clone());
       return i;
     });
-    a(this, "multiply", (e) => e instanceof C ? s(this, St).call(this, e) : e instanceof u ? s(this, dt).call(this, e) : e instanceof M ? s(this, Bt).call(this, e) : Number.isSafeInteger(e) && typeof e == "number" ? s(this, Pt).call(this, e) : this);
-    a(this, "one", () => (h(this, y, []), s(this, y).push(new M().one()), this));
+    a(this, "multiply", (e) => e instanceof $ ? s(this, Rt).call(this, e) : e instanceof u ? s(this, pt).call(this, e) : e instanceof k ? s(this, St).call(this, e) : Number.isSafeInteger(e) && typeof e == "number" ? s(this, Bt).call(this, e) : this);
+    a(this, "one", () => (h(this, y, []), s(this, y).push(new k().one()), this));
     // ------------------------------------------
     a(this, "opposite", () => (h(this, y, s(this, y).map((e) => e.opposite())), this));
     a(this, "pow", (e) => {
@@ -1584,14 +1623,14 @@ const C = class C {
       if (e < 0)
         return this.zero();
       if (e === 0)
-        return new C();
+        return new $();
       const t = this.clone();
       for (let i = 1; i < e; i++)
         this.multiply(t);
       return this.reduce();
     });
     a(this, "primitive", (e) => {
-      const t = new C();
+      const t = new $();
       for (const i of s(this, y))
         t.add(i.primitive(e));
       return t;
@@ -1600,14 +1639,15 @@ const C = class C {
       let e = 0;
       for (; e < s(this, y).length; ) {
         for (let t = e + 1; t < s(this, y).length; t++)
-          s(this, y)[e].isSameAs(s(this, y)[t]) && (s(this, y)[e].add(s(this, y)[t]), s(this, y).splice(t, 1), s(this, y)[e].isZero() && (s(this, y)[e] = new M().zero()), t--);
+          s(this, y)[e].isSameAs(s(this, y)[t]) && (s(this, y)[e].add(s(this, y)[t]), s(this, y).splice(t, 1), s(this, y)[e].isZero() && (s(this, y)[e] = new k().zero()), t--);
         e++;
       }
       h(this, y, s(this, y).filter((t) => !t.coefficient.isZero()));
       for (const t of s(this, y))
         t.coefficient.reduce();
-      return this.length === 0 ? new C().zero() : this.reorder();
+      return this.length === 0 ? new $().zero() : this.reorder();
     });
+    // ------------------------------------------
     a(this, "reorder", (e = "x", t) => {
       t === void 0 && (t = !1);
       const i = this.variables.filter((n) => n !== e);
@@ -1616,15 +1656,14 @@ const C = class C {
         if (l !== c)
           return t ? l - c : c - l;
         if (i.length > 0)
-          for (const d of i) {
-            const f = n.degree(d).value, g = r.degree(d).value;
-            if (f !== g)
-              return t ? f - g : g - f;
+          for (const f of i) {
+            const d = n.degree(f).value, p = r.degree(f).value;
+            if (d !== p)
+              return t ? d - p : p - d;
           }
         return 0;
       }), this;
     });
-    // ------------------------------------------
     /**
      * Replace a variable (letter) by a polynom.
      * @param letter
@@ -1632,22 +1671,22 @@ const C = class C {
      */
     a(this, "replaceBy", (e, t) => {
       let i;
-      const n = new C().zero();
+      const n = new $().zero();
       for (const r of this.monoms)
         !r.hasVariable(e) || r.literal[e].isZero() ? n.add(r.clone()) : (i = r.literal[e].clone(), r.removeVariable(e), n.add(t.clone().pow(Math.abs(i.numerator)).multiply(r)));
       return h(this, y, n.reduce().monoms), this;
     });
     a(this, "subtract", (...e) => {
       for (const t of e)
-        t instanceof C ? this.add(t.clone().opposite()) : t instanceof M ? s(this, y).push(t.clone().opposite()) : s(this, y).push(new M(t).opposite());
+        t instanceof $ ? this.add(t.clone().opposite()) : t instanceof k ? s(this, y).push(t.clone().opposite()) : s(this, y).push(new k(t).opposite());
       return this.reduce();
     });
     /**
      * Set the polynom to zero.
      * @returns {this}
      */
-    a(this, "zero", () => (h(this, y, []), s(this, y).push(new M().zero()), this));
-    p(this, Ke, (e, t) => {
+    a(this, "zero", () => (h(this, y, []), s(this, y).push(new k().zero()), this));
+    m(this, Je, (e, t) => {
       t === void 0 && (t = "=");
       const i = this.clone().reduce().reorder(), n = e.clone().reduce().reorder();
       switch (t) {
@@ -1659,47 +1698,47 @@ const C = class C {
           return !1;
       }
     });
-    p(this, ft, (e) => {
+    m(this, dt, (e) => {
       for (const t of s(this, y))
         t.coefficient.divide(e);
       return this;
     });
-    p(this, Mt, (e) => {
+    m(this, Ct, (e) => {
       const t = new u(e);
       for (const i of s(this, y))
         i.coefficient.divide(t);
       return this;
     });
-    p(this, Ct, (e) => {
+    m(this, $t, (e) => {
       let t = 0;
       return s(this, y).forEach((i) => {
         t += i.evaluate(e, !0);
       }), t;
     });
-    p(this, $t, (e) => {
-      var N;
-      let t, i, n, r, l, c, d, f, g;
+    m(this, Pt, (e) => {
+      var w;
+      let t, i, n, r, l, c, f, d, p;
       if (this.numberOfVars === 1)
-        return n = this.monomByDegree(2, e).coefficient, r = this.monomByDegree(1, e).coefficient, l = this.monomByDegree(0, e).coefficient, c = r.clone().pow(2).subtract(n.clone().multiply(l).multiply(4)), c.isZero() ? (d = r.clone().opposite().divide(n.clone().multiply(2)), t = new C(e).subtract(d.display).multiply(d.denominator), i = new C(e).subtract(d.display).multiply(d.denominator), g = n.divide(d.denominator).divide(d.denominator), g.isOne() ? [t, i] : [new C(g.display), t, i]) : c.isPositive() && c.isSquare() ? (d = r.clone().opposite().add(c.clone().sqrt()).divide(n.clone().multiply(2)), f = r.clone().opposite().subtract(c.clone().sqrt()).divide(n.clone().multiply(2)), g = n.divide(d.denominator).divide(f.denominator), g.isOne() ? [
-          new C(e).subtract(d.display).multiply(d.denominator),
-          new C(e).subtract(f.display).multiply(f.denominator)
+        return n = this.monomByDegree(2, e).coefficient, r = this.monomByDegree(1, e).coefficient, l = this.monomByDegree(0, e).coefficient, c = r.clone().pow(2).subtract(n.clone().multiply(l).multiply(4)), c.isZero() ? (f = r.clone().opposite().divide(n.clone().multiply(2)), t = new $(e).subtract(f.display).multiply(f.denominator), i = new $(e).subtract(f.display).multiply(f.denominator), p = n.divide(f.denominator).divide(f.denominator), p.isOne() ? [t, i] : [new $(p.display), t, i]) : c.isPositive() && c.isSquare() ? (f = r.clone().opposite().add(c.clone().sqrt()).divide(n.clone().multiply(2)), d = r.clone().opposite().subtract(c.clone().sqrt()).divide(n.clone().multiply(2)), p = n.divide(f.denominator).divide(d.denominator), p.isOne() ? [
+          new $(e).subtract(f.display).multiply(f.denominator),
+          new $(e).subtract(d.display).multiply(d.denominator)
         ] : [
-          new C(g.display),
-          new C(e).subtract(d.display).multiply(d.denominator),
-          new C(e).subtract(f.display).multiply(f.denominator)
+          new $(p.display),
+          new $(e).subtract(f.display).multiply(f.denominator),
+          new $(e).subtract(d.display).multiply(d.denominator)
         ]) : [this.clone()];
       if (n = this.monomByDegree(2, e), r = this.monomByDegree(1, e), l = this.monomByDegree(0, e), n.isLiteralSquare() && l.isLiteralSquare() && r.clone().pow(2).isSameAs(n.clone().multiply(l))) {
-        const V = new C("x", n.coefficient, r.coefficient, l.coefficient), pe = s(N = V, $t).call(N, "x"), he = [];
-        let ae;
-        if (pe.length >= 2) {
-          for (const ge of pe)
-            ge.degree().isZero() ? he.push(ge.clone()) : (ae = ge.clone(), ae.monoms[0].literal = n.literalSqrt, ae.monoms[1].literal = l.literalSqrt, he.push(ae.clone()));
-          return he;
+        const V = new $("x", n.coefficient, r.coefficient, l.coefficient), ce = s(w = V, Pt).call(w, "x"), ue = [];
+        let fe;
+        if (ce.length >= 2) {
+          for (const ye of ce)
+            ye.degree().isZero() ? ue.push(ye.clone()) : (fe = ye.clone(), fe.monoms[0].literal = n.literalSqrt, fe.monoms[1].literal = l.literalSqrt, ue.push(fe.clone()));
+          return ue;
         }
       }
       return [this.clone()];
     });
-    p(this, Je, (e, t, i, n) => {
+    m(this, _e, (e, t, i, n) => {
       let r = "";
       for (const l of s(this, y)) {
         if (l.coefficient.value === 0)
@@ -1709,56 +1748,56 @@ const C = class C {
       }
       return i === !0 && this.length > 1 && (e === "tex" ? r = `\\left( ${r} \\right)` : r = `(${r})`), r === "" && (r = "0"), r;
     });
-    p(this, It, (e, t, i) => {
+    m(this, It, (e, t, i) => {
       const n = e.monoms[0].dividers, r = e.monoms[e.monoms.length - 1].dividers, l = [];
       return n.forEach((c) => {
-        c.degree(i).isLeq(t) && r.forEach((d) => {
-          c.degree(i).isNotEqual(d.degree(i)) && (l.push(new C(c, d)), l.push(new C(c, d.clone().opposite())));
+        c.degree(i).isLeq(t) && r.forEach((f) => {
+          c.degree(i).isNotEqual(f.degree(i)) && (l.push(new $(c, f)), l.push(new $(c, f.clone().opposite())));
         });
       }), l;
     });
-    p(this, dt, (e) => {
+    m(this, pt, (e) => {
       for (const t of s(this, y))
         t.coefficient.multiply(e);
       return this.reduce();
     });
-    p(this, Pt, (e) => s(this, dt).call(this, new u(e)));
-    p(this, Bt, (e) => {
+    m(this, Bt, (e) => s(this, pt).call(this, new u(e)));
+    m(this, St, (e) => {
       for (const t of s(this, y))
         t.multiply(e);
       return this.reduce();
     });
-    p(this, St, (e) => {
+    m(this, Rt, (e) => {
       const t = [];
       for (const i of s(this, y))
         for (const n of e.monoms)
-          t.push(M.xMultiply(i, n));
+          t.push(k.xMultiply(i, n));
       return h(this, y, t), this.reduce();
     });
     /**
      * Main parse using a shutting yard class
      * @param inputStr
      */
-    p(this, zt, (e) => {
-      const i = new Qt().parse(e).rpn;
+    m(this, Lt, (e) => {
+      const i = new Kt().parse(e).rpn;
       this.zero();
       const n = [];
       for (const r of i)
-        s(this, Lt).call(this, n, r);
+        s(this, Vt).call(this, n, r);
       return n.length === 1 && this.add(n[0]), this.reorder();
     });
-    p(this, Lt, (e, t) => {
+    m(this, Vt, (e, t) => {
       switch (t.tokenType) {
-        case m.COEFFICIENT:
-          e.push(new C(t.token));
+        case g.COEFFICIENT:
+          e.push(new $(t.token));
           break;
-        case m.VARIABLE:
-          e.push(new C().add(new M(t.token)));
+        case g.VARIABLE:
+          e.push(new $().add(new k(t.token)));
           break;
-        case m.CONSTANT:
+        case g.CONSTANT:
           console.log("Actually, not supported - will be added later !");
           break;
-        case m.OPERATION:
+        case g.OPERATION:
           if (e.length >= 2) {
             const i = e.pop(), n = e.pop();
             if (n === void 0 || i === void 0)
@@ -1789,18 +1828,36 @@ const C = class C {
           } else
             throw new Error("Error parsing the polynom");
           break;
-        case m.MONOM:
+        case g.MONOM:
           console.error("The monom token should not appear here");
           break;
-        case m.FUNCTION:
+        case g.FUNCTION:
           console.error("The function token should not appear here - might be introduced later.");
           break;
       }
     });
-    return h(this, y, []), h(this, Ue, []), h(this, Qe, []), e !== void 0 && this.parse(e, ...t), this;
+    return h(this, y, []), h(this, je, []), h(this, Ke, []), e !== void 0 && this.parse(e, ...t), this;
+  }
+  get tex() {
+    return s(this, _e).call(this, "tex");
   }
   get display() {
-    return s(this, Je).call(this);
+    return s(this, _e).call(this);
+  }
+  fromCoefficients(...e) {
+    h(this, y, []);
+    const t = "x";
+    return e.reverse().forEach((i, n) => {
+      const r = new k();
+      r.coefficient = new u(i), r.setLetter(t, n), s(this, y).push(r);
+    }), this;
+  }
+  getCoefficients() {
+    const e = this.clone().reorder(), t = this.degree().value + 1, i = new Array(t).fill(new u(0));
+    return e.monoms.forEach((n) => {
+      const r = t - n.degree().value - 1;
+      i[r] = n.coefficient.clone();
+    }), i;
   }
   hasVariable(e) {
     return this.variables.includes(e);
@@ -1830,49 +1887,38 @@ const C = class C {
     return this.variables.length;
   }
   get plotFunction() {
-    return s(this, Je).call(this, "tex", !1, !1, !0);
+    return s(this, _e).call(this, "tex", !1, !1, !0);
   }
   root() {
     throw new Error("Cannot take the root from a polynom");
   }
   get roots() {
-    return s(this, ut) ? s(this, Qe) : this.getZeroes();
+    return s(this, ft) ? s(this, Ke) : this.getZeroes();
   }
   set roots(e) {
-    h(this, ut, !0), h(this, Qe, e);
+    h(this, ft, !0), h(this, Ke, e);
   }
   sqrt() {
     throw new Error("Cannot take the square root from a polynom");
   }
-  tableOfSigns(e) {
-    const t = e ?? this.getZeroes(), i = this.roots;
-    if (e && i.some((r) => t.findIndex((l) => l.value === r.value) === -1))
-      throw new Error("Some roots cannot be found !");
-    let n = [""];
-    if (t.forEach(() => n.push("t", "")), t.length === 0)
-      n = Re(n, "", this.monomsByDegree()[0].coefficient.isPositive() ? "+" : "-");
-    else if (this.degree().isOne()) {
-      const r = this.monomsByDegree(1)[0].coefficient.sign(), l = t.findIndex((c) => c.value === i[0].value) * 2 + 1;
-      n[l] = "z", n = Re(n, "", r === 1 ? "-" : "+", 0, l), n = Re(n, "", r === 1 ? "+" : "-", l);
-    } else {
-      let r = 0, l, c = "+";
-      i.forEach((d, f) => {
-        const g = t.findIndex((V) => V.value === i[f].value);
-        f === 0 && (l = g * 2 + 1, n[l] = "z", n = Re(
-          n,
-          "",
-          this.evaluate(d.value - 1, !0) < 0 ? "-" : "+",
-          0,
-          l
-        )), r = g * 2 + 1, l = f === i.length - 1 ? n.length : t.findIndex((V) => V.value === i[f + 1].value) * 2 + 1;
-        const N = f === i.length - 1 ? d.value + 1 : (d.value + i[f + 1].value) / 2;
-        c = this.evaluate(N, !0) < 0 ? "-" : "+", n[r] = "z", n = Re(n, "", c, r, l), r = +l, l = n.length;
+  tableOfSigns() {
+    const e = this.roots;
+    let t = new Array(2 * e.length + 1).fill("").map((i, n) => n % 2 === 0 ? "" : "z");
+    if (t.length === 1) {
+      const [i] = this.getCoefficients().map((n) => n.value);
+      t = Et(t, "", i > 0 ? "+" : "-");
+    } else if (this.degree().isOne()) {
+      const [i] = this.getCoefficients().map((n) => n.value);
+      t[0] = i > 0 ? "-" : "+", t[1] = "z", t[2] = i > 0 ? "+" : "-";
+    } else
+      [
+        e[0].value - 1,
+        ...e.map((n, r) => r === e.length - 1 ? e[r].value + 1 : (e[r].value + e[r + 1].value) / 2)
+      ].forEach((n, r) => {
+        const l = this.evaluate({ x: n }, !0);
+        t[r * 2] = l > 0 ? "+" : "-";
       });
-    }
-    return { roots: t, signs: n };
-  }
-  get tex() {
-    return s(this, Je).call(this, "tex");
+    return { roots: e, signs: t };
   }
   get variables() {
     let e = [];
@@ -1884,60 +1930,62 @@ const C = class C {
     return this.getZeroes();
   }
 };
-Ue = new WeakMap(), y = new WeakMap(), Qe = new WeakMap(), ut = new WeakMap(), Ke = new WeakMap(), ft = new WeakMap(), Mt = new WeakMap(), Ct = new WeakMap(), $t = new WeakMap(), Je = new WeakMap(), It = new WeakMap(), dt = new WeakMap(), Pt = new WeakMap(), Bt = new WeakMap(), St = new WeakMap(), Rt = new WeakSet(), Oi = function(e, ...t) {
+je = new WeakMap(), y = new WeakMap(), Ke = new WeakMap(), ft = new WeakMap(), Je = new WeakMap(), dt = new WeakMap(), Ct = new WeakMap(), $t = new WeakMap(), Pt = new WeakMap(), _e = new WeakMap(), It = new WeakMap(), pt = new WeakMap(), Bt = new WeakMap(), St = new WeakMap(), Rt = new WeakMap(), zt = new WeakSet(), Mi = function(e, ...t) {
   if (t.length === 0) {
     if (e = "" + e, e !== "" && !isNaN(Number(e))) {
       this.empty();
-      const i = new M(e);
+      const i = new k(e);
       return this.add(i), this;
     }
-    return s(this, zt).call(this, e);
+    return s(this, Lt).call(this, e);
   } else if (/^[a-z]+/.test(e)) {
     this.empty();
     const i = t.map((n) => new u(n));
     if (e.length > 1) {
       const n = e.split("");
+      if (n.length < t.length - 2)
+        throw new Error("Too many factors for too few variables !");
       let r = 0;
       for (const l of i) {
-        const c = new M();
+        const c = new k();
         c.coefficient = l.clone(), c.literalStr = n[r] || "", this.add(c), r++;
       }
     } else {
       let n = i.length - 1;
       for (const r of i) {
-        const l = new M();
+        const l = new k();
         l.coefficient = r.clone(), l.literalStr = `${e}^${n}`, this.add(l), n--;
       }
     }
     return this;
   } else
     return this.zero();
-}, zt = new WeakMap(), Lt = new WeakMap();
-let O = C;
-var A, k, se, Dt, _e, Vt;
-const Ee = class Ee {
+}, Lt = new WeakMap(), Vt = new WeakMap();
+let O = $;
+var q, M, re, Dt, et, Zt;
+const Ne = class Ne {
   constructor(e, t, i) {
     // Left part of the equation
-    p(this, A);
+    m(this, q);
     // Right part of the equation
-    p(this, k);
+    m(this, M);
     // Signe of the equation
-    p(this, se);
+    m(this, re);
     // ------------------------------------------
     a(this, "parse", (e) => {
       const t = s(this, Dt).call(this, e);
       if (t === !1)
         throw new Error("The equation is not valid (no sign found)");
       const i = e.split(t);
-      return this.create(new O(i[0]), new O(i[1]), s(this, _e).call(this, t));
+      return this.create(new O(i[0]), new O(i[1]), s(this, et).call(this, t));
     });
-    a(this, "create", (e, t, i) => (h(this, A, e), h(this, k, t), h(this, se, s(this, _e).call(this, i ?? "=")), this));
-    a(this, "clone", () => new Ee(s(this, A).clone(), s(this, k).clone(), s(this, se)));
+    a(this, "create", (e, t, i) => (h(this, q, e), h(this, M, t), h(this, re, s(this, et).call(this, i ?? "=")), this));
+    a(this, "clone", () => new Ne(s(this, q).clone(), s(this, M).clone(), s(this, re)));
     /**
      * Get the degree of the equation
      * @param letter
      */
-    a(this, "degree", (e) => u.max(s(this, A).degree(e), s(this, k).degree(e)));
+    a(this, "degree", (e) => u.max(s(this, q).degree(e), s(this, M).degree(e)));
     /**
      * divide an equation by a given value (transformed as a fraction)
      *
@@ -1968,7 +2016,7 @@ const Ee = class Ee {
     /**
      * Determine if the equation contains more than one letter/variable.
      */
-    a(this, "isMultiVariable", () => s(this, A).isMultiVariable || s(this, k).isMultiVariable);
+    a(this, "isMultiVariable", () => s(this, q).isMultiVariable || s(this, M).isMultiVariable);
     // -----------------------------------------------
     // Equations helpers
     a(this, "isEqualTo", (e) => {
@@ -1983,47 +2031,47 @@ const Ee = class Ee {
       if (!this.degree(e).isOne() || this.isMultiVariable())
         return !1;
       let t;
-      s(this, A).subtract(s(this, k)), s(this, k).zero();
-      const i = [...s(this, A).monoms];
+      s(this, q).subtract(s(this, M)), s(this, M).zero();
+      const i = [...s(this, q).monoms];
       for (const r of i)
-        r.hasVariable(e) || (t = r.clone(), s(this, A).subtract(t), s(this, k).subtract(t));
-      if (s(this, A).length !== 1)
+        r.hasVariable(e) || (t = r.clone(), s(this, q).subtract(t), s(this, M).subtract(t));
+      if (s(this, q).length !== 1)
         return !1;
-      const n = s(this, A).monoms[0].coefficient.clone();
-      return s(this, A).divide(n), s(this, k).divide(n), this;
+      const n = s(this, q).monoms[0].coefficient.clone();
+      return s(this, q).divide(n), s(this, M).divide(n), this;
     });
     // -----------------------------------------------
     // Equations operations
     // -----------------------------------------------
-    a(this, "letters", () => [.../* @__PURE__ */ new Set([...s(this, A).letters(), ...s(this, k).letters()])]);
+    a(this, "letters", () => [.../* @__PURE__ */ new Set([...s(this, q).letters(), ...s(this, M).letters()])]);
     // -----------------------------------------------
     /**
      * Reorder will move all monoms containing a letter on the left, all the other on the right.
      */
-    a(this, "moveLeft", () => (s(this, A).subtract(s(this, k)), s(this, k).zero(), this));
+    a(this, "moveLeft", () => (s(this, q).subtract(s(this, M)), s(this, M).zero(), this));
     /**
      * Multiple an equation by a fraction value.
      * @param value
      */
     a(this, "multiply", (e) => {
       const t = new u(e);
-      return s(this, A).multiply(t), s(this, k).multiply(t), s(this, se) !== "=" && t.sign() === -1 && s(this, Vt).call(this), this;
+      return s(this, q).multiply(t), s(this, M).multiply(t), s(this, re) !== "=" && t.sign() === -1 && s(this, Zt).call(this), this;
     });
-    a(this, "opposite", () => (h(this, A, s(this, A).opposite()), h(this, k, s(this, k).opposite()), this));
-    a(this, "reorder", (e) => (s(this, A).subtract(s(this, k)), s(this, k).zero(), s(this, A).reorder(), e ? this : (s(this, A).monoms.filter((t) => t.degree().isZero()).forEach((t) => {
+    a(this, "opposite", () => (h(this, q, s(this, q).opposite()), h(this, M, s(this, M).opposite()), this));
+    a(this, "reorder", (e) => (s(this, q).subtract(s(this, M)), s(this, M).zero(), s(this, q).reorder(), e ? this : (s(this, q).monoms.filter((t) => t.degree().isZero()).forEach((t) => {
       const i = t.clone();
-      s(this, A).subtract(i), s(this, k).subtract(i);
-    }), s(this, A).reorder(), s(this, k).reorder(), this)));
+      s(this, q).subtract(i), s(this, M).subtract(i);
+    }), s(this, q).reorder(), s(this, M).reorder(), this)));
     // ------------------------------------------
-    a(this, "replaceBy", (e, t) => (s(this, A).replaceBy(e, t), s(this, k).replaceBy(e, t), this));
+    a(this, "replaceBy", (e, t) => (s(this, q).replaceBy(e, t), s(this, M).replaceBy(e, t), this));
     /**
      * Multiply by the lcm denominator and divide by the gcm numerators.
      */
-    a(this, "simplify", () => (this.multiply(G.lcm(...s(this, A).getDenominators(), ...s(this, k).getDenominators())), this.divide(G.gcd(...s(this, A).getNumerators(), ...s(this, k).getNumerators())), this));
+    a(this, "simplify", () => (this.multiply(U.lcm(...s(this, q).getDenominators(), ...s(this, M).getDenominators())), this.divide(U.gcd(...s(this, q).getNumerators(), ...s(this, M).getNumerators())), this));
     // -----------------------------------------------
-    a(this, "solve", () => new Nt(this.clone()).solve());
+    a(this, "solve", () => new At(this.clone()).solve());
     a(this, "test", (e) => this.left.evaluate(e).isEqual(this.right.evaluate(e)));
-    p(this, Dt, (e) => {
+    m(this, Dt, (e) => {
       if (e.includes("geq"))
         return e.includes("\\geq") ? "\\geq" : "geq";
       if (e.includes("leq"))
@@ -2046,10 +2094,10 @@ const Ee = class Ee {
     });
     // -----------------------------------------------
     // Equations solving algorithms
-    p(this, _e, (e) => e === void 0 ? "=" : e.includes("geq") || e.includes(">=") || e.includes("=>") ? ">=" : e.includes(">") ? ">" : e.includes("leq") || e.includes("<=") || e.includes("=<") ? "<=" : e.includes("<") ? "<" : "=");
-    p(this, Vt, () => s(this, se) === "=" ? this : s(this, se).includes("<") ? (s(this, se).replace("<", ">"), this) : s(this, se).includes(">") ? (s(this, se).replace(">", "<"), this) : this);
-    if (h(this, A, new O().zero()), h(this, k, new O().zero()), h(this, se, "="), e !== void 0 && t === void 0) {
-      if (e instanceof Ee)
+    m(this, et, (e) => e === void 0 ? "=" : e.includes("geq") || e.includes(">=") || e.includes("=>") ? ">=" : e.includes(">") ? ">" : e.includes("leq") || e.includes("<=") || e.includes("=<") ? "<=" : e.includes("<") ? "<" : "=");
+    m(this, Zt, () => s(this, re) === "=" ? this : s(this, re).includes("<") ? (s(this, re).replace("<", ">"), this) : s(this, re).includes(">") ? (s(this, re).replace(">", "<"), this) : this);
+    if (h(this, q, new O().zero()), h(this, M, new O().zero()), h(this, re, "="), e !== void 0 && t === void 0) {
+      if (e instanceof Ne)
         return e.clone();
       typeof e == "string" && this.parse(e);
     } else e !== void 0 && t !== void 0 && (this.left = new O(e), this.right = new O(t));
@@ -2064,12 +2112,12 @@ const Ee = class Ee {
    * @param value | Polynom | Monom | Fraction | string | monom
    */
   add(e) {
-    if (e instanceof Ee)
-      return s(this, A).add(e.left), s(this, k).add(e.right), this;
-    if (typeof e == "string" && !Ee.isEquationString(e))
-      return this.add(new Ee(e));
+    if (e instanceof Ne)
+      return s(this, q).add(e.left), s(this, M).add(e.right), this;
+    if (typeof e == "string" && !Ne.isEquationString(e))
+      return this.add(new Ne(e));
     const t = new O(e);
-    return s(this, A).add(t), s(this, k).add(t), this;
+    return s(this, q).add(t), s(this, M).add(t), this;
   }
   /**
    * Create an Equation using two polynoms.
@@ -2078,29 +2126,29 @@ const Ee = class Ee {
    * @param asNumeric
    */
   evaluate(e, t) {
-    const i = s(this, A).evaluate(e, t), n = s(this, k).evaluate(e, t);
+    const i = s(this, q).evaluate(e, t), n = s(this, M).evaluate(e, t);
     return t ? i === n : i.isEqual(n);
   }
   isEqual(e) {
-    const t = new Ee(e);
-    return t.left.isEqual(s(this, A)) && t.right.isEqual(s(this, k));
+    const t = new Ne(e);
+    return t.left.isEqual(s(this, q)) && t.right.isEqual(s(this, M));
   }
   pow(e) {
-    return s(this, A).pow(e), s(this, k).pow(e), this;
+    return s(this, q).pow(e), s(this, M).pow(e), this;
   }
   reduce() {
-    return this.moveLeft(), s(this, A).reduce(), this.simplify(), s(this, A).monoms[0].coefficient.isNegative() && this.multiply(-1), this;
+    return this.moveLeft(), s(this, q).reduce(), this.simplify(), s(this, q).monoms[0].coefficient.isNegative() && this.multiply(-1), this;
   }
   split() {
-    return [s(this, A).clone(), s(this, k).clone()];
+    return [s(this, q).clone(), s(this, M).clone()];
   }
   subtract(e) {
-    if (e instanceof Ee)
-      return s(this, A).subtract(e.left), s(this, k).subtract(e.right), this;
-    if (typeof e == "string" && !Ee.isEquationString(e))
-      return this.subtract(new Ee(e));
+    if (e instanceof Ne)
+      return s(this, q).subtract(e.left), s(this, M).subtract(e.right), this;
+    if (typeof e == "string" && !Ne.isEquationString(e))
+      return this.subtract(new Ne(e));
     const t = new O(e);
-    return s(this, A).subtract(t), s(this, k).subtract(t), this;
+    return s(this, q).subtract(t), s(this, M).subtract(t), this;
   }
   static isEquationString(e) {
     return e.includes("=") || e.includes("<") || e.includes(">") || e.includes("<=") || e.includes(">=");
@@ -2110,58 +2158,51 @@ const Ee = class Ee {
     return t === !0 && n.sort((r, l) => r.value - l.value), n;
   }
   get display() {
-    return `${s(this, A).display}${this.signAsTex}${s(this, k).display}`;
+    return `${s(this, q).display}${this.signAsTex}${s(this, M).display}`;
   }
   // Getter and setter
   get left() {
-    return s(this, A);
+    return s(this, q);
   }
   set left(e) {
-    h(this, A, e);
+    h(this, q, e);
   }
   get numberOfVars() {
     return this.variables.length;
   }
   get right() {
-    return s(this, k);
+    return s(this, M);
   }
   set right(e) {
-    h(this, k, e);
+    h(this, M, e);
   }
   // ------------------------------------------
   get sign() {
-    return s(this, se);
+    return s(this, re);
   }
   set sign(e) {
-    h(this, se, s(this, _e).call(this, e));
+    h(this, re, s(this, et).call(this, e));
   }
   get signAsTex() {
-    return s(this, se) === ">=" ? "\\geq" : s(this, se) === "<=" ? "\\leq" : s(this, se);
+    return s(this, re) === ">=" ? "\\geq" : s(this, re) === "<=" ? "\\leq" : s(this, re);
   }
   get tex() {
-    return `${s(this, A).tex}${this.signAsTex}${s(this, k).tex}`;
+    return `${s(this, q).tex}${this.signAsTex}${s(this, M).tex}`;
   }
   get variables() {
-    return [...new Set(s(this, k).variables.concat(s(this, A).variables))];
+    return [...new Set(s(this, M).variables.concat(s(this, q).variables))];
   }
 };
-A = new WeakMap(), k = new WeakMap(), se = new WeakMap(), Dt = new WeakMap(), _e = new WeakMap(), Vt = new WeakMap();
-let H = Ee;
-var Ne, Te, Ae, et;
+q = new WeakMap(), M = new WeakMap(), re = new WeakMap(), Dt = new WeakMap(), et = new WeakMap(), Zt = new WeakMap();
+let H = Ne;
+var Te, Ae, be, tt;
 const ve = class ve {
   constructor(e, t) {
-    p(this, Ne);
-    p(this, Te);
-    p(this, Ae);
-    p(this, et, !1);
-    if (e instanceof ve)
-      h(this, Te, e.polynom.clone()), h(this, Ae, e.power.clone());
-    else if (typeof e == "string" && t === void 0) {
-      const [i, n = "1"] = e.split("^");
-      h(this, Te, new O(i)), h(this, Ae, new u(n.replace("(", "").replace(")", "")));
-    } else
-      h(this, Te, new O(e)), h(this, Ae, new u(t ?? 1));
-    return h(this, Ne, 1), this;
+    m(this, Te);
+    m(this, Ae);
+    m(this, be);
+    m(this, tt, !1);
+    return e instanceof ve ? (h(this, Ae, e.polynom.clone()), h(this, be, e.power.clone()), t !== void 0 && s(this, be).multiply(new u(t))) : e !== void 0 ? (h(this, Ae, new O(e)), h(this, be, new u(t ?? 1))) : (h(this, Ae, new O()), h(this, be, new u(1))), h(this, Te, 1), this;
   }
   parse() {
     throw new Error("Method not implemented.");
@@ -2169,27 +2210,30 @@ const ve = class ve {
   clone() {
     return new ve(this);
   }
+  fromPolynom(e) {
+    return h(this, Ae, new O(e)), h(this, be, new u(1)), this;
+  }
   get tex() {
     const e = this.power.numerator, t = this.power.denominator;
     let i, n;
-    return s(this, Ne) === 0 && t > 1 ? (i = `\\sqrt${t === 2 ? "" : `[ ${t} ]`}{ ${this.polynom.tex} }`, n = e === 1 ? "" : `^{ ${e} }`) : (i = s(this, et) && this.power.isOne() ? this.polynom.tex : wi(this.polynom.tex), n = t === 1 && e === 1 ? "" : `^{ ${this.power.tex} }`), i = `${i}${n}`, s(this, Ne) === 0 && e < 0 && (i = `\\frac{ 1 }{ ${i} }`), i;
+    return s(this, Te) === 0 && t > 1 ? (i = `\\sqrt${t === 2 ? "" : `[ ${t} ]`}{ ${this.polynom.tex} }`, n = e === 1 ? "" : `^{ ${e} }`) : (i = s(this, tt) && this.power.isOne() ? this.polynom.tex : yi(this.polynom.tex), n = t === 1 && e === 1 ? "" : `^{ ${this.power.tex} }`), i = `${i}${n}`, s(this, Te) === 0 && e < 0 && (i = `\\frac{ 1 }{ ${i} }`), i;
   }
   get display() {
     const e = this.power.numerator, t = this.power.denominator;
     let i, n;
-    return s(this, Ne) === 0 && t > 1 ? (i = `${t === 2 ? "sqrt" : `root(${t})`}(${this.polynom.display})`, n = e === 1 ? "" : `^(${e})`) : (i = s(this, et) && this.power.isOne() ? this.polynom.display : wi(this.polynom.display, !1), n = t === 1 && e === 1 ? "" : `^(${this.power.display})`), i = `${i}${n}`, s(this, Ne) === 0 && e < 0 && (i = `1/(${i})`), i;
+    return s(this, Te) === 0 && t > 1 ? (i = `${t === 2 ? "sqrt" : `root(${t})`}(${this.polynom.display})`, n = e === 1 ? "" : `^(${e})`) : (i = s(this, tt) && this.power.isOne() ? this.polynom.display : yi(this.polynom.display, !1), n = t === 1 && e === 1 ? "" : `^(${this.power.display})`), i = `${i}${n}`, s(this, Te) === 0 && e < 0 && (i = `1/(${i})`), i;
   }
   add() {
     throw new Error("Adding two factors is not possible");
   }
   get asSingle() {
-    return h(this, et, !0), this;
+    return h(this, tt, !0), this;
   }
   degree(e) {
     return this.polynom.degree(e).multiply(this.power);
   }
   derivative() {
-    return this.power.isZero() ? [new ve("0", "1")] : this.power.isOne() ? [new ve(this.polynom.clone().derivative())] : [
+    return this.power.isZero() ? [new ve("0")] : this.power.isOne() ? [new ve(this.polynom.clone().derivative())] : [
       new ve(this.power.clone()),
       new ve(this.polynom.clone().derivative()),
       new ve(this.polynom.clone(), this.power.clone().subtract(1))
@@ -2239,25 +2283,25 @@ const ve = class ve {
     throw new Error("The two factors must be the same");
   }
   one() {
-    return s(this, Te).one(), s(this, Ae).one(), this;
+    return s(this, Ae).one(), s(this, be).one(), this;
   }
   opposite() {
     throw new Error("Method not implemented.");
   }
   get polynom() {
-    return s(this, Te);
+    return s(this, Ae);
   }
   set polynom(e) {
-    h(this, Te, e);
+    h(this, Ae, e);
   }
   pow(e) {
     return this.power.multiply(e), this;
   }
   get power() {
-    return s(this, Ae);
+    return s(this, be);
   }
   set power(e) {
-    h(this, Ae, new u(e));
+    h(this, be, new u(e));
   }
   primitive() {
     throw new Error("Method not implemented.");
@@ -2274,48 +2318,48 @@ const ve = class ve {
   subtract() {
     throw new Error("Subtracting two factors is not possible");
   }
-  tableOfSigns(e) {
-    const t = this.power.clone().reduce(), i = this.polynom.tableOfSigns(e);
-    return t.isStrictlyNegative() && (i.signs = Re(i.signs, "z", "d")), t.denominator % 2 === 0 ? i.signs = Re(i.signs, "-", "h") : t.numerator % 2 === 0 && (i.signs = Re(i.signs, "-", "+")), { roots: i.roots, signs: i.signs };
+  tableOfSigns() {
+    const e = this.power.clone().reduce(), t = this.polynom.tableOfSigns();
+    return e.isStrictlyNegative() && (t.signs = Et(t.signs, "z", "d")), e.denominator % 2 === 0 ? t.signs = Et(t.signs, "-", "h") : e.numerator % 2 === 0 && (t.signs = Et(t.signs, "-", "+")), { roots: t.roots, signs: t.signs };
   }
   get variables() {
     return this.polynom.variables;
   }
   get withPower() {
-    return h(this, Ne, 1), this;
+    return h(this, Te, 1), this;
   }
   get withRoot() {
-    return h(this, Ne, 0), this;
+    return h(this, Te, 0), this;
   }
   zero() {
-    return s(this, Te).zero(), s(this, Ae).one(), this;
+    return s(this, Ae).zero(), s(this, be).one(), this;
   }
 };
-Ne = new WeakMap(), Te = new WeakMap(), Ae = new WeakMap(), et = new WeakMap();
-let le = ve;
-var at = /* @__PURE__ */ ((o) => (o[o.ROOT = 0] = "ROOT", o[o.POWER = 1] = "POWER", o))(at || {}), D, $e, Zt, Ft;
-const He = class He {
+Te = new WeakMap(), Ae = new WeakMap(), be = new WeakMap(), tt = new WeakMap();
+let ae = ve;
+var lt = /* @__PURE__ */ ((o) => (o[o.ROOT = 0] = "ROOT", o[o.POWER = 1] = "POWER", o))(lt || {}), D, $e, Ft, jt;
+const Xe = class Xe {
   constructor(...e) {
-    p(this, D);
+    m(this, D);
     // Determine the letters in the linear system, usually ['x', 'y']
-    p(this, $e);
-    a(this, "parse", (...e) => (h(this, D, e.map((t) => new H(t))), s(this, Zt).call(this), this));
-    a(this, "clone", () => new He().parse(...s(this, D).map((e) => e.clone())));
+    m(this, $e);
+    a(this, "parse", (...e) => (h(this, D, e.map((t) => new H(t))), s(this, Ft).call(this), this));
+    a(this, "clone", () => new Xe().parse(...s(this, D).map((e) => e.clone())));
     a(this, "buildTex", (e, t) => {
       let i, n, r = [];
       const l = [];
-      for (const d of e)
-        r = r.concat(d.letters());
+      for (const f of e)
+        r = r.concat(f.letters());
       r = [...new Set(r)], r.sort();
-      for (let d = 0; d < e.length; d++) {
-        const f = e[d];
+      for (let f = 0; f < e.length; f++) {
+        const d = e[f];
         i = [];
-        for (const g of r)
-          n = f.left.monomByLetter(g), i.length === 0 ? i.push(n.isZero() ? "" : n.tex) : i.push(n.isZero() ? "" : (n.coefficient.sign() === 1 ? "+" : "") + n.tex);
-        if (i.push("="), i.push(f.right.tex), (t == null ? void 0 : t[d]) !== void 0) {
+        for (const p of r)
+          n = d.left.monomByLetter(p), i.length === 0 ? i.push(n.isZero() ? "" : n.tex) : i.push(n.isZero() ? "" : (n.coefficient.sign() === 1 ? "+" : "") + n.tex);
+        if (i.push("="), i.push(d.right.tex), (t == null ? void 0 : t[f]) !== void 0) {
           i[i.length - 1] = i[i.length - 1] + " \\phantom{\\quad}";
-          for (const g of t[d])
-            i.push(`\\ \\cdot\\ ${g.startsWith("-") ? "\\left(" + g + "\\right)" : g}`);
+          for (const p of t[f])
+            i.push(`\\ \\cdot\\ ${p.startsWith("-") ? "\\left(" + p + "\\right)" : p}`);
         }
         l.push(i.join("&"));
       }
@@ -2341,16 +2385,16 @@ const He = class He {
           if (l === n)
             continue;
           const c = i[l][n].clone().opposite();
-          for (let d = 0; d < i[l].length; d++)
-            i[l][d].add(i[n][d].clone().multiply(c));
-          if (i[l].slice(0, i[l].length - 1).every((d) => d.isZero()))
+          for (let f = 0; f < i[l].length; f++)
+            i[l][f].add(i[n][f].clone().multiply(c));
+          if (i[l].slice(0, i[l].length - 1).every((f) => f.isZero()))
             return i[l][i[l].length - 1].isZero() ? [new u().infinite()] : [];
         }
       }
       return i.map((n) => n[n.length - 1]);
     });
-    p(this, Zt, () => (h(this, $e, s(this, D).reduce((e, t) => [.../* @__PURE__ */ new Set([...e, ...t.variables])], [])), s(this, $e).sort(), this));
-    p(this, Ft, () => {
+    m(this, Ft, () => (h(this, $e, s(this, D).reduce((e, t) => [.../* @__PURE__ */ new Set([...e, ...t.variables])], [])), s(this, $e).sort(), this));
+    m(this, jt, () => {
       const e = [], t = [];
       for (const i of s(this, D)) {
         const n = [], r = i.clone().reorder();
@@ -2369,7 +2413,7 @@ const He = class He {
     if (e.some((r) => r.length !== i))
       throw new Error("All rows must have the same number of columns");
     const n = t.split("").splice(0, i - 1);
-    return new He(
+    return new Xe(
       ...e.map((r) => {
         const l = new O(n.join(""), ...r);
         return new H(l, 0);
@@ -2377,7 +2421,7 @@ const He = class He {
     );
   }
   add(e, t) {
-    if (e instanceof He) {
+    if (e instanceof Xe) {
       const i = e.equations.length;
       if (i !== s(this, D).length)
         throw new Error("The number of equations must be the same");
@@ -2417,7 +2461,7 @@ const He = class He {
     return this.variables.length === s(this, D).length;
   }
   get matrix() {
-    return s(this, Ft).call(this);
+    return s(this, jt).call(this);
   }
   multiply(e, t) {
     if (Array.isArray(e)) {
@@ -2438,7 +2482,7 @@ const He = class He {
     return [];
   }
   subtract(e, t) {
-    if (e instanceof He) {
+    if (e instanceof Xe) {
       const i = e.equations.length;
       if (i !== s(this, D).length)
         throw new Error("The number of equations must be the same");
@@ -2464,18 +2508,18 @@ const He = class He {
     t.sort(), h(this, $e, t);
   }
 };
-D = new WeakMap(), $e = new WeakMap(), Zt = new WeakMap(), Ft = new WeakMap();
-let ii = He;
-var Oe, pt, si;
-class es {
+D = new WeakMap(), $e = new WeakMap(), Ft = new WeakMap(), jt = new WeakMap();
+let ti = Xe;
+var Oe, mt, ii;
+class ss {
   /**
    *
    * @param {string} value (optional) Default polynom to parse on class creation
    */
   constructor(e) {
-    p(this, pt);
-    p(this, Oe);
-    a(this, "parse", (e) => (h(this, Oe, new Qt(Fe.SET).parse(e).rpn), this));
+    m(this, mt);
+    m(this, Oe);
+    a(this, "parse", (e) => (h(this, Oe, new Kt(Ze.SET).parse(e).rpn), this));
     return h(this, Oe, []), e !== void 0 && this.parse(e), this;
   }
   evaluate(e) {
@@ -2552,20 +2596,20 @@ class es {
     return s(this, Oe).filter((e) => e.tokenType === "variable").map((e) => e.token);
   }
   vennAB() {
-    return q(this, pt, si).call(this, {
+    return A(this, mt, ii).call(this, {
       A: ["A", "AB"],
       B: ["B", "AB"]
     }, ["A", "B", "AB", "E"]);
   }
   vennABC() {
-    return q(this, pt, si).call(this, {
+    return A(this, mt, ii).call(this, {
       A: ["A", "AB", "AC", "ABC"],
       B: ["B", "AB", "BC", "ABC"],
       C: ["C", "AC", "BC", "ABC"]
     }, ["A", "B", "C", "AB", "AC", "BC", "ABC", "E"]);
   }
 }
-Oe = new WeakMap(), pt = new WeakSet(), si = function(e, t) {
+Oe = new WeakMap(), mt = new WeakSet(), ii = function(e, t) {
   const i = [];
   let n;
   if (t === void 0) {
@@ -2585,7 +2629,7 @@ Oe = new WeakMap(), pt = new WeakSet(), si = function(e, t) {
         case "&":
           if (i.length >= 2) {
             const l = i.pop(), c = i.pop();
-            c && l && i.push(new Set([...c].filter((d) => l.has(d))));
+            c && l && i.push(new Set([...c].filter((f) => l.has(f))));
           }
           break;
         case "|":
@@ -2597,7 +2641,7 @@ Oe = new WeakMap(), pt = new WeakSet(), si = function(e, t) {
         case "-":
           if (i.length >= 2) {
             const l = i.pop(), c = i.pop();
-            c && l && i.push(new Set([...c].filter((d) => !l.has(d))));
+            c && l && i.push(new Set([...c].filter((f) => !l.has(f))));
           }
           break;
         case "!":
@@ -2609,111 +2653,152 @@ Oe = new WeakMap(), pt = new WeakSet(), si = function(e, t) {
       }
   return [...i[0]].sort();
 };
-var tt, E, mt, ri, gt, oi;
-const ue = class ue {
+var it, N, We, ni, ki, gt, ri;
+const X = class X {
   constructor(...e) {
-    p(this, gt);
-    p(this, tt, at.POWER);
-    p(this, E, []);
+    m(this, gt);
+    m(this, it, lt.POWER);
+    m(this, N, []);
     return this.parse(...e), this;
   }
   parse(...e) {
-    return e.length === 0 ? this : (h(this, E, []), e.forEach((t) => {
-      if (typeof t == "string") {
-        const i = t.split(")(").join(")*(").split("*");
-        s(this, E).push(...i.map((n) => new le(n)));
-      } else t instanceof ue ? s(this, E).push(...t.factors.map((i) => i.clone())) : s(this, E).push(new le(t));
+    return h(this, N, []), e.length === 0 ? this : (e.forEach((t) => {
+      t instanceof X ? s(this, N).push(...t.factors.map((i) => i.clone())) : s(this, N).push(new ae(t));
     }), this);
   }
   clone() {
-    return new ue(...s(this, E).map((e) => e.clone()));
+    return new X(...s(this, N).map((e) => e.clone()));
   }
   get tex() {
-    const { num: e, den: t } = q(this, gt, oi).call(this);
+    const { num: e, den: t } = A(this, gt, ri).call(this);
     if (t.length === 0)
       return e.length === 1 ? e[0].asSingle.tex : e.map((r) => r.tex).join("");
     const i = e.length === 1 ? e[0].asSingle.tex : e.map((r) => r.tex).join(""), n = t.length === 1 ? t[0].asSingle.tex : t.map((r) => r.tex).join("");
     return `\\frac{ ${i} }{ ${n} }`;
   }
   get display() {
-    const { num: e, den: t } = q(this, gt, oi).call(this);
+    const { num: e, den: t } = A(this, gt, ri).call(this);
     if (t.length === 0)
-      return e.length === 1 ? e[0].asSingle.display : e.map((r) => r.display).join("");
+      return e.length === 1 ? e[0].asSingle.display : e.map(
+        (r, l) => l === 0 && r.polynom.monoms.length === 1 ? r.asSingle.display : r.display
+      ).join("");
     const i = e.length === 1 ? e[0].asSingle.display : e.map((r) => r.display).join(""), n = t.length === 1 ? t[0].asSingle.display : t.map((r) => r.display).join("");
     return `(${i})/(${n})`;
   }
   static gcd(...e) {
     var i;
     if (e.length === 0)
-      return new ue().one();
+      return new X().one();
     if (e.length === 1)
       return e[0];
     if (e.length === 2)
-      return q(i = ue, mt, ri).call(i, e[0], e[1]);
+      return A(i = X, We, ni).call(i, e[0], e[1]);
     let t = e[0];
     return e.shift(), e.forEach((n) => {
       var r;
-      return t = q(r = ue, mt, ri).call(r, t, n);
+      return t = A(r = X, We, ni).call(r, t, n);
+    }), t;
+  }
+  static lcm(...e) {
+    if (e.length === 0)
+      return new X().one();
+    if (e.length === 1)
+      return e[0];
+    let t = e[0];
+    return e.shift(), e.forEach((i) => {
+      var n;
+      return t = A(n = X, We, ki).call(n, t, i);
     }), t;
   }
   add(...e) {
-    let t = [this, ...e];
-    const i = ue.gcd(...t);
-    t = t.map((r) => r.divide(i).reduce());
-    const n = new O("0");
-    return t.forEach((r) => n.add(r.develop())), h(this, E, [
-      ...i.factors,
-      new le(n)
-    ]), this;
+    const t = [this.numerator, ...e.map((c) => c.numerator)], i = [this.denominator, ...e.map((c) => c.denominator)];
+    let n;
+    if (i.some((c) => c.factors.length > 0)) {
+      const c = X.lcm(...i);
+      t.forEach((f, d) => {
+        f.multiply(c.clone().divide(i[d]));
+      }), n = c;
+    }
+    const r = X.gcd(...t), l = new O(0).add(
+      ...t.map((c) => c.divide(r).reduce().develop().factors[0].polynom)
+    ).reduce();
+    return h(this, N, [
+      ...r.factors,
+      new ae(l)
+    ]), n && this.divide(n), h(this, N, s(this, N).filter((c) => !c.power.isZero())), this;
   }
   get asPower() {
-    return h(this, tt, at.POWER), this;
+    return h(this, it, lt.POWER), this;
   }
   get asRoot() {
-    return h(this, tt, at.ROOT), this;
+    return h(this, it, lt.ROOT), this;
   }
   degree(e) {
-    return s(this, E).reduce((t, i) => t.add(i.degree(e)), new u("0"));
+    return s(this, N).reduce((t, i) => t.add(i.degree(e)), new u("0"));
   }
   get denominator() {
-    return new ue(...s(this, E).filter((e) => e.power.isNegative()));
+    return new X(
+      ...s(this, N).filter((e) => e.power.isNegative()).map((e) => e.clone().inverse())
+    );
   }
   derivative() {
-    const e = [], t = s(this, E).length;
+    const e = [], t = s(this, N).length;
     for (let n = 0; n < t; n++) {
-      const r = s(this, E).slice(), l = r.splice(n, 1)[0];
-      e.push(new ue(...r).multiply(new ue(...l.derivative())));
+      const r = s(this, N).slice(), l = r.splice(n, 1)[0].derivative();
+      e.push(
+        new X(...r, ...l)
+      );
     }
     e.forEach((n) => n.reduce());
     const i = e.shift();
-    return i !== void 0 && h(this, E, i.factors), this.add(...e);
+    return i !== void 0 && h(this, N, i.factors), this.add(...e);
   }
   develop() {
-    const e = new O("1");
-    return s(this, E).forEach((t) => {
-      e.multiply(t.develop());
-    }), e;
+    const e = new O("1"), t = new O("1");
+    return this.numerator.factors.forEach((i) => {
+      e.multiply(i.develop());
+    }), this.denominator.factors.forEach((i) => {
+      t.multiply(i.develop());
+    }), new X().fromPolynom(e, t);
   }
   divide(e) {
-    return h(this, E, s(this, E).concat(e.clone().factors.map((t) => t.inverse()))), this;
+    return h(this, N, s(this, N).concat(e.clone().factors.map((t) => t.inverse()))), this;
   }
   evaluate(e, t) {
-    return t ? s(this, E).reduce((i, n) => i * n.evaluate(e, t), 1) : s(this, E).reduce((i, n) => i.multiply(n.evaluate(e)), new u("1"));
+    return t ? s(this, N).reduce((i, n) => i * n.evaluate(e, t), 1) : s(this, N).reduce((i, n) => i.multiply(n.evaluate(e)), new u("1"));
+  }
+  factorize(e) {
+    const t = [];
+    s(this, N).forEach((l) => {
+      const c = l.polynom.factorize(e);
+      if (c.length > 1) {
+        const f = l.power.clone();
+        t.push(...c.map((d) => new ae(d, f)));
+      } else
+        t.push(l.clone());
+    });
+    const i = new X(...t), n = i.numerator.reduce(), r = i.denominator.reduce();
+    return n.divide(r);
   }
   get factors() {
-    return s(this, E);
+    return s(this, N);
   }
   set factors(e) {
-    h(this, E, e);
+    h(this, N, e);
   }
-  fromPolynom(e, t, i) {
-    return h(this, E, new O(e).factorize(i).map((n) => new le(n))), t && new O(t).factorize(i).forEach((n) => s(this, E).push(new le(n, -1))), this;
-  }
-  getFactors() {
-    return s(this, E);
+  fromPolynom(e, t) {
+    if (h(this, N, [new ae(new O(e))]), t) {
+      const i = new O(t);
+      if (i.isOne())
+        return this;
+      if (i.isZero())
+        throw new Error("Cannot divide by zero");
+      s(this, N).push(new ae(i, -1));
+    }
+    return this;
   }
   getZeroes() {
-    const e = [].concat(...s(this, E).map((t) => t.polynom.getZeroes()));
+    const e = [].concat(...s(this, N).map((t) => t.polynom.getZeroes()));
     return e.sort((t, i) => t.value - i.value), e.filter(
       (t, i, n) => i === n.findIndex(
         (r) => r.value === t.value
@@ -2721,113 +2806,146 @@ const ue = class ue {
     );
   }
   hasVariable(e) {
-    return s(this, E).some((t) => t.hasVariable(e));
+    return s(this, N).some((t) => t.hasVariable(e));
   }
   inverse() {
-    return h(this, E, s(this, E).map((e) => e.inverse())), this;
+    return h(this, N, s(this, N).map((e) => e.inverse())), this;
   }
   isEqual(e) {
-    const t = ue.gcd(this, e), i = this.clone().divide(t).reduce(), n = e.clone().divide(t).reduce();
+    const t = X.gcd(this, e), i = this.clone().divide(t).reduce(), n = e.clone().divide(t).reduce();
     return i.isOne() && n.isOne();
   }
   isOne() {
-    return s(this, E).every((e) => e.isOne());
+    return s(this, N).every((e) => e.isOne());
   }
   isZero() {
-    return s(this, E).every((e) => e.isZero());
+    return s(this, N).every((e) => e.isZero());
   }
   multiply(...e) {
     return e.forEach((t) => {
-      h(this, E, s(this, E).concat(t.clone().factors));
+      h(this, N, s(this, N).concat(t.clone().factors));
     }), this;
   }
   get numerator() {
-    return new ue(...s(this, E).filter((e) => e.power.isPositive()));
+    return new X(...s(this, N).filter((e) => e.power.isPositive()));
   }
   one() {
-    return h(this, E, [new le("1", "1")]), this;
+    return h(this, N, [new ae("1", "1")]), this;
   }
   opposite() {
-    const e = s(this, E).findIndex((t) => t.display === "(-1)");
-    return e >= 0 ? s(this, E).splice(e, 1) : s(this, E).push(new le("-1", "1")), this;
+    const e = s(this, N).findIndex((t) => t.display === "(-1)");
+    return e >= 0 ? s(this, N).splice(e, 1) : s(this, N).push(new ae("-1", "1")), this;
   }
   pow(e) {
-    return h(this, E, s(this, E).map((t) => t.pow(e))), this;
+    return h(this, N, s(this, N).map((t) => t.pow(e))), this;
   }
   primitive() {
     throw new Error("Method not implemented.");
   }
   reduce() {
-    const e = ei(this);
-    return h(this, E, Object.values(e).map((t) => {
+    const e = ht(this);
+    return h(this, N, Object.values(e).map((t) => {
       const i = t[0].polynom, n = t.reduce((r, l) => r.add(l.power), new u("0"));
-      return new le(i, n.reduce());
+      return new ae(i, n.reduce());
     }).filter((t) => !t.power.isZero())), this;
   }
   root(e) {
-    return h(this, E, s(this, E).map((t) => t.root(e))), this;
+    return h(this, N, s(this, N).map((t) => t.root(e))), this;
   }
-  sort() {
-    return h(this, E, s(this, E).sort((e, t) => e.degree().isLeq(t.degree()) ? -1 : 1)), this;
+  /**
+   * Reoarder the factors using :
+   * 1. number of monoms
+   * 2. degree of polynom
+   * 3. power of polyfactor
+   */
+  sort(e) {
+    return s(this, N).sort((t, i) => {
+      const n = t.power.value, r = i.power.value;
+      if (n * r < 0)
+        return -n;
+      const l = t.polynom.monoms.length, c = i.polynom.monoms.length;
+      if (l !== c)
+        return l - c;
+      const f = t.polynom.degree(e).value, d = i.polynom.degree(e).value;
+      return f !== d ? f - d : n !== r ? n - r : t.degree().isLeq(i.degree()) ? -1 : 1;
+    }), this;
   }
   sqrt() {
-    return h(this, E, s(this, E).map((e) => e.sqrt())), this;
+    return h(this, N, s(this, N).map((e) => e.sqrt())), this;
   }
   subtract(...e) {
     return this.add(...e.map((t) => t.opposite()));
   }
   tableOfSigns() {
-    const e = this.getZeroes(), t = this.factors.map((n) => ({ factor: new le(n), ...n.tableOfSigns(e) }));
-    return { signs: t.map((n) => n.signs).reduce((n, r) => (n.length === 0 ? n = r : r.forEach((l, c) => {
-      switch (l) {
+    const e = this.getZeroes(), t = e.map((r) => r.value), i = this.factors.map((r) => ({ factor: new ae(r), ...r.tableOfSigns() }));
+    return i.forEach((r) => {
+      const l = new Array(2 * e.length + 1).fill("");
+      let c = r.signs.shift(), f = r.roots.shift();
+      const d = l.map((p, w) => {
+        if (w % 2 === 0)
+          return c;
+        if (f === void 0 || f.value !== t[(w - 1) / 2])
+          return "t";
+        const V = r.signs.shift();
+        return c = r.signs.shift(), f = r.roots.shift(), V;
+      });
+      r.roots = e, r.signs = d;
+    }), { signs: i.map((r) => r.signs).reduce((r, l) => r.length === 0 ? l : (l.forEach((c, f) => {
+      switch (c) {
         case "d":
-          n[c] = "d";
+          r[f] = "d";
           break;
         case "z":
-          n[c] = n[c] === "d" ? "d" : "z";
+          r[f] = r[f] === "d" ? "d" : "z";
           break;
         case "h":
-          n[c] = "h";
+          r[f] = "h";
           break;
         case "-":
-          n[c] = n[c] === "h" ? "h" : n[c] === "-" ? "+" : "-";
+          r[f] = r[f] === "h" ? "h" : r[f] === "-" ? "+" : "-";
           break;
       }
-    }), n), []), roots: e, factors: t };
+    }), r), []), roots: e, factors: i };
   }
   get variables() {
-    return s(this, E).reduce((e, t) => e.concat(t.variables), []);
+    return s(this, N).reduce((e, t) => e.concat(t.variables), []);
   }
   zero() {
-    return h(this, E, [new le("0", "1")]), this;
+    return h(this, N, [new ae("0", "1")]), this;
   }
 };
-tt = new WeakMap(), E = new WeakMap(), mt = new WeakSet(), ri = function(e, t) {
-  const i = ei(e), n = ei(t), l = Object.keys(i).filter((c) => Object.hasOwn(n, c)).map((c) => {
-    const d = i[c].reduce((g, N) => g.add(N.power), new u("0")), f = n[c].reduce((g, N) => g.add(N.power), new u("0"));
-    return new le(c, u.min(d, f));
+it = new WeakMap(), N = new WeakMap(), We = new WeakSet(), ni = function(e, t) {
+  const i = ht(e), n = ht(t), l = Object.keys(i).filter((c) => Object.hasOwn(n, c)).map((c) => {
+    const f = i[c].reduce((p, w) => p.add(w.power), new u("0")), d = n[c].reduce((p, w) => p.add(w.power), new u("0"));
+    return new ae(c, u.min(f, d));
   });
-  return new ue(...l);
-}, gt = new WeakSet(), oi = function() {
+  return new X(...l);
+}, ki = function(e, t) {
+  const i = ht(e), n = ht(t), l = [.../* @__PURE__ */ new Set([...Object.keys(i), ...Object.keys(n)])].map((c) => {
+    const f = Object.hasOwn(i, c) ? i[c].reduce((p, w) => p.add(w.power), new u("0")) : new u(0), d = Object.hasOwn(n, c) ? n[c].reduce((p, w) => p.add(w.power), new u("0")) : new u(0);
+    return new ae(c, u.max(f, d));
+  });
+  return new X(...l);
+}, gt = new WeakSet(), ri = function() {
   let e, t = [];
-  return s(this, tt) === at.ROOT ? (e = this.numerator.factors, t = this.denominator.factors.map((i) => i.clone().inverse())) : e = s(this, E), e.length === 0 && (e = [new le("1")]), { num: e, den: t };
-}, p(ue, mt);
-let ni = ue;
-function ei(o) {
+  return s(this, it) === lt.ROOT ? (e = this.numerator.factors, t = this.denominator.factors) : e = s(this, N), e.length === 0 && (e = [new ae("1")]), { num: e, den: t };
+}, m(X, We);
+let si = X;
+function ht(o) {
   const e = new u().one(), t = o.factors.reduce((i, n) => {
     if (n.polynom.degree().isZero())
       return n.polynom.monoms.length > 0 && e.multiply(n.polynom.monoms[0].coefficient), i;
     const r = n.polynom.display;
     return Object.hasOwn(i, r) ? i[r].push(n) : i[r] = [n], i;
   }, {});
-  return e.isOne() || (t[e.display] = [new le(e.display, 1)]), t;
+  return e.isOne() || (t[e.display] = [new ae(e.display, 1)]), t;
 }
-function ts(o, e) {
+function ns(o, e) {
   return o.dimension === e.dimension && o.array.every(
     (t, i) => e.array[i].isEqual(t)
   );
 }
-function is(o, e) {
+function rs(o, e) {
   if (o.dimension !== e.dimension)
     return !1;
   const t = e.array[0].value / o.array[0].value;
@@ -2835,14 +2953,18 @@ function is(o, e) {
     (i, n) => e.array[n].value === i.value * t
   );
 }
-function ss(o, e) {
+function os(o, e) {
   return o.dimension !== e.dimension ? new u().invalid() : o.array.reduce(
     (t, i, n) => t.add(i.clone().multiply(e.array[n])),
     new u(0)
   );
 }
-function ns(...o) {
-  return o.some((e) => e.dimension !== o[0].dimension) ? new u().invalid() : o[0].dimension === 2 && o.length !== 2 ? new u().invalid() : o[0].dimension === 3 && o.length !== 3 ? new u().invalid() : o[0].dimension === 2 ? o[0].array[0].clone().multiply(o[1].array[1]).subtract(o[0].array[1].clone().multiply(o[1].array[0])) : o[0].array[0].clone().multiply(
+function hs(...o) {
+  if (o.some((e) => e.dimension !== o[0].dimension))
+    throw new Error("All vectors must have the same dimension");
+  if (o[0].dimension !== o.length)
+    throw new Error(`The determinant of dimension ${o[0].dimension} must have the same number of vectors (${o.length} given)`);
+  return o[0].dimension === 2 ? o[0].array[0].clone().multiply(o[1].array[1]).subtract(o[0].array[1].clone().multiply(o[1].array[0])) : o[0].array[0].clone().multiply(
     o[1].array[1].clone().multiply(o[2].array[2]).subtract(o[1].array[2].clone().multiply(o[2].array[1]))
   ).subtract(
     o[0].array[1].clone().multiply(
@@ -2850,11 +2972,11 @@ function ns(...o) {
     )
   ).add(o[0].array[2].clone().multiply(o[1].array[0].clone().multiply(o[2].array[1]).subtract(o[1].array[1].clone().multiply(o[2].array[0]))));
 }
-var F, De;
+var F, Le;
 const Se = class Se {
   constructor(...e) {
-    p(this, F, []);
-    p(this, De, !1);
+    m(this, F, []);
+    m(this, Le, !1);
     a(this, "zero", () => (s(this, F).forEach((e) => e.zero()), this));
     a(this, "one", () => (this.zero(), this.x.one(), this));
     a(this, "opposite", () => (s(this, F).forEach((e) => e.opposite()), this));
@@ -2864,15 +2986,15 @@ const Se = class Se {
       const e = this.norm;
       return e === 0 ? this : this.divideByScalar(e);
     });
-    a(this, "dot", (e) => ss(this, e));
+    a(this, "dot", (e) => os(this, e));
     a(this, "normal", () => {
       if (this.dimension >= 3)
         throw new Error("Normal vector can only be determined in 2D");
       const e = this.x.clone().opposite(), t = this.y.clone();
       return s(this, F)[0] = t, s(this, F)[1] = e, this;
     });
-    a(this, "isEqual", (e) => ts(this, e));
-    a(this, "isColinearTo", (e) => is(this, e));
+    a(this, "isEqual", (e) => ns(this, e));
+    a(this, "isColinearTo", (e) => rs(this, e));
     a(this, "isNormalTo", (e) => this.dot(e).isZero());
     a(this, "multiplyByScalar", (e) => {
       const t = new u(e);
@@ -2880,9 +3002,9 @@ const Se = class Se {
     });
     a(this, "divideByScalar", (e) => this.multiplyByScalar(new u(e).inverse()));
     a(this, "simplify", () => this.multiplyByScalar(
-      G.lcm(...this.array.map((e) => e.denominator))
+      U.lcm(...this.array.map((e) => e.denominator))
     ).divideByScalar(
-      G.gcd(...this.array.map((e) => e.numerator))
+      U.gcd(...this.array.map((e) => e.numerator))
     ).multiplyByScalar(
       this.x.isNegative() ? -1 : 1
     ));
@@ -2929,10 +3051,10 @@ const Se = class Se {
     s(this, F)[2] = new u(e);
   }
   get asPoint() {
-    return s(this, De);
+    return s(this, Le);
   }
   set asPoint(e) {
-    h(this, De, e);
+    h(this, Le, e);
   }
   get normSquare() {
     return this.array.reduce((e, t) => e.add(t.clone().pow(2)), new u(0));
@@ -2941,10 +3063,10 @@ const Se = class Se {
     return Math.sqrt(this.normSquare.value);
   }
   get tex() {
-    return s(this, De) ? `\\left(${this.array.map((e) => e.tex).join(";")}\\right)` : `\\begin{pmatrix} ${this.array.map((e) => e.tex).join(" \\\\ ")} \\end{pmatrix}`;
+    return s(this, Le) ? `\\left(${this.array.map((e) => e.tex).join(";")}\\right)` : `\\begin{pmatrix} ${this.array.map((e) => e.tex).join(" \\\\ ")} \\end{pmatrix}`;
   }
   get display() {
-    return s(this, De) ? `(${this.array.map((e) => e.display).join(";")})` : `((${this.array.map((e) => e.display).join(",")}))`;
+    return s(this, Le) ? `(${this.array.map((e) => e.display).join(";")})` : `((${this.array.map((e) => e.display).join(",")}))`;
   }
   setDimension(e = 2) {
     if (e < 2)
@@ -2972,7 +3094,7 @@ const Se = class Se {
     return `((${e.join(",")}))`;
   }
   defineAsPoint(e) {
-    return h(this, De, e !== !1), this;
+    return h(this, Le, e !== !1), this;
   }
   parse(...e) {
     if (e.length === 0)
@@ -3035,40 +3157,40 @@ const Se = class Se {
     };
   }
 };
-F = new WeakMap(), De = new WeakMap();
-let x = Se;
-var we = /* @__PURE__ */ ((o) => (o.None = "none", o.Parallel = "parallel", o.Perpendicular = "perpendicular", o.Tangent = "tangent", o))(we || {}), hi = /* @__PURE__ */ ((o) => (o.None = "none", o.Parallel = "parallel", o.Perpendicular = "perpendicular", o.Tangent = "tangent", o))(hi || {});
-function qi(o = 0.5) {
+F = new WeakMap(), Le = new WeakMap();
+let E = Se;
+var we = /* @__PURE__ */ ((o) => (o.None = "none", o.Parallel = "parallel", o.Perpendicular = "perpendicular", o.Tangent = "tangent", o))(we || {}), oi = /* @__PURE__ */ ((o) => (o.None = "none", o.Parallel = "parallel", o.Perpendicular = "perpendicular", o.Tangent = "tangent", o))(oi || {});
+function Ci(o = 0.5) {
   return Math.random() < o;
 }
-function de(o, e, t) {
+function pe(o, e, t) {
   if (e === void 0)
-    return o >= 0 ? de(0, o) : de(o, 0);
+    return o >= 0 ? pe(0, o) : pe(o, 0);
   if (o === e)
     return o;
   if (t === void 0)
     return Math.floor(Math.random() * (e - o + 1) + o);
   if (Math.abs(e - o) <= t.length)
     throw new Error("The number of excluded values is too high.");
-  let i = de(o, e);
+  let i = pe(o, e);
   for (; t.includes(i); )
-    i = de(o, e);
+    i = pe(o, e);
   return i;
 }
-function Y(o, e) {
-  return e === !1 ? qi() ? de(1, o) : -de(1, o) : de(-o, o);
+function Q(o, e) {
+  return e === !1 ? Ci() ? pe(1, o) : -pe(1, o) : pe(-o, o);
 }
-function rs(o) {
-  let e = G.primes();
-  return o !== void 0 && (e = e.filter((t) => t < o)), pi(e);
+function as(o) {
+  let e = U.primes();
+  return o !== void 0 && (e = e.filter((t) => t < o)), di(e);
 }
-function os(o, e) {
-  return e === void 0 && (e = 1), o.length <= 0 ? Object.values(o) : ki(o).slice(0, e);
+function ls(o, e) {
+  return e === void 0 && (e = 1), o.length <= 0 ? Object.values(o) : $i(o).slice(0, e);
 }
-function pi(o) {
-  return o.length === 0 ? null : o[de(0, o.length - 1)];
+function di(o) {
+  return o.length === 0 ? null : o[pe(0, o.length - 1)];
 }
-function ki(o) {
+function $i(o) {
   const e = Object.values(o);
   for (let t = e.length - 1; t > 0; t--) {
     const i = Math.floor(Math.random() * (t + 1)), n = e[t];
@@ -3076,19 +3198,19 @@ function ki(o) {
   }
   return e;
 }
-class B extends x {
+class B extends E {
   constructor(...e) {
     super(), e.length > 0 && this.parse(...e);
   }
   parse(...e) {
     if (this.asPoint = !0, e.length === 1) {
-      if (e[0] instanceof x)
+      if (e[0] instanceof E)
         return this.array = e[0].copy(), this;
       if (typeof e[0] == "string")
         return this.fromString(e[0]), this;
     }
     if (e.length > 1) {
-      if (e.some((i) => i instanceof x))
+      if (e.some((i) => i instanceof E))
         throw new Error("Creating a point with  multiple argument requires an input fraction");
       const t = e.map((i) => new u(i));
       if (t.some((i) => i.isNaN()))
@@ -3102,32 +3224,32 @@ class B extends x {
     return e.array = this.copy(), e.asPoint = !0, e;
   }
 }
-var Ie, S, z, U, ne, Q, Pe, be;
-const Ze = class Ze {
+var Pe, S, z, W, oe, J, Ie, xe;
+const De = class De {
   /**
    * Value can be a mix of:
    *
    * @param values
    */
   constructor(...e) {
-    p(this, Ie);
+    m(this, Pe);
     // ax + by + c = 0
-    p(this, S);
-    p(this, z);
-    p(this, U);
-    p(this, ne);
-    p(this, Q);
-    p(this, Pe);
-    p(this, be, "canonical");
+    m(this, S);
+    m(this, z);
+    m(this, W);
+    m(this, oe);
+    m(this, J);
+    m(this, Ie);
+    m(this, xe, "canonical");
     a(this, "randomPoint", (e) => {
-      const t = s(this, Q).clone().multiplyByScalar(Y(e === void 0 || e <= 1 ? 3 : e, !1)).add(s(this, ne));
+      const t = s(this, J).clone().multiplyByScalar(Q(e === void 0 || e <= 1 ? 3 : e, !1)).add(s(this, oe));
       return t.asPoint = !0, t;
     });
     a(this, "randomNearPoint", (e) => {
       const t = this.randomPoint(e);
       let i = 10;
       for (; this.isOnLine(t) && i > 0; )
-        t.x.add(Y(1, !1)), t.y.add(Y(1, !1)), i--;
+        t.x.add(Q(1, !1)), t.y.add(Q(1, !1)), i--;
       return t;
     });
     // ------------------------------------------
@@ -3142,7 +3264,7 @@ const Ze = class Ze {
       if (e.length === 0)
         return this;
       if (e.length === 1) {
-        if (e[0] instanceof Ze)
+        if (e[0] instanceof De)
           return this.fromCoefficient(e[0].a, e[0].b, e[0].c);
         if (e[0] instanceof H)
           return this.fromEquation(e[0]);
@@ -3154,21 +3276,21 @@ const Ze = class Ze {
             return this;
           }
       }
-      if (e.length === 2 && e.every((t) => t instanceof x)) {
+      if (e.length === 2 && e.every((t) => t instanceof E)) {
         const t = e;
         if (t[0].asPoint && t[1].asPoint)
-          return this.fromPointAndDirection(t[0], new x(t[0], t[1]));
+          return this.fromPointAndDirection(t[0], new E(t[0], t[1]));
         if (t[0].asPoint && !t[1].asPoint)
           return this.fromPointAndDirection(t[0], t[1]);
       }
       if (e.length === 3) {
-        if (e[0] instanceof x && e[1] instanceof x) {
+        if (e[0] instanceof E && e[1] instanceof E) {
           if (e[2] === we.Perpendicular)
             return this.fromPointAndNormal(e[0], e[1]);
           if (e[2] === we.Parallel)
             return this.fromPointAndDirection(e[0], e[1]);
         }
-        return e[0] instanceof x && e[1] instanceof Ze ? e[2] === we.Parallel || e[2] === null ? this.fromPointAndLine(e[0], e[1], we.Parallel) : this.fromPointAndLine(e[0], e[1], we.Perpendicular) : this.fromCoefficient(
+        return e[0] instanceof E && e[1] instanceof De ? e[2] === we.Parallel || e[2] === null ? this.fromPointAndLine(e[0], e[1], we.Parallel) : this.fromPointAndLine(e[0], e[1], we.Perpendicular) : this.fromCoefficient(
           e[0],
           e[1],
           e[2]
@@ -3189,42 +3311,42 @@ const Ze = class Ze {
         e.left.monomByDegree(0).coefficient
       );
     });
-    a(this, "fromCoefficient", (e, t, i) => (h(this, S, new u(e)), h(this, z, new u(t)), h(this, U, new u(i)), h(this, Q, new x(s(this, z).clone(), s(this, S).clone().opposite())), h(this, ne, new x(new u().zero(), s(this, U).clone())), h(this, Pe, s(this, Q).clone().normal()), this));
+    a(this, "fromCoefficient", (e, t, i) => (h(this, S, new u(e)), h(this, z, new u(t)), h(this, W, new u(i)), h(this, J, new E(s(this, z).clone(), s(this, S).clone().opposite())), h(this, oe, new E(new u().zero(), s(this, W).clone())), h(this, Ie, s(this, J).clone().normal()), this));
     a(this, "fromPointAndDirection", (e, t) => (this.fromCoefficient(
       t.y,
       t.x.clone().opposite(),
       e.x.clone().multiply(t.y).subtract(e.y.clone().multiply(t.x)).opposite()
-    ), h(this, ne, e.clone()), h(this, Q, t.clone()), h(this, Pe, s(this, Q).clone().normal()), this));
+    ), h(this, oe, e.clone()), h(this, J, t.clone()), h(this, Ie, s(this, J).clone().normal()), this));
     a(this, "fromPointAndNormal", (e, t) => this.fromCoefficient(
       t.x,
       t.y,
       e.x.clone().multiply(t.x).add(e.y.clone().multiply(t.y)).opposite()
     ));
     a(this, "fromPointAndLine", (e, t, i) => (i === void 0 && (i = we.Parallel), i === we.Parallel ? this.fromPointAndNormal(e, t.normal) : i === we.Perpendicular ? this.fromPointAndNormal(e, t.director) : this));
-    a(this, "clone", () => (h(this, S, s(this, S).clone()), h(this, z, s(this, z).clone()), h(this, U, s(this, U).clone()), h(this, Q, s(this, Q).clone()), h(this, ne, s(this, ne).clone()), h(this, Pe, s(this, Pe).clone()), this));
+    a(this, "clone", () => (h(this, S, s(this, S).clone()), h(this, z, s(this, z).clone()), h(this, W, s(this, W).clone()), h(this, J, s(this, J).clone()), h(this, oe, s(this, oe).clone()), h(this, Ie, s(this, Ie).clone()), this));
     // ------------------------------------------
     // Mathematical operations
     // ------------------------------------------
     a(this, "isOnLine", (e) => s(this, S).clone().multiply(e.x).add(
       s(this, z).clone().multiply(e.y)
-    ).add(s(this, U)).isZero());
+    ).add(s(this, W)).isZero());
     a(this, "isParallelTo", (e) => this.slope.isEqual(e.slope) && this.height.isNotEqual(e.height));
     a(this, "isSameAs", (e) => this.slope.isEqual(e.slope) && this.height.isEqual(e.height));
     a(this, "isPerpendicularTo", (e) => this.d.isNormalTo(e.d));
     a(this, "isVertical", () => this.slope.isInfinity());
     a(this, "simplify", () => {
-      const e = G.lcm(s(this, S).denominator, s(this, z).denominator, s(this, U).denominator), t = G.gcd(s(this, S).numerator, s(this, z).numerator, s(this, U).numerator);
+      const e = U.lcm(s(this, S).denominator, s(this, z).denominator, s(this, W).denominator), t = U.gcd(s(this, S).numerator, s(this, z).numerator, s(this, W).numerator);
       return this.fromCoefficient(
         s(this, S).clone().multiply(e).divide(t),
         s(this, z).clone().multiply(e).divide(t),
-        s(this, U).clone().multiply(e).divide(t)
+        s(this, W).clone().multiply(e).divide(t)
       ), this;
     });
-    a(this, "simplifyDirection", () => (s(this, Q).simplify(), this));
+    a(this, "simplifyDirection", () => (s(this, J).simplify(), this));
     a(this, "intersection", (e) => {
       const t = new B();
       let i = !1, n = !1;
-      return s(this, z).isZero() || e.b.isZero(), this.isParallelTo(e) ? (t.x = new u().invalid(), t.y = new u().invalid(), i = !0) : this.isSameAs(e) ? (t.x = new u().invalid(), t.y = new u().invalid(), n = !0) : (t.x = s(this, z).clone().multiply(e.c).subtract(s(this, U).clone().multiply(e.b)).divide(s(this, S).clone().multiply(e.b).subtract(s(this, z).clone().multiply(e.a))), t.y = s(this, S).clone().multiply(e.c).subtract(s(this, U).clone().multiply(e.a)).divide(s(this, z).clone().multiply(e.a).subtract(s(this, S).clone().multiply(e.b)))), {
+      return s(this, z).isZero() || e.b.isZero(), this.isParallelTo(e) ? (t.x = new u().invalid(), t.y = new u().invalid(), i = !0) : this.isSameAs(e) ? (t.x = new u().invalid(), t.y = new u().invalid(), n = !0) : (t.x = s(this, z).clone().multiply(e.c).subtract(s(this, W).clone().multiply(e.b)).divide(s(this, S).clone().multiply(e.b).subtract(s(this, z).clone().multiply(e.a))), t.y = s(this, S).clone().multiply(e.c).subtract(s(this, W).clone().multiply(e.a)).divide(s(this, z).clone().multiply(e.a).subtract(s(this, S).clone().multiply(e.b)))), {
         point: t,
         hasIntersection: !(i || n),
         isParallel: i,
@@ -3239,7 +3361,7 @@ const Ze = class Ze {
       const t = this.getEquation().isolate("x"), i = new u(e);
       return t instanceof H ? t.right.evaluate({ y: i }) : new u().invalid();
     });
-    return h(this, S, new u().zero()), h(this, z, new u().zero()), h(this, U, new u().zero()), h(this, ne, new x()), h(this, Q, new x()), h(this, Pe, new x()), h(this, Ie, !0), e.length > 0 && this.parse(...e), this;
+    return h(this, S, new u().zero()), h(this, z, new u().zero()), h(this, W, new u().zero()), h(this, oe, new E()), h(this, J, new E()), h(this, Ie, new E()), h(this, Pe, !0), e.length > 0 && this.parse(...e), this;
   }
   get a() {
     return s(this, S);
@@ -3256,30 +3378,30 @@ const Ze = class Ze {
     h(this, z, e);
   }
   get c() {
-    return s(this, U);
+    return s(this, W);
   }
   set c(e) {
-    h(this, U, e);
+    h(this, W, e);
   }
   get OA() {
-    return s(this, ne);
+    return s(this, oe);
   }
   set OA(e) {
-    h(this, ne, e);
+    h(this, oe, e);
   }
   get d() {
-    return s(this, Q);
+    return s(this, J);
   }
   set d(e) {
-    h(this, Q, e);
+    h(this, J, e);
   }
   get n() {
-    return s(this, Pe);
+    return s(this, Ie);
   }
   // ------------------------------------------
   getEquation() {
-    const e = new H(new O().parse("xy", s(this, S), s(this, z), s(this, U)), new O("0"));
-    return s(this, Ie) ? e.simplify() : e;
+    const e = new H(new O().parse("xy", s(this, S), s(this, z), s(this, W)), new O("0"));
+    return s(this, Pe) ? e.simplify() : e;
   }
   // get system(): { x: Equation, y: Equation } {
   //     const e1 = new Equation(
@@ -3295,33 +3417,33 @@ const Ze = class Ze {
   //     return { x: e1, y: e2 }
   // }
   get canonical() {
-    return h(this, be, "canonical"), this;
+    return h(this, xe, "canonical"), this;
   }
   get equation() {
-    return h(this, be, "equation"), this;
+    return h(this, xe, "equation"), this;
   }
   get mxh() {
-    return h(this, be, "mxh"), this;
+    return h(this, xe, "mxh"), this;
   }
   get parametric() {
-    return h(this, be, "parametric"), this;
+    return h(this, xe, "parametric"), this;
   }
   get system() {
-    return h(this, be, "system"), this;
+    return h(this, xe, "system"), this;
   }
   get tex() {
-    const e = s(this, be);
-    switch (h(this, be, "canonical"), e) {
+    const e = s(this, xe);
+    switch (h(this, xe, "canonical"), e) {
       case "equation":
         return this.getEquation().reorder().tex;
       case "mxh":
         return this.slope.isInfinity() ? "x=" + this.OA.x.tex : "y=" + new O().parse("x", this.slope, this.height).tex;
       case "parametric":
       case "system": {
-        const t = s(this, Q).clone();
-        return s(this, Ie) && t.simplify(), e === "parametric" ? `${x.asTex("x", "y")} = ${x.asTex(s(this, ne).x.tex, s(this, ne).y.tex)} + k\\cdot ${x.asTex(t.x.tex, t.y.tex)}` : `\\left\\{\\begin{aligned}
-            x &= ${new O(s(this, ne).x).add(new M(s(this, Q).x).multiply(new M("k"))).reorder("k", !0).tex}\\\\ 
-            y &= ${new O(s(this, ne).y).add(new M(s(this, Q).y).multiply(new M("k"))).reorder("k", !0).tex}
+        const t = s(this, J).clone();
+        return s(this, Pe) && t.simplify(), e === "parametric" ? `${E.asTex("x", "y")} = ${E.asTex(s(this, oe).x.tex, s(this, oe).y.tex)} + k\\cdot ${E.asTex(t.x.tex, t.y.tex)}` : `\\left\\{\\begin{aligned}
+            x &= ${new O(s(this, oe).x).add(new k(s(this, J).x).multiply(new k("k"))).reorder("k", !0).tex}\\\\ 
+            y &= ${new O(s(this, oe).y).add(new k(s(this, J).y).multiply(new k("k"))).reorder("k", !0).tex}
             \\end{aligned}\\right.`;
       }
       default: {
@@ -3331,21 +3453,21 @@ const Ze = class Ze {
     }
   }
   get reduceBeforeDisplay() {
-    return s(this, Ie);
+    return s(this, Pe);
   }
   set reduceBeforeDisplay(e) {
-    h(this, Ie, e);
+    h(this, Pe, e);
   }
   get display() {
-    const e = s(this, be);
-    switch (h(this, be, "canonical"), e) {
+    const e = s(this, xe);
+    switch (h(this, xe, "canonical"), e) {
       case "equation":
         return this.getEquation().reorder().display;
       case "mxh":
         return this.slope.isInfinity() ? "x=" + this.OA.x.display : "y=" + new O().parse("x", this.slope, this.height).display;
       case "parametric": {
-        const t = s(this, Q).clone();
-        return s(this, Ie) && t.simplify(), `((x,y))=((${s(this, ne).x.display},${s(this, ne).y.display}))+k((${t.x.display},${t.y.display}))`;
+        const t = s(this, J).clone();
+        return s(this, Pe) && t.simplify(), `((x,y))=((${s(this, oe).x.display},${s(this, oe).y.display}))+k((${t.x.display},${t.y.display}))`;
       }
       default: {
         const t = this.getEquation();
@@ -3354,22 +3476,22 @@ const Ze = class Ze {
     }
   }
   get normal() {
-    return new x(s(this, S), s(this, z));
+    return new E(s(this, S), s(this, z));
   }
   get director() {
-    return s(this, Q).clone();
+    return s(this, J).clone();
   }
   get slope() {
     return s(this, S).clone().opposite().divide(s(this, z));
   }
   get height() {
-    return s(this, U).clone().opposite().divide(s(this, z));
+    return s(this, W).clone().opposite().divide(s(this, z));
   }
   fromPoints(e, t) {
-    return this.fromPointAndDirection(e, new x(e, t));
+    return this.fromPointAndDirection(e, new E(e, t));
   }
   distanceTo(e) {
-    const t = e.x.clone().multiply(s(this, S)).add(e.y.clone().multiply(s(this, z))).add(s(this, U)).abs(), i = this.normal.normSquare;
+    const t = e.x.clone().multiply(s(this, S)).add(e.y.clone().multiply(s(this, z))).add(s(this, W)).abs(), i = this.normal.normSquare;
     if (i.isZero())
       return {
         value: NaN,
@@ -3389,7 +3511,7 @@ const Ze = class Ze {
   }
   hitSegment(e, t) {
     const i = this.intersection(
-      new Ze().fromPoints(e, t)
+      new De().fromPoints(e, t)
     );
     return i.hasIntersection ? i.point.x.value >= Math.min(e.x.value, t.x.value) && i.point.x.value <= Math.max(e.x.value, t.x.value) && i.point.y.value >= Math.min(e.y.value, t.y.value) && i.point.y.value <= Math.max(e.y.value, t.y.value) : !1;
   }
@@ -3399,28 +3521,28 @@ const Ze = class Ze {
   canonicalAsFloatCoefficient(e) {
     e === void 0 && (e = 2);
     let t = "";
-    return s(this, S).isZero() || (s(this, S).isOne() ? t = "x" : s(this, S).clone().opposite().isOne() ? t = "-x" : t = s(this, S).value.toFixed(e) + "x"), s(this, z).isZero() || (s(this, z).isPositive() && (t += "+"), t += s(this, z).value.toFixed(e) + "y"), s(this, U).isZero() || (s(this, U).isPositive() && (t += "+"), t += s(this, U).value.toFixed(e)), t + "=0";
+    return s(this, S).isZero() || (s(this, S).isOne() ? t = "x" : s(this, S).clone().opposite().isOne() ? t = "-x" : t = s(this, S).value.toFixed(e) + "x"), s(this, z).isZero() || (s(this, z).isPositive() && (t += "+"), t += s(this, z).value.toFixed(e) + "y"), s(this, W).isZero() || (s(this, W).isPositive() && (t += "+"), t += s(this, W).value.toFixed(e)), t + "=0";
   }
 };
-Ie = new WeakMap(), S = new WeakMap(), z = new WeakMap(), U = new WeakMap(), ne = new WeakMap(), Q = new WeakMap(), Pe = new WeakMap(), be = new WeakMap(), // A line is defined as the canonical form
-a(Ze, "PERPENDICULAR", we.Perpendicular), a(Ze, "PARALLEL", we.Parallel);
-let R = Ze;
-var re, j, qe, jt, Ut, Gt, oe, Mi, Et, Ci, $i, Ii, ai;
-const Wt = class Wt {
+Pe = new WeakMap(), S = new WeakMap(), z = new WeakMap(), W = new WeakMap(), oe = new WeakMap(), J = new WeakMap(), Ie = new WeakMap(), xe = new WeakMap(), // A line is defined as the canonical form
+a(De, "PERPENDICULAR", we.Perpendicular), a(De, "PARALLEL", we.Parallel);
+let R = De;
+var he, j, qe, Ut, Wt, Gt, le, Pi, Nt, Ii, Bi, Si, hi;
+const Ht = class Ht {
   constructor(...e) {
-    p(this, oe);
-    p(this, re);
-    p(this, j);
-    p(this, qe);
+    m(this, le);
+    m(this, he);
+    m(this, j);
+    m(this, qe);
     /**
      * Get the relative position between circle and line. It corresponds to the number of intersection.
      * @param {Line} L
      * @returns {number}
      */
     a(this, "relativePosition", (e) => {
-      if (s(this, re) === void 0 || s(this, j) === void 0)
+      if (s(this, he) === void 0 || s(this, j) === void 0)
         throw new Error("Circle not defined");
-      const t = e.distanceTo(s(this, re)), i = Math.sqrt(s(this, j).value);
+      const t = e.distanceTo(s(this, he)), i = Math.sqrt(s(this, j).value);
       return t.value - i > 1e-10 ? 0 : Math.abs(t.value - i) < 1e-10 ? 1 : 2;
     });
     a(this, "lineIntersection", (e) => {
@@ -3430,13 +3552,13 @@ const Wt = class Wt {
       const i = s(this, qe).clone(), n = e.getEquation().clone().isolate("x"), r = e.getEquation().clone().isolate("y");
       return n instanceof H && r instanceof H && (i.replaceBy("y", r.right).simplify(), i.solve()), t;
     });
-    a(this, "tangents", (e) => e instanceof u ? s(this, Gt).call(this, e) : this.isPointOnCircle(e) ? s(this, jt).call(this, e) : s(this, re) !== void 0 && s(this, re).distanceTo(e).value > this.radius.value ? s(this, Ut).call(this, e) : (console.log("No tangents as the point is inside !"), []));
+    a(this, "tangents", (e) => e instanceof u ? s(this, Gt).call(this, e) : this.isPointOnCircle(e) ? s(this, Ut).call(this, e) : s(this, he) !== void 0 && s(this, he).distanceTo(e).value > this.radius.value ? s(this, Wt).call(this, e) : (console.log("No tangents as the point is inside !"), []));
     a(this, "isPointOnCircle", (e) => {
       var t;
       return ((t = s(this, qe)) == null ? void 0 : t.test({ x: e.x, y: e.y })) ?? !1;
     });
     a(this, "getPointsOnCircle", (e) => {
-      const t = G.pythagoreanTripletsWithTarget(this.squareRadius.value, !0), i = [];
+      const t = U.pythagoreanTripletsWithTarget(this.squareRadius.value, !0), i = [];
       return t.forEach((n) => {
         for (const r of [[1, 1], [-1, 1], [-1, -1], [1, -1]])
           i.push(
@@ -3447,26 +3569,26 @@ const Wt = class Wt {
           );
       }), i;
     });
-    p(this, jt, (e) => {
-      const t = new x(this.center, e);
+    m(this, Ut, (e) => {
+      const t = new E(this.center, e);
       return [new R(e, t, we.Perpendicular)];
     });
-    p(this, Ut, (e) => {
+    m(this, Wt, (e) => {
       const t = this.center.x.clone().subtract(e.x), i = this.center.y.clone().subtract(e.y), n = new O("x"), r = new O("x^2+1");
-      return n.multiply(t).subtract(i).pow(2), r.multiply(this.squareRadius), new H(n, r).solve().map((d) => {
-        let f;
-        const g = new H("y", "x");
-        return d.exact instanceof u ? (f = e.x.clone().opposite().multiply(d.exact).add(e.y), g.right.multiply(d.exact).add(f)) : (f = e.x.clone().opposite().multiply(d.value).add(e.y), g.right.multiply(d.value).add(f)), new R(g);
+      return n.multiply(t).subtract(i).pow(2), r.multiply(this.squareRadius), new H(n, r).solve().map((f) => {
+        let d;
+        const p = new H("y", "x");
+        return f.exact instanceof u ? (d = e.x.clone().opposite().multiply(f.exact).add(e.y), p.right.multiply(f.exact).add(d)) : (d = e.x.clone().opposite().multiply(f.value).add(e.y), p.right.multiply(f.value).add(d)), new R(p);
       });
     });
-    p(this, Gt, (e) => {
-      const t = e.numerator, i = -e.denominator, n = this.center.x.clone(), r = this.center.y.clone(), l = this.squareRadius.clone().multiply(e.numerator ** 2 + e.denominator ** 2), c = n.clone().multiply(t).opposite().subtract(r.clone().multiply(i)).add(l.clone().sqrt()), d = n.clone().multiply(t).opposite().subtract(r.clone().multiply(i)).subtract(l.clone().sqrt());
-      return [new R(t, i, c), new R(t, i, d)];
+    m(this, Gt, (e) => {
+      const t = e.numerator, i = -e.denominator, n = this.center.x.clone(), r = this.center.y.clone(), l = this.squareRadius.clone().multiply(e.numerator ** 2 + e.denominator ** 2), c = n.clone().multiply(t).opposite().subtract(r.clone().multiply(i)).add(l.clone().sqrt()), f = n.clone().multiply(t).opposite().subtract(r.clone().multiply(i)).subtract(l.clone().sqrt());
+      return [new R(t, i, c), new R(t, i, f)];
     });
     e.length > 0 && this.parse(...e);
   }
   get center() {
-    return s(this, re) ?? new B();
+    return s(this, he) ?? new B();
   }
   get squareRadius() {
     return s(this, j) ?? new u(0);
@@ -3499,17 +3621,17 @@ const Wt = class Wt {
     return this.center.x.isZero() ? e = "x^2" : e = `(x${this.center.x.isNegative() ? "+" : "-"}${this.center.x.clone().abs().tex})^2`, this.center.y.isZero() ? t = "y^2" : t = `(y${this.center.y.isNegative() ? "+" : "-"}${this.center.y.clone().abs().tex})^2`, `${e}+${t}=${this.squareRadius.display}`;
   }
   clone() {
-    return new Wt(
+    return new Ht(
       this.center.clone(),
       this.squareRadius.clone(),
       !0
     );
   }
   setRadius(e, t) {
-    return t ? h(this, j, new u(e)) : h(this, j, new u(e).pow(2)), q(this, oe, Et).call(this), this;
+    return t ? h(this, j, new u(e)) : h(this, j, new u(e).pow(2)), A(this, le, Nt).call(this), this;
   }
   parse(...e) {
-    return q(this, oe, Mi).call(this), typeof e[0] == "string" ? q(this, oe, ai).call(this, new H(e[0])) : e[0] instanceof H ? q(this, oe, ai).call(this, e[0]) : e[0] instanceof Wt ? q(this, oe, Ci).call(this, e[0]) : e[0] instanceof B && e.length > 1 && (e[1] instanceof B ? e[2] instanceof B || q(this, oe, Ii).call(this, e[0], e[1]) : (e[1] instanceof u || typeof e[1] == "number") && q(this, oe, $i).call(this, e[0], e[1], typeof e[2] == "boolean" ? e[2] : !1)), q(this, oe, Et).call(this), this;
+    return A(this, le, Pi).call(this), typeof e[0] == "string" ? A(this, le, hi).call(this, new H(e[0])) : e[0] instanceof H ? A(this, le, hi).call(this, e[0]) : e[0] instanceof Ht ? A(this, le, Ii).call(this, e[0]) : e[0] instanceof B && e.length > 1 && (e[1] instanceof B ? e[2] instanceof B || A(this, le, Si).call(this, e[0], e[1]) : (e[1] instanceof u || typeof e[1] == "number") && A(this, le, Bi).call(this, e[0], e[1], typeof e[2] == "boolean" ? e[2] : !1)), A(this, le, Nt).call(this), this;
   }
   // private _parseThroughtThreePoints(A: Point, B: Point, C: Point): this {
   //     const T = new Triangle(A, B, C), mAB = T.remarquables.mediators.AB.clone(),
@@ -3518,35 +3640,35 @@ const Wt = class Wt {
   //     return this
   // }
 };
-re = new WeakMap(), j = new WeakMap(), qe = new WeakMap(), jt = new WeakMap(), Ut = new WeakMap(), Gt = new WeakMap(), oe = new WeakSet(), Mi = function() {
-  return h(this, re, void 0), h(this, j, void 0), h(this, qe, void 0), this;
-}, Et = function() {
+he = new WeakMap(), j = new WeakMap(), qe = new WeakMap(), Ut = new WeakMap(), Wt = new WeakMap(), Gt = new WeakMap(), le = new WeakSet(), Pi = function() {
+  return h(this, he, void 0), h(this, j, void 0), h(this, qe, void 0), this;
+}, Nt = function() {
   h(this, qe, new H(
     new O(`(x-(${this.center.x.display}))^2+(y-(${this.center.y.display}))^2`),
     new O(this.squareRadius.display)
   ).moveLeft());
-}, Ci = function(e) {
-  return h(this, re, e.center.clone()), h(this, j, e.squareRadius.clone()), q(this, oe, Et).call(this), this;
-}, $i = function(e, t, i) {
-  return h(this, re, e.clone()), i ? h(this, j, new u(t)) : h(this, j, new u(t).pow(2)), this;
-}, Ii = function(e, t) {
-  return h(this, re, e.clone()), h(this, j, new x(s(this, re), t).normSquare), this;
-}, ai = function(e) {
+}, Ii = function(e) {
+  return h(this, he, e.center.clone()), h(this, j, e.squareRadius.clone()), A(this, le, Nt).call(this), this;
+}, Bi = function(e, t, i) {
+  return h(this, he, e.clone()), i ? h(this, j, new u(t)) : h(this, j, new u(t).pow(2)), this;
+}, Si = function(e, t) {
+  return h(this, he, e.clone()), h(this, j, new E(s(this, he), t).normSquare), this;
+}, hi = function(e) {
   if (e.moveLeft(), e.degree("x").value === 2 && e.degree("y").value === 2) {
     const t = e.left.monomByDegree(2, "x"), i = e.left.monomByDegree(2, "y");
     let n, r, l;
-    t.coefficient.isEqual(i.coefficient) ? (e.divide(t.coefficient), n = e.left.monomByDegree(1, "x"), r = e.left.monomByDegree(1, "y"), l = e.left.monomByDegree(0), h(this, re, new B(n.coefficient.clone().divide(2).opposite(), r.coefficient.clone().divide(2).opposite())), h(this, j, l.coefficient.clone().opposite().add(s(this, re).x.clone().pow(2)).add(s(this, re).y.clone().pow(2)))) : (h(this, re, void 0), h(this, j, void 0));
+    t.coefficient.isEqual(i.coefficient) ? (e.divide(t.coefficient), n = e.left.monomByDegree(1, "x"), r = e.left.monomByDegree(1, "y"), l = e.left.monomByDegree(0), h(this, he, new B(n.coefficient.clone().divide(2).opposite(), r.coefficient.clone().divide(2).opposite())), h(this, j, l.coefficient.clone().opposite().add(s(this, he).x.clone().pow(2)).add(s(this, he).y.clone().pow(2)))) : (h(this, he, void 0), h(this, j, void 0));
   }
   return this;
 };
-let Tt = Wt;
-var L, W;
-const lt = class lt {
+let Tt = Ht;
+var L, G;
+const ct = class ct {
   constructor(e, t) {
     // ax + by + c = 0
-    p(this, L, new B());
-    p(this, W, new x());
-    a(this, "clone", () => (h(this, W, s(this, W).clone()), h(this, L, s(this, L).clone()), this));
+    m(this, L, new B());
+    m(this, G, new E());
+    a(this, "clone", () => (h(this, G, s(this, G).clone()), h(this, L, s(this, L).clone()), this));
     // ------------------------------------------
     // Mathematical operations
     // ------------------------------------------
@@ -3586,14 +3708,14 @@ const lt = class lt {
     //     return new Fraction().invalid()
     // }
     a(this, "randomPoint", (e = 5) => {
-      const t = s(this, L).clone(), i = new u(Y(e, !1));
+      const t = s(this, L).clone(), i = new u(Q(e, !1));
       return new B(
-        t.x.clone().add(s(this, W).x.clone().multiply(i)),
-        t.y.clone().add(s(this, W).y.clone().multiply(i)),
-        t.z.clone().add(s(this, W).z.clone().multiply(i))
+        t.x.clone().add(s(this, G).x.clone().multiply(i)),
+        t.y.clone().add(s(this, G).y.clone().multiply(i)),
+        t.z.clone().add(s(this, G).z.clone().multiply(i))
       );
     });
-    return h(this, L, e.clone()), h(this, W, t.asPoint ? new x(e, t) : t.clone()), this;
+    return h(this, L, e.clone()), h(this, G, t.asPoint ? new E(e, t) : t.clone()), this;
   }
   get OA() {
     return s(this, L);
@@ -3605,18 +3727,18 @@ const lt = class lt {
     return s(this, L).clone();
   }
   get d() {
-    return s(this, W);
+    return s(this, G);
   }
   set d(e) {
-    h(this, W, e);
+    h(this, G, e);
   }
   get tex() {
     return {
-      parametric: `${x.asTex("x", "y", "z")} = ${x.asTex(s(this, L).x.tex, s(this, L).y.tex, s(this, L).z.tex)} + k\\cdot ${x.asTex(s(this, W).x.tex, s(this, W).y.tex, s(this, W).z.tex)}`,
+      parametric: `${E.asTex("x", "y", "z")} = ${E.asTex(s(this, L).x.tex, s(this, L).y.tex, s(this, L).z.tex)} + k\\cdot ${E.asTex(s(this, G).x.tex, s(this, G).y.tex, s(this, G).z.tex)}`,
       system: `\\left\\{\\begin{aligned}
-    x &= ${new O(s(this, L).x).add(new M(s(this, W).x).multiply(new M("k"))).reorder("k", !0).tex}\\\\ 
-    y &= ${new O(s(this, L).y).add(new M(s(this, W).y).multiply(new M("k"))).reorder("k", !0).tex}\\\\
-    z &= ${new O(s(this, L).z).add(new M(s(this, W).z).multiply(new M("k"))).reorder("k", !0).tex}
+    x &= ${new O(s(this, L).x).add(new k(s(this, G).x).multiply(new k("k"))).reorder("k", !0).tex}\\\\ 
+    y &= ${new O(s(this, L).y).add(new k(s(this, G).y).multiply(new k("k"))).reorder("k", !0).tex}\\\\
+    z &= ${new O(s(this, L).z).add(new k(s(this, G).z).multiply(new k("k"))).reorder("k", !0).tex}
 \\end{aligned}\\right.`,
       cartesian: `\\frac{ ${new O("x", 1, s(this, L).x.clone().opposite()).tex} }{ ${this.direction.x.tex} } = \\frac{ ${new O("y", 1, s(this, L).y.clone().opposite()).tex} }{ ${this.direction.y.tex} } = \\frac{ ${new O("z", 1, s(this, L).z.clone().opposite()).tex} }{ ${this.direction.z.tex} }`
     };
@@ -3624,16 +3746,16 @@ const lt = class lt {
   get display() {
     const e = s(this, L).x.display, t = s(this, L).y.display, i = s(this, L).z.display, n = this.direction.simplify(), r = n.x.display, l = n.y.display, c = n.z.display;
     return {
-      parametric: `${x.asDisplay("x", "y", "z")} = ${x.asDisplay(s(this, L).x.display, s(this, L).y.display, s(this, L).z.display)} + k\\cdot ${x.asDisplay(s(this, W).x.display, s(this, W).y.display, s(this, W).z.display)}`,
+      parametric: `${E.asDisplay("x", "y", "z")} = ${E.asDisplay(s(this, L).x.display, s(this, L).y.display, s(this, L).z.display)} + k\\cdot ${E.asDisplay(s(this, G).x.display, s(this, G).y.display, s(this, G).z.display)}`,
       system: "",
       cartesian: `(x-${e})/${r} = (y-${t})/${l} = (z-${i})/${c}`
     };
   }
   get direction() {
-    return s(this, W).clone();
+    return s(this, G).clone();
   }
   distanceTo(e) {
-    const t = new x(s(this, L), e), i = this.direction, n = this.direction.normSquare, r = t.cross(i).normSquare, l = r.clone().divide(n), c = l.clone().sqrt();
+    const t = new E(s(this, L), e), i = this.direction, n = this.direction.normSquare, r = t.cross(i).normSquare, l = r.clone().divide(n), c = l.clone().sqrt();
     return console.log("CROSS", t.cross(i).display), {
       value: Math.sqrt(l.value),
       fraction: l.clone().sqrt(),
@@ -3642,69 +3764,69 @@ const lt = class lt {
   }
   hitSegment(e, t) {
     const i = this.intersection(
-      new lt(e, t)
+      new ct(e, t)
     );
     return i.hasIntersection ? i.point.x.value >= Math.min(e.x.value, t.x.value) && i.point.x.value <= Math.max(e.x.value, t.x.value) && i.point.y.value >= Math.min(e.y.value, t.y.value) && i.point.y.value <= Math.max(e.y.value, t.y.value) && i.point.z.value >= Math.min(e.z.value, t.z.value) && i.point.z.value <= Math.max(e.z.value, t.z.value) : !1;
   }
 };
-L = new WeakMap(), W = new WeakMap(), // A line is defined as the canonical form
-a(lt, "PERPENDICULAR", hi.Perpendicular), a(lt, "PARALLEL", hi.Parallel);
-let At = lt;
-var ke;
-class hs {
+L = new WeakMap(), G = new WeakMap(), // A line is defined as the canonical form
+a(ct, "PERPENDICULAR", oi.Perpendicular), a(ct, "PARALLEL", oi.Parallel);
+let Ot = ct;
+var Me;
+class cs {
   constructor(...e) {
-    p(this, ke, []);
-    return h(this, ke, e), this;
+    m(this, Me, []);
+    return h(this, Me, e), this;
   }
   get values() {
-    return s(this, ke);
+    return s(this, Me);
   }
   get array() {
-    return s(this, ke).map((e) => e.array);
+    return s(this, Me).map((e) => e.array);
   }
   get dimension() {
-    return [s(this, ke).length, s(this, ke)[0].dimension];
+    return [s(this, Me).length, s(this, Me)[0].dimension];
   }
   isSquare() {
-    return s(this, ke).length === s(this, ke)[0].dimension;
+    return s(this, Me).length === s(this, Me)[0].dimension;
   }
   determinant() {
     if (!this.isSquare())
       throw new Error("Matrix is not square");
-    return ns(...this.values);
+    return hs(...this.values);
   }
 }
-ke = new WeakMap();
-var Me, Ge;
-const gi = class gi {
+Me = new WeakMap();
+var ke, Ue;
+const mi = class mi {
   constructor(e) {
-    p(this, Me, new x(0, 0, 1));
-    p(this, Ge, new B(0, 0, 0));
+    m(this, ke, new E(0, 0, 1));
+    m(this, Ue, new B(0, 0, 0));
     return e && this.parse(e), this;
   }
   get normal() {
-    return s(this, Me);
+    return s(this, ke);
   }
   set normal(e) {
-    h(this, Me, e), s(this, Me).asPoint = !1;
+    h(this, ke, e), s(this, ke).asPoint = !1;
   }
   get point() {
-    return s(this, Ge);
+    return s(this, Ue);
   }
   set point(e) {
-    h(this, Ge, e), s(this, Ge).asPoint = !0;
+    h(this, Ue, e), s(this, Ue).asPoint = !0;
   }
   get a() {
-    return s(this, Me).x;
+    return s(this, ke).x;
   }
   get b() {
-    return s(this, Me).y;
+    return s(this, ke).y;
   }
   get c() {
-    return s(this, Me).z;
+    return s(this, ke).z;
   }
   get d() {
-    return s(this, Me).dot(s(this, Ge)).opposite();
+    return s(this, ke).dot(s(this, Ue)).opposite();
   }
   get tex() {
     return new H(
@@ -3725,26 +3847,26 @@ const gi = class gi {
       return;
     }
     if (e.equation) {
-      const r = e.equation.moveLeft().reduce().left, l = r.monomByLetter("x").coefficient, c = r.monomByLetter("y").coefficient, d = r.monomByLetter("z").coefficient, f = r.monomByDegree(0).coefficient;
-      this.normal = new x(l, c, d), l.isNotZero() ? this.point = new B(f.clone().divide(l).opposite(), 0, 0) : c.isNotZero() ? this.point = new B(0, f.clone().divide(c).opposite(), 0) : this.point = new B(0, 0, f.clone().divide(d).opposite());
+      const r = e.equation.moveLeft().reduce().left, l = r.monomByLetter("x").coefficient, c = r.monomByLetter("y").coefficient, f = r.monomByLetter("z").coefficient, d = r.monomByDegree(0).coefficient;
+      this.normal = new E(l, c, f), l.isNotZero() ? this.point = new B(d.clone().divide(l).opposite(), 0, 0) : c.isNotZero() ? this.point = new B(0, d.clone().divide(c).opposite(), 0) : this.point = new B(0, 0, d.clone().divide(f).opposite());
       return;
     }
-    if (((i = e.points) == null ? void 0 : i.length) === 3 && e.points.every((r) => r instanceof x)) {
-      const r = e.points[0], l = e.points[1], c = e.points[2], d = new x(r, l), f = new x(r, c);
-      this.normal = d.cross(f), this.point = r;
+    if (((i = e.points) == null ? void 0 : i.length) === 3 && e.points.every((r) => r instanceof E)) {
+      const r = e.points[0], l = e.points[1], c = e.points[2], f = new E(r, l), d = new E(r, c);
+      this.normal = f.cross(d), this.point = r;
       return;
     }
     if (((n = e.coefficients) == null ? void 0 : n.length) === 4) {
-      const [r, l, c, d] = e.coefficients;
-      this.normal = new x(r, l, c), this.point = new B(0, 0, -d);
+      const [r, l, c, f] = e.coefficients;
+      this.normal = new E(r, l, c), this.point = new B(0, 0, -f);
       return;
     }
   }
   angle(e, t, i) {
-    if (e instanceof gi)
+    if (e instanceof mi)
       return this.normal.angle(e.normal, t, i);
     let n;
-    if (e instanceof x) {
+    if (e instanceof E) {
       if (e.dimension !== 3)
         throw new Error("Vector is not 3D");
       n = e;
@@ -3766,25 +3888,25 @@ const gi = class gi {
     return this.normal.dot(e).add(this.d).isZero();
   }
 };
-Me = new WeakMap(), Ge = new WeakMap();
-let li = gi;
-var X, K, J, it, Ce, yt, Ht, wt, Be, Xt, st;
-const Yt = class Yt {
+ke = new WeakMap(), Ue = new WeakMap();
+let ai = mi;
+var Y, _, ee, st, Ce, yt, Xt, wt, Be, Yt, nt;
+const Qt = class Qt {
   constructor(...e) {
-    p(this, X, new B());
-    p(this, K, new B());
-    p(this, J, new B());
-    p(this, it, {
+    m(this, Y, new B());
+    m(this, _, new B());
+    m(this, ee, new B());
+    m(this, st, {
       AB: new R(),
       AC: new R(),
       BC: new R()
     });
-    p(this, Ce, {
+    m(this, Ce, {
       AB: new B(),
       AC: new B(),
       BC: new B()
     });
-    p(this, yt, null);
+    m(this, yt, null);
     // ------------------------------------------
     // Creation / parsing functions
     // ------------------------------------------
@@ -3792,7 +3914,6 @@ const Yt = class Yt {
      * Parse values to a triangle. Supported formats:
      * Vector2D, Vector2D, Vector2D
      * x1, y1, x2, y2, x3, y3
-     * TODO: Something else ?
      * @param values
      */
     a(this, "parse", (...e) => {
@@ -3801,9 +3922,9 @@ const Yt = class Yt {
         if (t.some((i) => i.isNaN()))
           throw new Error("One of the values is not a valid number");
         return this.parse(
-          new x(t[0], t[1]),
-          new x(t[2], t[3]),
-          new x(t[4], t[5])
+          new E(t[0], t[1]),
+          new E(t[2], t[3]),
+          new E(t[4], t[5])
         );
       } else if (e.length === 3) {
         if (e.every((t) => typeof t == "string"))
@@ -3812,36 +3933,36 @@ const Yt = class Yt {
           );
         if (e.every((t) => t instanceof R)) {
           const t = e[0].clone(), i = e[1].clone(), n = e[2].clone();
-          h(this, it, { AB: t, BC: i, AC: n });
+          h(this, st, { AB: t, BC: i, AC: n });
           let r = t.intersection(i);
           if (r.hasIntersection)
-            h(this, K, r.point.clone());
+            h(this, _, r.point.clone());
           else
             throw new Error("Lines do not intersect !");
           if (r = i.intersection(n), r.hasIntersection)
-            h(this, J, r.point.clone());
+            h(this, ee, r.point.clone());
           else
             throw new Error("Lines do not intersect !");
           if (r = n.intersection(t), r.hasIntersection)
-            h(this, X, r.point.clone());
+            h(this, Y, r.point.clone());
           else
             throw new Error("Lines do not intersect !");
-        } else e.every((t) => t instanceof B) && (h(this, X, e[0].clone()), h(this, K, e[1].clone()), h(this, J, e[2].clone()), h(this, it, {
-          AB: new R(s(this, X), s(this, K)),
-          BC: new R(s(this, K), s(this, J)),
-          AC: new R(s(this, X), s(this, J))
+        } else e.every((t) => t instanceof B) && (h(this, Y, e[0].clone()), h(this, _, e[1].clone()), h(this, ee, e[2].clone()), h(this, st, {
+          AB: new R(s(this, Y), s(this, _)),
+          BC: new R(s(this, _), s(this, ee)),
+          AC: new R(s(this, Y), s(this, ee))
         }));
-      } else if (e.length === 1 && e[0] instanceof Yt)
+      } else if (e.length === 1 && e[0] instanceof Qt)
         return e[0].clone();
-      return s(this, Ht).call(this), this;
+      return s(this, Xt).call(this), this;
     });
     /**
      * Clone the Triangle class
      */
-    a(this, "clone", () => new Yt(
-      s(this, X).clone(),
-      s(this, K).clone(),
-      s(this, J).clone()
+    a(this, "clone", () => new Qt(
+      s(this, Y).clone(),
+      s(this, _).clone(),
+      s(this, ee).clone()
     ));
     // ------------------------------------------
     // Triangle operations and properties
@@ -3849,79 +3970,79 @@ const Yt = class Yt {
     /**
      * Generate the Line object for the three segments of the triangle
      */
-    p(this, Ht, () => {
-      s(this, X).asPoint = !0, s(this, K).asPoint = !0, s(this, J).asPoint = !0, h(this, Ce, {
-        AB: new B().middleOf(s(this, X), s(this, K)),
-        AC: new B().middleOf(s(this, X), s(this, J)),
-        BC: new B().middleOf(s(this, K), s(this, J))
-      }), h(this, yt, s(this, Xt).call(this));
+    m(this, Xt, () => {
+      s(this, Y).asPoint = !0, s(this, _).asPoint = !0, s(this, ee).asPoint = !0, h(this, Ce, {
+        AB: new B().middleOf(s(this, Y), s(this, _)),
+        AC: new B().middleOf(s(this, Y), s(this, ee)),
+        BC: new B().middleOf(s(this, _), s(this, ee))
+      }), h(this, yt, s(this, Yt).call(this));
     });
     /**
      * Get the Vector2D class for the given name
      * @param ptName
      */
-    p(this, wt, (e) => {
+    m(this, wt, (e) => {
       switch (e.toUpperCase()) {
         case "A":
-          return s(this, X);
+          return s(this, Y);
         case "B":
-          return s(this, K);
+          return s(this, _);
         case "C":
-          return s(this, J);
+          return s(this, ee);
       }
-      return s(this, X);
+      return s(this, Y);
     });
     /**
      * Get the vector for the segment given by name.
      * @param ptName1
      * @param ptName2
      */
-    p(this, Be, (e, t) => new x(
+    m(this, Be, (e, t) => new E(
       s(this, wt).call(this, e),
       s(this, wt).call(this, t)
     ));
-    p(this, Xt, () => {
+    m(this, Yt, () => {
       const e = {
-        A: new R().fromPoints(s(this, X), s(this, Ce).BC),
-        B: new R().fromPoints(s(this, K), s(this, Ce).AC),
-        C: new R().fromPoints(s(this, J), s(this, Ce).AB),
+        A: new R().fromPoints(s(this, Y), s(this, Ce).BC),
+        B: new R().fromPoints(s(this, _), s(this, Ce).AC),
+        C: new R().fromPoints(s(this, ee), s(this, Ce).AB),
         intersection: null
       }, t = {
-        AB: new R().fromPointAndNormal(s(this, Ce).AB, new x(s(this, X), s(this, K)).normal()),
-        AC: new R().fromPointAndNormal(s(this, Ce).AC, new x(s(this, X), s(this, J)).normal()),
-        BC: new R().fromPointAndNormal(s(this, Ce).BC, new x(s(this, K), s(this, J)).normal()),
+        AB: new R().fromPointAndNormal(s(this, Ce).AB, new E(s(this, Y), s(this, _)).normal()),
+        AC: new R().fromPointAndNormal(s(this, Ce).AC, new E(s(this, Y), s(this, ee)).normal()),
+        BC: new R().fromPointAndNormal(s(this, Ce).BC, new E(s(this, _), s(this, ee)).normal()),
         intersection: null
       }, i = {
-        A: new R().fromPointAndNormal(s(this, X), new x(s(this, K), s(this, J)).normal()),
-        B: new R().fromPointAndNormal(s(this, K), new x(s(this, X), s(this, J)).normal()),
-        C: new R().fromPointAndNormal(s(this, J), new x(s(this, X), s(this, K)).normal()),
+        A: new R().fromPointAndNormal(s(this, Y), new E(s(this, _), s(this, ee)).normal()),
+        B: new R().fromPointAndNormal(s(this, _), new E(s(this, Y), s(this, ee)).normal()),
+        C: new R().fromPointAndNormal(s(this, ee), new E(s(this, Y), s(this, _)).normal()),
         intersection: null
-      }, n = s(this, st).call(this, "A"), r = s(this, st).call(this, "B"), l = s(this, st).call(this, "C"), c = {
+      }, n = s(this, nt).call(this, "A"), r = s(this, nt).call(this, "B"), l = s(this, nt).call(this, "C"), c = {
         A: n.internal,
         B: r.internal,
         C: r.internal,
         intersection: null
-      }, d = {
+      }, f = {
         A: n.external,
         B: r.external,
         C: l.external,
         intersection: null
-      }, f = {
+      }, d = {
         medians: e,
         mediators: t,
         heights: i,
         bisectors: c,
-        externalBisectors: d
+        externalBisectors: f
       };
-      return f.medians.intersection = f.medians.A.intersection(f.medians.B).point, f.mediators.intersection = f.mediators.AB.intersection(f.mediators.BC).point, f.heights.intersection = f.heights.A.intersection(f.heights.B).point, f.bisectors.intersection = f.bisectors.A.intersection(f.bisectors.B).point, f;
+      return d.medians.intersection = d.medians.A.intersection(d.medians.B).point, d.mediators.intersection = d.mediators.AB.intersection(d.mediators.BC).point, d.heights.intersection = d.heights.A.intersection(d.heights.B).point, d.bisectors.intersection = d.bisectors.A.intersection(d.bisectors.B).point, d;
     });
-    p(this, st, (e) => {
+    m(this, nt, (e) => {
       const t = this.lines;
       let i, n;
       if (e === "A" ? (i = t.AB, n = t.AC) : e === "B" ? (i = t.AB, n = t.BC) : e === "C" && (i = t.BC, n = t.AC), i === void 0 || n === void 0)
         throw new Error(`The point ${e} does not exist`);
-      const r = i.n.simplify().norm, l = n.n.simplify().norm, c = i.getEquation().multiply(l), d = n.getEquation().multiply(r), f = new R(c.clone().subtract(d).simplify()), g = new R(d.clone().subtract(c).simplify());
-      return e === "A" ? f.hitSegment(this.B, this.C) ? { internal: f, external: g } : { internal: g, external: f } : e === "B" ? f.hitSegment(this.A, this.C) ? { internal: f, external: g } : { internal: g, external: f } : e === "C" ? f.hitSegment(this.B, this.A) ? { internal: f, external: g } : { internal: g, external: f } : { internal: f, external: g };
+      const r = i.n.simplify().norm, l = n.n.simplify().norm, c = i.getEquation().multiply(l), f = n.getEquation().multiply(r), d = new R(c.clone().subtract(f).simplify()), p = new R(f.clone().subtract(c).simplify());
+      return e === "A" ? d.hitSegment(this.B, this.C) ? { internal: d, external: p } : { internal: p, external: d } : e === "B" ? d.hitSegment(this.A, this.C) ? { internal: d, external: p } : { internal: p, external: d } : e === "C" ? d.hitSegment(this.B, this.A) ? { internal: d, external: p } : { internal: p, external: d } : { internal: d, external: p };
     });
     return e.length > 0 && this.parse(...e), this;
   }
@@ -3929,13 +4050,13 @@ const Yt = class Yt {
   // Getter and setters
   // ------------------------------------------
   get A() {
-    return s(this, X);
+    return s(this, Y);
   }
   get B() {
-    return s(this, K);
+    return s(this, _);
   }
   get C() {
-    return s(this, J);
+    return s(this, ee);
   }
   get AB() {
     return s(this, Be).call(this, "A", "B");
@@ -3965,15 +4086,15 @@ const Yt = class Yt {
     return this.AB.normSquare.isEqual(this.BC.normSquare) || this.AB.normSquare.isEqual(this.AC.normSquare) || this.BC.normSquare.isEqual(this.AC.normSquare);
   }
   get lines() {
-    return s(this, it);
+    return s(this, st);
   }
   get remarquables() {
     return s(this, yt);
   }
 };
-X = new WeakMap(), K = new WeakMap(), J = new WeakMap(), it = new WeakMap(), Ce = new WeakMap(), yt = new WeakMap(), Ht = new WeakMap(), wt = new WeakMap(), Be = new WeakMap(), Xt = new WeakMap(), st = new WeakMap();
-let ci = Yt;
-function Ot(o) {
+Y = new WeakMap(), _ = new WeakMap(), ee = new WeakMap(), st = new WeakMap(), Ce = new WeakMap(), yt = new WeakMap(), Xt = new WeakMap(), wt = new WeakMap(), Be = new WeakMap(), Yt = new WeakMap(), nt = new WeakMap();
+let li = Qt;
+function qt(o) {
   const e = Object.assign(
     {
       negative: !0,
@@ -3984,16 +4105,16 @@ function Ot(o) {
     },
     o
   ), t = new u();
-  if (e.negative ? t.numerator = Y(e.max, e.zero) : t.numerator = de(e.zero ? 0 : 1, e.max), e.natural)
+  if (e.negative ? t.numerator = Q(e.max, e.zero) : t.numerator = pe(e.zero ? 0 : 1, e.max), e.natural)
     t.denominator = 1;
   else {
     let i = 0;
     for (; t.isRelative() && i < 10; )
-      t.denominator = de(1, e.max), i++;
+      t.denominator = pe(1, e.max), i++;
   }
   return e.reduced ? t.reduce() : t;
 }
-function Pi(o) {
+function Ri(o) {
   const e = Object.assign(
     {
       letters: "x",
@@ -4002,8 +4123,8 @@ function Pi(o) {
       zero: !1
     },
     o
-  ), t = new M();
-  if (t.coefficient = Ot({
+  ), t = new k();
+  if (t.coefficient = qt({
     zero: e.zero,
     reduced: !0,
     natural: !e.fraction
@@ -4011,14 +4132,14 @@ function Pi(o) {
     for (const i of e.letters.split(""))
       t.setLetter(i, 0);
     for (let i = 0; i < e.degree; i++) {
-      const n = pi(e.letters.split(""));
+      const n = di(e.letters.split(""));
       t.setLetter(n, t.degree(n).clone().add(1));
     }
   } else
     t.setLetter(e.letters, e.degree);
   return t;
 }
-const as = {
+const us = {
   letters: "x",
   degree: 2,
   fraction: !1,
@@ -4029,14 +4150,14 @@ const as = {
   numberOfMonoms: 0,
   positive: !0
 };
-function Bi(o) {
+function zi(o) {
   const e = Object.assign(
-    as,
+    us,
     o
   ), t = new O().empty();
   let i;
   for (let n = e.degree; n >= 0; n--)
-    i = Pi({
+    i = Ri({
       letters: e.letters,
       degree: n,
       fraction: e.fraction,
@@ -4044,12 +4165,12 @@ function Bi(o) {
     }), e.unit && e.degree === n && i.coefficient.one(), t.add(i);
   if (e.positive && t.monomByDegree().coefficient.isNegative() && t.monomByDegree().coefficient.opposite(), e.numberOfMonoms && e.numberOfMonoms > 0 && e.numberOfMonoms < t.length)
     for (; t.length > e.numberOfMonoms; ) {
-      const n = de(1, t.length - 1);
+      const n = pe(1, t.length - 1);
       t.monoms.splice(n, 1);
     }
   return t.reduce();
 }
-function ls(o) {
+function fs(o) {
   const e = Object.assign(
     {
       letters: "x",
@@ -4071,7 +4192,7 @@ function ls(o) {
     o
   ), t = new O().one();
   for (let i = 0; i < e.degree; i++) {
-    const n = Bi({
+    const n = zi({
       degree: 1,
       unit: e.unit,
       fraction: e.fraction,
@@ -4082,7 +4203,7 @@ function ls(o) {
   }
   return new H(t, 0);
 }
-function ui(o) {
+function ci(o) {
   const e = Object.assign(
     {
       axis: !0,
@@ -4091,10 +4212,10 @@ function ui(o) {
       quadrant: null
     },
     o
-  ), t = e.axis === "x", i = e.axis === "y", n = e.fraction ? Ot({ max: e.max, zero: t }) : new u(Y(e.max, t)), r = e.fraction ? Ot({ max: e.max, zero: i }) : new u(Y(e.max, i));
+  ), t = e.axis === "x", i = e.axis === "y", n = e.fraction ? qt({ max: e.max, zero: t }) : new u(Q(e.max, t)), r = e.fraction ? qt({ max: e.max, zero: i }) : new u(Q(e.max, i));
   return Number(e.quadrant) === 1 && (n.abs(), r.abs()), Number(e.quadrant) === 2 && (n.isPositive() && n.opposite(), r.isNegative() && r.opposite()), Number(e.quadrant) === 3 && (n.isPositive() && n.opposite(), r.isPositive() && r.opposite()), Number(e.quadrant) === 4 && (n.isNegative() && n.opposite(), r.isPositive() && r.opposite()), new B(n, r);
 }
-function cs(o) {
+function ds(o) {
   const e = Object.assign(
     {
       center: {
@@ -4104,116 +4225,116 @@ function cs(o) {
       pointsOnCircle: 8
     },
     o
-  ), t = ui(e.center);
+  ), t = ci(e.center);
   let i, n;
-  return e.pointsOnCircle === 8 ? (i = de(1, 3), n = i ** 2 + (i + 1) ** 2) : n = de(1, 20), new Tt(t, n, !0);
+  return e.pointsOnCircle === 8 ? (i = pe(1, 3), n = i ** 2 + (i + 1) ** 2) : n = pe(1, 20), new Tt(t, n, !0);
 }
-function us(o) {
+function ps(o) {
   const e = Object.assign(
     {
       A: {
-        x: Y(10),
-        y: Y(10)
+        x: Q(10),
+        y: Q(10)
       }
     },
     o
-  ), t = new x(
-    Y(10),
-    Y(10)
+  ), t = new E(
+    Q(10),
+    Q(10)
   );
   for (; t.isNull; )
-    t.x = Y(10), t.y = Y(10);
-  return e.slope === 1 ? t.x.sign() !== t.y.sign() && t.y.opposite() : e.slope === -1 && t.x.sign() !== t.y.sign() && t.y.opposite(), new R(new x(e.A.x, e.A.y), t);
+    t.x = Q(10), t.y = Q(10);
+  return e.slope === 1 ? t.x.sign() !== t.y.sign() && t.y.opposite() : e.slope === -1 && t.x.sign() !== t.y.sign() && t.y.opposite(), new R(new E(e.A.x, e.A.y), t);
 }
-function fs(o) {
+function ms(o) {
   const e = Object.assign(
     {
       A: {
-        x: Y(10),
-        y: Y(10),
-        z: Y(10)
+        x: Q(10),
+        y: Q(10),
+        z: Q(10)
       },
       direction: {
-        x: Y(10),
-        y: Y(10),
-        z: Y(10)
+        x: Q(10),
+        y: Q(10),
+        z: Q(10)
       }
     },
     o
-  ), t = new B(e.A.x, e.A.y, e.A.z), i = new x(e.direction.x, e.direction.y, e.direction.z);
-  return new At(t, i);
+  ), t = new B(e.A.x, e.A.y, e.A.z), i = new E(e.direction.x, e.direction.y, e.direction.z);
+  return new Ot(t, i);
 }
-const ds = {
-  equation: (o) => ls(o),
-  polynom: (o) => Bi(o),
-  monom: (o) => Pi(o),
-  fraction: (o) => Ot(o),
-  number: (o, e, t) => de(o, e, t),
-  numberSym: (o, e) => Y(o, e),
-  prime: (o) => rs(o),
-  bool: (o) => qi(o),
-  array: (o, e) => os(o, e),
-  item: (o) => pi(o),
-  shuffle: (o) => ki(o),
-  line: (o) => us(o),
-  line3: (o) => fs(o),
-  vector: (o) => ui(o),
+const gs = {
+  equation: (o) => fs(o),
+  polynom: (o) => zi(o),
+  monom: (o) => Ri(o),
+  fraction: (o) => qt(o),
+  number: (o, e, t) => pe(o, e, t),
+  numberSym: (o, e) => Q(o, e),
+  prime: (o) => as(o),
+  bool: (o) => Ci(o),
+  array: (o, e) => ls(o, e),
+  item: (o) => di(o),
+  shuffle: (o) => $i(o),
+  line: (o) => ps(o),
+  line3: (o) => ms(o),
+  vector: (o) => ci(o),
   point: (o) => {
-    const e = ui(o);
+    const e = ci(o);
     return e.asPoint = !0, e;
   },
-  circle: (o) => cs(o)
-}, ms = {
-  Numeric: G,
+  circle: (o) => ds(o)
+}, ws = {
+  Numeric: U,
   Fraction: u,
   Root: vt,
-  Monom: M,
+  Monom: k,
   Polynom: O,
   Equation: H,
-  Matrix: hs,
-  LinearSystem: ii,
-  Factor: le,
-  PolyFactor: ni,
-  LogicalSet: es,
-  Random: ds,
+  Matrix: cs,
+  LinearSystem: ti,
+  Factor: ae,
+  PolyFactor: si,
+  LogicalSet: ss,
+  Random: gs,
   Geometry: {
-    Vector: x,
+    Vector: E,
     Point: B,
     Line: R,
-    Triangle: ci,
+    Triangle: li,
     Circle: Tt,
-    Line3: At,
-    Plane3: li
+    Line3: Ot,
+    Plane3: ai
   },
-  NumExp: _i
+  NumExp: is
 };
 export {
   Tt as Circle,
   H as Equation,
-  Nt as EquationSolver,
-  at as FACTOR_DISPLAY,
-  le as Factor,
+  At as EquationSolver,
+  lt as FACTOR_DISPLAY,
+  ae as Factor,
   u as Fraction,
   R as Line,
-  At as Line3,
-  ii as LinearSystem,
-  es as LogicalSet,
-  hs as Matrix,
-  M as Monom,
+  Ot as Line3,
+  ti as LinearSystem,
+  ss as LogicalSet,
+  cs as Matrix,
+  k as Monom,
   vt as NthRoot,
-  _i as NumExp,
-  G as Numeric,
-  li as Plane3,
+  is as NumExp,
+  U as Numeric,
+  ai as Plane3,
   B as Point,
-  ni as PolyFactor,
+  si as PolyFactor,
   O as Polynom,
-  ds as Random,
-  ci as Triangle,
-  x as Vector,
-  is as areVectorsColinears,
-  ts as areVectorsEquals,
-  ms as default,
-  ns as determinant,
-  ss as dotProduct
+  gs as Random,
+  li as Triangle,
+  E as Vector,
+  rs as areVectorsColinears,
+  ns as areVectorsEquals,
+  ws as default,
+  hs as determinant,
+  os as dotProduct
 };
 //# sourceMappingURL=pimath.js.map