pimath 0.1.20 → 0.1.22

package/dist/pimath.js

@@ -4,7 +4,7 @@ var yi = (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), O = (o, e, t) => (Jt(o, e, "access private method"), 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 = [];
   let i, n;
@@ -391,37 +391,37 @@ class vt {
   }
 }
 Z = new WeakMap(), ee = new WeakMap(), ie = new WeakMap(), Xe = new WeakMap();
-var ce, je, P, ot, ye, bi, ti, xi, Ei, Ni;
+var fe, je, I, ot, ye, bi, ti, xi, Ei, Ni;
 const mi = class mi {
   constructor(e, t, i = "x") {
-    p(this, P);
-    p(this, ce);
+    p(this, I);
+    p(this, fe);
     p(this, je);
     if (h(this, je, i), Object.hasOwn(e, "moveLeft")) {
       const n = e;
-      h(this, ce, n.left.clone().subtract(n.right));
+      h(this, fe, n.left.clone().subtract(n.right));
     } else
-      h(this, ce, e.clone().subtract(t ?? 0));
+      h(this, fe, e.clone().subtract(t ?? 0));
   }
   solve() {
-    if (s(this, ce).degree().isOne())
-      return O(this, P, xi).call(this);
-    if (s(this, ce).degree().value === 2)
-      return O(this, P, Ei).call(this);
-    const e = O(this, P, bi).call(this);
+    if (s(this, fe).degree().isOne())
+      return q(this, I, xi).call(this);
+    if (s(this, fe).degree().value === 2)
+      return q(this, I, Ei).call(this);
+    const e = q(this, I, bi).call(this);
     if (e.length > 0)
       return e;
-    if (s(this, ce).degree().value === 3)
-      return O(this, P, ti).call(this);
+    if (s(this, fe).degree().value === 3)
+      return q(this, I, ti).call(this);
     throw new Error("The equation degree is too high.");
   }
   solveAsCardan() {
-    if (s(this, ce).degree().value !== 3)
+    if (s(this, fe).degree().value !== 3)
       throw new Error("The equation is not cubic.");
-    return O(this, P, ti).call(this);
+    return q(this, I, ti).call(this);
   }
 };
-ce = new WeakMap(), je = new WeakMap(), P = new WeakSet(), ot = function(e, t) {
+fe = new WeakMap(), je = new WeakMap(), I = new WeakSet(), ot = function(e, t) {
   return {
     variable: s(this, je),
     exact: !1,
@@ -431,7 +431,7 @@ ce = new WeakMap(), je = new WeakMap(), P = new WeakSet(), ot = function(e, t) {
   };
 }, ye = function(e) {
   if (e instanceof u && e.isApproximative())
-    return O(this, P, ot).call(this, e.value);
+    return q(this, I, ot).call(this, e.value);
   const t = new u(e);
   return {
     variable: s(this, je),
@@ -441,23 +441,23 @@ ce = new WeakMap(), je = new WeakMap(), P = new WeakSet(), ot = function(e, t) {
     display: t.display
   };
 }, bi = function() {
-  let e = s(this, ce).clone(), t = [];
+  let e = s(this, fe).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(O(this, P, ye).call(this, 0)), e = e.divide("x"), r = e.monomByDegree(0).coefficient;
+    t.length === 0 && t.push(q(this, I, 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 E = new u(g, f);
-      e.evaluate(E).isZero() && !t.find((V) => V.value === E.value) && t.push(O(this, P, ye).call(this, E)), E.opposite(), e.evaluate(E).isZero() && !t.find((V) => V.value === E.value) && t.push(O(this, P, ye).call(this, E));
+      const N = new u(g, f);
+      e.evaluate(N).isZero() && !t.find((V) => V.value === N.value) && t.push(q(this, I, ye).call(this, N)), N.opposite(), e.evaluate(N).isZero() && !t.find((V) => V.value === N.value) && t.push(q(this, I, ye).call(this, N));
     }
   for (const f of t) {
     if (f.exact !== !1 && f.exact.isZero())
       continue;
-    const g = s(this, ce).clone().parse("x", f.exact.denominator, -f.exact.numerator);
+    const g = s(this, fe).clone().parse("x", f.exact.denominator, -f.exact.numerator);
     for (; e.isDividableBy(g); )
       e = e.divide(g);
   }
@@ -468,61 +468,61 @@ ce = new WeakMap(), je = new WeakMap(), P = new WeakSet(), ot = function(e, t) {
   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, ce), 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)), E = g.clone().opposite(), V = f.clone().opposite().pow(3).divide(27), de = E.clone().pow(2).subtract(V.clone().multiply(4)).opposite();
-  if (de.isNegative()) {
-    const he = g.clone().opposite().add(de.clone().opposite().sqrt()).divide(2).root(3), ae = g.clone().opposite().subtract(de.clone().opposite().sqrt()).divide(2).root(3), ge = he.clone().add(ae).subtract(l.clone().divide(3));
-    return [O(this, P, ye).call(this, ge)];
+  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, I, ye).call(this, ge)];
   }
-  if (de.isZero()) {
+  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) ? [O(this, P, ye).call(this, ae)] : [
-      O(this, P, ye).call(this, ge),
-      O(this, P, ye).call(this, ae)
-    ].sort((pe, le) => pe.value - le.value);
+    return ae.isEqual(ge) ? [q(this, I, ye).call(this, ae)] : [
+      q(this, I, ye).call(this, ge),
+      q(this, I, ye).call(this, ae)
+    ].sort((me, ce) => me.value - ce.value);
   }
-  if (de.isPositive()) {
-    const he = [], ae = f.value, ge = g.value, pe = l.value;
-    for (let le = 0; le < 3; le++)
-      he.push(2 * Math.sqrt(-ae / 3) * Math.cos(Math.acos(3 * ge / (2 * ae) * Math.sqrt(-3 / ae)) / 3 + 2 * Math.PI * le / 3) - pe / 3);
-    return he.map((le) => O(this, P, ot).call(this, le)).sort((le, nt) => le.value - nt.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, I, ot).call(this, ce)).sort((ce, nt) => ce.value - nt.value);
   }
   return [];
 }, xi = function() {
-  const e = s(this, ce).monomByDegree(0).coefficient.clone().opposite().divide(s(this, ce).monomByDegree(1).coefficient);
+  const e = s(this, fe).monomByDegree(0).coefficient.clone().opposite().divide(s(this, fe).monomByDegree(1).coefficient);
   return [
-    O(this, P, ye).call(this, e)
+    q(this, I, ye).call(this, e)
   ];
 }, Ei = function() {
-  const e = s(this, ce), 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));
+  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));
   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() ? [O(this, P, ye).call(this, c)] : [
-      O(this, P, ye).call(this, c),
-      O(this, P, ye).call(this, d)
+    return l.isZero() ? [q(this, I, ye).call(this, c)] : [
+      q(this, I, ye).call(this, c),
+      q(this, I, ye).call(this, d)
     ].sort((f, g) => f.value - g.value);
   }
-  return O(this, P, Ni).call(this, t, i, r);
+  return q(this, I, Ni).call(this, t, i, r);
 }, Ni = function(e, t, i) {
-  const n = G.dividers(i.value).filter((pe) => Math.sqrt(pe) % 1 === 0).map((pe) => Math.sqrt(pe)).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} `, E = n === 1 ? "+" : `+${f} `;
-  function V(pe, le, nt, Kt) {
-    return `\\frac{ ${le} ${nt}\\sqrt{ ${Kt} } }{ ${pe} }`;
+  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 de(pe, le, nt, Kt) {
-    return `(${le}${nt}sqrt(${Kt}))/${pe}`;
+  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);
   return [
-    O(this, P, ot).call(this, ae, {
+    q(this, I, ot).call(this, ae, {
       tex: V(c.tex, l.tex, g.toString(), d.tex),
-      display: de(c.display, l.display, g.toString(), d.display)
+      display: pe(c.display, l.display, g.toString(), d.display)
     }),
-    O(this, P, ot).call(this, ge, {
-      tex: V(c.tex, l.tex, E.toString(), d.tex),
-      display: de(c.display, l.display, E.toString(), d.display)
+    q(this, I, ot).call(this, ge, {
+      tex: V(c.tex, l.tex, N.toString(), d.tex),
+      display: pe(c.display, l.display, N.toString(), d.display)
     })
-  ].sort((pe, le) => pe.value - le.value);
+  ].sort((me, ce) => me.value - ce.value);
 };
 let Nt = mi;
 var Gi = Object.defineProperty, Ti = (o) => {
@@ -536,22 +536,22 @@ var m = /* @__PURE__ */ ((o) => (o.VARIABLE = "variable", o.COEFFICIENT = "coeff
 function Hi(o, e) {
   if (o.length <= 1)
     return o;
-  const t = Object.keys(e).filter((E) => e[E].type === m.FUNCTION).map((E) => E);
-  t.sort((E, V) => V.length - E.length);
+  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((E, V) => V.length - E.length);
+  n.sort((N, V) => V.length - N.length);
   const r = new RegExp(`^(${n.join("|")})`), l = /^(\d+(\.\d+)?)/;
   let c = "", d, f, g;
   for (; o.length > 0; ) {
     if (d = f, g = void 0, t.length > 0 && i.exec(o)) {
-      const E = t.find((V) => o.startsWith(V));
-      E && (g = E + "(", o = o.slice(E.length + 1), f = m.FUNCTION);
+      const N = t.find((V) => o.startsWith(V));
+      N && (g = N + "(", o = o.slice(N.length + 1), f = m.FUNCTION);
     } else if (n.length > 0 && r.exec(o)) {
-      const E = n.find((V) => o.startsWith(V));
-      E && (g = E, o = o.slice(E.length), f = m.CONSTANT);
+      const N = n.find((V) => o.startsWith(V));
+      N && (g = N, o = o.slice(N.length), f = m.CONSTANT);
     } else if (l.exec(o)) {
-      const E = l.exec(o);
-      E && (g = E[0], o = o.slice(E[0].length), f = m.COEFFICIENT);
+      const N = l.exec(o);
+      N && (g = N[0], o = o.slice(N[0].length), f = m.COEFFICIENT);
     } else
       switch (g = o[0], o = o.slice(1), g) {
         case "(":
@@ -703,17 +703,17 @@ class Qt {
           break;
         case m.OPERATION:
           if (n.length > 0) {
-            let E = n[n.length - 1];
-            for (g = +d; E.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)[E.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)[E.token].precedence); ) {
+            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) {
                 console.log("SECURITY LEVEL 2 OPERATION EXIT");
                 break;
               }
               if (i.push(n.pop() ?? { token: "", tokenType: m.OPERATION }), n.length === 0)
                 break;
-              E = n[n.length - 1];
+              N = n[n.length - 1];
             }
           }
           n.push({ token: r, tokenType: c });
@@ -852,7 +852,7 @@ class _i {
   }
 }
 var T, v, Ye, xt, Le, qt, kt;
-const I = class I {
+const P = class P {
   constructor(e) {
     p(this, Ye);
     p(this, T);
@@ -861,7 +861,7 @@ const I = class I {
      * Clone the current Monom.
      */
     a(this, "clone", () => {
-      const e = new I();
+      const e = new P();
       e.coefficient = s(this, T).clone();
       for (const t in s(this, v))
         e.setLetter(t, s(this, v)[t].clone());
@@ -873,8 +873,8 @@ 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() && O(this, Ye, xt).call(this, i), s(this, T).add(i.coefficient)) : console.log("Add monom: " + this.display + " is not similar with ", i.display);
+        const i = t instanceof P ? t : new P(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);
       }
       return this;
     });
@@ -893,7 +893,7 @@ const I = class I {
         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;
       } else
-        return new I().zero();
+        return new P().zero();
     });
     /**
      * Divide the current monoms by multiple monoms
@@ -901,7 +901,7 @@ const I = class I {
      */
     a(this, "divide", (...e) => {
       for (const t of e) {
-        const i = t instanceof I ? t : new I(t);
+        const i = t instanceof P ? t : new P(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);
@@ -1008,7 +1008,7 @@ const I = class I {
      */
     a(this, "multiply", (...e) => {
       for (const t of e) {
-        const i = t instanceof I ? t : new I(t);
+        const i = t instanceof P ? t : new P(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();
@@ -1075,8 +1075,8 @@ 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() && O(this, Ye, xt).call(this, i), s(this, T).add(i.clone().coefficient.opposite())) : console.log("Subtract: Is not similar: ", i.display);
+        const i = t instanceof P ? t : new P(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);
       }
       return this;
     });
@@ -1122,23 +1122,23 @@ const I = class I {
       var d;
       let i, n, r, l, c;
       if (t.tokenType === m.COEFFICIENT)
-        e.push(new I(new u(t.token)));
+        e.push(new P(new u(t.token)));
       else if (t.tokenType === m.VARIABLE) {
-        const f = new I().one();
+        const f = new P().one();
         f.setLetter(t.token, 1), e.push(f.clone());
       } else if (t.tokenType === m.OPERATION)
         switch (t.token) {
           case "-":
-            n = e.pop() ?? new I().zero(), i = e.pop() ?? new I().zero(), e.push(i.subtract(n));
+            n = e.pop() ?? new P().zero(), i = e.pop() ?? new P().zero(), e.push(i.subtract(n));
             break;
           case "*":
-            n = e.pop() ?? new I().one(), i = e.pop() ?? new I().one(), e.push(i.multiply(n));
+            n = e.pop() ?? new P().one(), i = e.pop() ?? new P().one(), e.push(i.multiply(n));
             break;
           case "/":
-            n = e.pop() ?? new I().one(), i = e.pop() ?? new I().one(), e.push(i.divide(n));
+            n = e.pop() ?? new P().one(), i = e.pop() ?? new P().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 = ((d = e.pop()) == null ? void 0 : d.coefficient) ?? new u().one(), r = e.pop() ?? new P().one(), l = r.variables[0], l && r.setLetter(l, c), e.push(r);
             break;
           }
         }
@@ -1151,7 +1151,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()), O(this, Ye, xt).call(this, e)), this;
+    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 P && (h(this, T, s(e, T).clone()), q(this, Ye, xt).call(this, e)), this;
   }
   /**
    * Get the coefficient \\(k\\) of the Monom \\(k\\cdot x^{n}\\)
@@ -1193,20 +1193,20 @@ const I = class I {
     if (t.length > 0 && e.length > 0)
       for (const n of e)
         for (const r of t) {
-          const l = new I();
+          const l = new P();
           l.coefficient = new u(n), l.literal = r, i.push(l);
         }
     else if (e.length === 0)
       for (const n of t) {
-        const r = new I();
+        const r = new P();
         r.coefficient = new u().one(), r.literal = n, i.push(r);
       }
     else
       for (const n of e) {
-        const r = new I();
+        const r = new P();
         r.coefficient = new u(n), i.push(r);
       }
-    return i.length === 0 ? [new I().one()] : i;
+    return i.length === 0 ? [new P().one()] : i;
   }
   integrate(e, t, i) {
     const n = this.primitive(i);
@@ -1301,11 +1301,11 @@ const I = class I {
 T = new WeakMap(), v = new WeakMap(), Ye = 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) => {
+}, Le = new WeakMap(), qt = new WeakMap(), kt = new WeakMap(), a(P, "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));
+      return new P().zero();
+  const t = new P(), i = G.gcd(...e.map((r) => r.coefficient.numerator)), n = G.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)
@@ -1318,20 +1318,20 @@ T = new WeakMap(), v = new WeakMap(), Ye = new WeakSet(), xt = function(e) {
  * Multiply two monoms and return a NEW monom.
  * @param monoms
  */
-a(I, "xMultiply", (...e) => {
-  const t = new I().one();
+a(P, "xMultiply", (...e) => {
+  const t = new P().one();
   for (const i of e)
     t.multiply(i);
   return t;
 });
-let M = I;
+let M = P;
 function wi(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);
 }
-var Ue, y, Qe, ut, Ke, ft, Mt, Ct, $t, Je, It, dt, Pt, Bt, St, Rt, Oi, zt, Lt;
+var Ue, y, Qe, ut, Ke, ft, Mt, Ct, $t, Je, Pt, dt, It, Bt, St, Rt, Oi, zt, Lt;
 const C = class C {
   constructor(e, ...t) {
     p(this, Rt);
@@ -1346,7 +1346,7 @@ const C = class C {
      */
     a(this, "parse", (e, ...t) => {
       if (h(this, y, []), h(this, Ue, []), typeof e == "string")
-        return O(this, Rt, Oi).call(this, e, ...t);
+        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)
@@ -1462,7 +1462,7 @@ const C = class C {
           t.push(i.clone()), i.one();
           break;
         } else {
-          let c = s(this, It).call(this, i, l, e ?? "x");
+          let c = s(this, Pt).call(this, i, l, e ?? "x");
           for (l = i.degree(e).value; c.length > 0; ) {
             const d = c[0];
             if (!i.isDividableBy(d))
@@ -1470,8 +1470,8 @@ const C = class C {
             else {
               const f = i.euclidean(d);
               t.push(d), i = f.quotient.clone(), c = c.filter((g) => {
-                const E = i.monoms[0], V = i.monoms[i.monoms.length - 1], de = g.monoms[0], he = g.monoms[g.monoms.length - 1];
-                return V.isDivisible(he) ? E.isDivisible(de) : !1;
+                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;
               });
             }
           }
@@ -1574,7 +1574,7 @@ 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, "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, It).call(this, e) : this);
     a(this, "one", () => (h(this, y, []), s(this, y).push(new M().one()), this));
     // ------------------------------------------
     a(this, "opposite", () => (h(this, y, s(this, y).map((e) => e.opposite())), this));
@@ -1677,7 +1677,7 @@ const C = class C {
       }), t;
     });
     p(this, $t, (e) => {
-      var E;
+      var N;
       let t, i, n, r, l, c, d, f, g;
       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() ? [
@@ -1689,10 +1689,10 @@ const C = class C {
           new C(e).subtract(f.display).multiply(f.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), de = s(E = V, $t).call(E, "x"), he = [];
+        const V = new C("x", n.coefficient, r.coefficient, l.coefficient), pe = s(N = V, $t).call(N, "x"), he = [];
         let ae;
-        if (de.length >= 2) {
-          for (const ge of de)
+        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;
         }
@@ -1709,7 +1709,7 @@ 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) => {
+    p(this, Pt, (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) => {
@@ -1722,7 +1722,7 @@ const C = class C {
         t.coefficient.multiply(e);
       return this.reduce();
     });
-    p(this, Pt, (e) => s(this, dt).call(this, new u(e)));
+    p(this, It, (e) => s(this, dt).call(this, new u(e)));
     p(this, Bt, (e) => {
       for (const t of s(this, y))
         t.multiply(e);
@@ -1865,8 +1865,8 @@ const C = class C {
           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 E = f === i.length - 1 ? d.value + 1 : (d.value + i[f + 1].value) / 2;
-        c = this.evaluate(E, !0) < 0 ? "-" : "+", n[r] = "z", n = Re(n, "", c, r, l), r = +l, l = n.length;
+        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;
       });
     }
     return { roots: t, signs: n };
@@ -1884,7 +1884,7 @@ 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) {
+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(), Pt = new WeakMap(), dt = new WeakMap(), It = new WeakMap(), Bt = new WeakMap(), St = new WeakMap(), Rt = new WeakSet(), Oi = function(e, ...t) {
   if (t.length === 0) {
     if (e = "" + e, e !== "" && !isNaN(Number(e))) {
       this.empty();
@@ -1913,7 +1913,7 @@ Ue = new WeakMap(), y = new WeakMap(), Qe = new WeakMap(), ut = new WeakMap(), K
   } else
     return this.zero();
 }, zt = new WeakMap(), Lt = new WeakMap();
-let q = C;
+let O = C;
 var A, k, se, Dt, _e, Vt;
 const Ee = class Ee {
   constructor(e, t, i) {
@@ -1929,7 +1929,7 @@ const Ee = class Ee {
       if (t === !1)
         throw new Error("The equation is not valid (no sign found)");
       const i = e.split(t);
-      return this.create(new q(i[0]), new q(i[1]), s(this, _e).call(this, t));
+      return this.create(new O(i[0]), new O(i[1]), s(this, _e).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)));
@@ -2048,11 +2048,11 @@ 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 q().zero()), h(this, k, new q().zero()), h(this, se, "="), e !== void 0 && t === void 0) {
+    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)
         return e.clone();
       typeof e == "string" && this.parse(e);
-    } else e !== void 0 && t !== void 0 && (this.left = new q(e), this.right = new q(t));
+    } else e !== void 0 && t !== void 0 && (this.left = new O(e), this.right = new O(t));
     return i !== void 0 && (this.sign = i), this;
   }
   /**
@@ -2068,7 +2068,7 @@ const Ee = class 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));
-    const t = new q(e);
+    const t = new O(e);
     return s(this, A).add(t), s(this, k).add(t), this;
   }
   /**
@@ -2099,7 +2099,7 @@ const Ee = class 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));
-    const t = new q(e);
+    const t = new O(e);
     return s(this, A).subtract(t), s(this, k).subtract(t), this;
   }
   static isEquationString(e) {
@@ -2158,9 +2158,9 @@ const ve = class 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 q(i)), h(this, Ae, new u(n.replace("(", "").replace(")", "")));
+      h(this, Te, new O(i)), h(this, Ae, new u(n.replace("(", "").replace(")", "")));
     } else
-      h(this, Te, new q(e)), h(this, Ae, new u(t ?? 1));
+      h(this, Te, new O(e)), h(this, Ae, new u(t ?? 1));
     return h(this, Ne, 1), this;
   }
   parse() {
@@ -2203,7 +2203,7 @@ const ve = class ve {
   divide(e) {
     if (e instanceof ve && this.isSameAs(e))
       return this.power.subtract(e.power), this;
-    const t = new q(e);
+    const t = new O(e);
     if (this.isSameAs(t))
       return this.power.subtract(1), this;
     throw new Error("The two factors must be the same");
@@ -2225,7 +2225,7 @@ const ve = class ve {
   }
   isSameAs(e) {
     let t;
-    return e instanceof ve ? t = e.polynom : e instanceof q ? t = e : t = new q(e), this.polynom.isEqual(t);
+    return e instanceof ve ? t = e.polynom : e instanceof O ? t = e : t = new O(e), this.polynom.isEqual(t);
   }
   isZero() {
     return this.polynom.isZero();
@@ -2233,7 +2233,7 @@ const ve = class ve {
   multiply(e) {
     if (e instanceof ve && this.isSameAs(e))
       return this.power.add(e.power), this;
-    const t = new q(e);
+    const t = new O(e);
     if (this.isSameAs(t))
       return this.power.add(1), this;
     throw new Error("The two factors must be the same");
@@ -2292,7 +2292,7 @@ const ve = class ve {
   }
 };
 Ne = new WeakMap(), Te = new WeakMap(), Ae = new WeakMap(), et = new WeakMap();
-let ue = ve;
+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 {
   constructor(...e) {
@@ -2371,7 +2371,7 @@ const He = class He {
     const n = t.split("").splice(0, i - 1);
     return new He(
       ...e.map((r) => {
-        const l = new q(n.join(""), ...r);
+        const l = new O(n.join(""), ...r);
         return new H(l, 0);
       })
     );
@@ -2552,13 +2552,13 @@ class es {
     return s(this, Oe).filter((e) => e.tokenType === "variable").map((e) => e.token);
   }
   vennAB() {
-    return O(this, pt, si).call(this, {
+    return q(this, pt, si).call(this, {
       A: ["A", "AB"],
       B: ["B", "AB"]
     }, ["A", "B", "AB", "E"]);
   }
   vennABC() {
-    return O(this, pt, si).call(this, {
+    return q(this, pt, si).call(this, {
       A: ["A", "AB", "AC", "ABC"],
       B: ["B", "AB", "BC", "ABC"],
       C: ["C", "AC", "BC", "ABC"]
@@ -2609,34 +2609,34 @@ Oe = new WeakMap(), pt = new WeakSet(), si = function(e, t) {
       }
   return [...i[0]].sort();
 };
-var tt, N, mt, ri, gt, oi;
-const me = class me {
+var tt, E, mt, ri, gt, oi;
+const ue = class ue {
   constructor(...e) {
     p(this, gt);
     p(this, tt, at.POWER);
-    p(this, N, []);
+    p(this, E, []);
     return this.parse(...e), this;
   }
   parse(...e) {
-    return e.length === 0 ? this : (h(this, N, []), e.forEach((t) => {
+    return e.length === 0 ? this : (h(this, E, []), e.forEach((t) => {
       if (typeof t == "string") {
         const i = t.split(")(").join(")*(").split("*");
-        s(this, N).push(...i.map((n) => new ue(n)));
-      } else t instanceof me ? s(this, N).push(...t.factors.map((i) => i.clone())) : s(this, N).push(new ue(t));
+        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));
     }), this);
   }
   clone() {
-    return new me(...s(this, N).map((e) => e.clone()));
+    return new ue(...s(this, E).map((e) => e.clone()));
   }
   get tex() {
-    const { num: e, den: t } = O(this, gt, oi).call(this);
+    const { num: e, den: t } = q(this, gt, oi).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 } = O(this, gt, oi).call(this);
+    const { num: e, den: t } = q(this, gt, oi).call(this);
     if (t.length === 0)
       return e.length === 1 ? e[0].asSingle.display : e.map((r) => 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("");
@@ -2645,25 +2645,25 @@ const me = class me {
   static gcd(...e) {
     var i;
     if (e.length === 0)
-      return new me().one();
+      return new ue().one();
     if (e.length === 1)
       return e[0];
     if (e.length === 2)
-      return O(i = me, mt, ri).call(i, e[0], e[1]);
+      return q(i = ue, mt, ri).call(i, e[0], e[1]);
     let t = e[0];
     return e.shift(), e.forEach((n) => {
       var r;
-      return t = O(r = me, mt, ri).call(r, t, n);
+      return t = q(r = ue, mt, ri).call(r, t, n);
     }), t;
   }
   add(...e) {
     let t = [this, ...e];
-    const i = me.gcd(...t);
+    const i = ue.gcd(...t);
     t = t.map((r) => r.divide(i).reduce());
-    const n = new q("0");
-    return t.forEach((r) => n.add(r.develop())), h(this, N, [
+    const n = new O("0");
+    return t.forEach((r) => n.add(r.develop())), h(this, E, [
       ...i.factors,
-      new ue(n)
+      new le(n)
     ]), this;
   }
   get asPower() {
@@ -2673,47 +2673,47 @@ const me = class me {
     return h(this, tt, at.ROOT), this;
   }
   degree(e) {
-    return s(this, N).reduce((t, i) => t.add(i.degree(e)), new u("0"));
+    return s(this, E).reduce((t, i) => t.add(i.degree(e)), new u("0"));
   }
   get denominator() {
-    return s(this, N).filter((e) => e.power.isNegative());
+    return s(this, E).filter((e) => e.power.isNegative());
   }
   derivative() {
-    const e = [], t = s(this, N).length;
+    const e = [], t = s(this, E).length;
     for (let n = 0; n < t; n++) {
-      const r = s(this, N).slice(), l = r.splice(n, 1)[0];
-      e.push(new me(...r).multiply(new me(...l.derivative())));
+      const r = s(this, E).slice(), l = r.splice(n, 1)[0];
+      e.push(new ue(...r).multiply(new ue(...l.derivative())));
     }
     e.forEach((n) => n.reduce());
     const i = e.shift();
-    return i !== void 0 && h(this, N, i.factors), this.add(...e);
+    return i !== void 0 && h(this, E, i.factors), this.add(...e);
   }
   develop() {
-    const e = new q("1");
-    return s(this, N).forEach((t) => {
+    const e = new O("1");
+    return s(this, E).forEach((t) => {
       e.multiply(t.develop());
     }), e;
   }
   divide(e) {
-    return h(this, N, s(this, N).concat(e.clone().factors.map((t) => t.inverse()))), this;
+    return h(this, E, s(this, E).concat(e.clone().factors.map((t) => t.inverse()))), this;
   }
   evaluate(e, t) {
-    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"));
+    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"));
   }
   get factors() {
-    return s(this, N);
+    return s(this, E);
   }
   set factors(e) {
-    h(this, N, e);
+    h(this, E, e);
   }
-  fromPolynom(e, t) {
-    return h(this, N, new q(e).factorize(t).map((i) => new ue(i))), this;
+  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, N);
+    return s(this, E);
   }
   getZeroes() {
-    const e = [].concat(...s(this, N).map((t) => t.polynom.getZeroes()));
+    const e = [].concat(...s(this, E).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,63 +2721,69 @@ const me = class me {
     );
   }
   hasVariable(e) {
-    return s(this, N).some((t) => t.hasVariable(e));
+    return s(this, E).some((t) => t.hasVariable(e));
   }
   inverse() {
-    return h(this, N, s(this, N).map((e) => e.inverse())), this;
+    return h(this, E, s(this, E).map((e) => e.inverse())), this;
   }
   isEqual(e) {
-    const t = me.gcd(this, e), i = this.clone().divide(t).reduce(), n = e.clone().divide(t).reduce();
+    const t = ue.gcd(this, e), i = this.clone().divide(t).reduce(), n = e.clone().divide(t).reduce();
     return i.isOne() && n.isOne();
   }
   isOne() {
-    return s(this, N).every((e) => e.isOne());
+    return s(this, E).every((e) => e.isOne());
   }
   isZero() {
-    return s(this, N).every((e) => e.isZero());
+    return s(this, E).every((e) => e.isZero());
   }
   multiply(...e) {
     return e.forEach((t) => {
-      h(this, N, s(this, N).concat(t.clone().factors));
+      h(this, E, s(this, E).concat(t.clone().factors));
     }), this;
   }
   get numerator() {
-    return s(this, N).filter((e) => e.power.isPositive());
+    return s(this, E).filter((e) => e.power.isPositive());
   }
   one() {
-    return h(this, N, [new ue("1", "1")]), this;
+    return h(this, E, [new le("1", "1")]), this;
   }
   opposite() {
-    const e = s(this, N).findIndex((t) => t.display === "(-1)");
-    return e >= 0 ? s(this, N).splice(e, 1) : s(this, N).push(new ue("-1", "1")), this;
+    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;
   }
   pow(e) {
-    return h(this, N, s(this, N).map((t) => t.pow(e))), this;
+    return h(this, E, s(this, E).map((t) => t.pow(e))), this;
   }
   primitive() {
     throw new Error("Method not implemented.");
   }
   reduce() {
     const e = ei(this);
-    return h(this, N, Object.values(e).map((t) => {
+    return h(this, E, Object.values(e).map((t) => {
       const i = t[0].polynom, n = t.reduce((r, l) => r.add(l.power), new u("0"));
-      return new ue(i, n.reduce());
+      return new le(i, n.reduce());
     }).filter((t) => !t.power.isZero())), this;
   }
   root(e) {
-    return h(this, N, s(this, N).map((t) => t.root(e))), this;
+    return h(this, E, s(this, E).map((t) => t.root(e))), this;
   }
   sort() {
-    return h(this, N, s(this, N).sort((e, t) => e.degree().isLeq(t.degree()) ? -1 : 1)), this;
+    return h(this, E, s(this, E).sort((e, t) => e.degree().isLeq(t.degree()) ? -1 : 1)), this;
+  }
+  splitFactors() {
+    return {
+      numerator: new ue(...s(this, E).filter((e) => e.power.isPositive())),
+      denominator: new ue(...s(this, E).filter((e) => e.power.isStrictlyNegative()))
+    };
   }
   sqrt() {
-    return h(this, N, s(this, N).map((e) => e.sqrt())), this;
+    return h(this, E, s(this, E).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 ue(n), ...n.tableOfSigns(e) }));
+    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) {
         case "d":
@@ -2796,23 +2802,23 @@ const me = class me {
     }), n), []), roots: e, factors: t };
   }
   get variables() {
-    return s(this, N).reduce((e, t) => e.concat(t.variables), []);
+    return s(this, E).reduce((e, t) => e.concat(t.variables), []);
   }
   zero() {
-    return h(this, N, [new ue("0", "1")]), this;
+    return h(this, E, [new le("0", "1")]), this;
   }
 };
-tt = new WeakMap(), N = new WeakMap(), mt = new WeakSet(), ri = function(e, t) {
+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, E) => g.add(E.power), new u("0")), f = n[c].reduce((g, E) => g.add(E.power), new u("0"));
-    return new ue(c, u.min(d, f));
+    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));
   });
-  return new me(...l);
+  return new ue(...l);
 }, gt = new WeakSet(), oi = function() {
   let e, t = [];
-  return s(this, tt) === at.ROOT ? (e = this.numerator, t = this.denominator.map((i) => i.clone().inverse())) : e = s(this, N), e.length === 0 && (e = [new ue("1")]), { num: e, den: t };
-}, p(me, mt);
-let ni = me;
+  return s(this, tt) === at.ROOT ? (e = this.numerator, t = this.denominator.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) {
   const e = new u().one(), t = o.factors.reduce((i, n) => {
     if (n.polynom.degree().isZero())
@@ -2820,7 +2826,7 @@ function ei(o) {
     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 ue(e.display, 1)]), t;
+  return e.isOne() || (t[e.display] = [new le(e.display, 1)]), t;
 }
 function ts(o, e) {
   return o.dimension === e.dimension && o.array.every(
@@ -3041,22 +3047,22 @@ var we = /* @__PURE__ */ ((o) => (o.None = "none", o.Parallel = "parallel", o.Pe
 function qi(o = 0.5) {
   return Math.random() < o;
 }
-function fe(o, e, t) {
+function de(o, e, t) {
   if (e === void 0)
-    return o >= 0 ? fe(0, o) : fe(o, 0);
+    return o >= 0 ? de(0, o) : de(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 = fe(o, e);
+  let i = de(o, e);
   for (; t.includes(i); )
-    i = fe(o, e);
+    i = de(o, e);
   return i;
 }
 function Y(o, e) {
-  return e === !1 ? qi() ? fe(1, o) : -fe(1, o) : fe(-o, o);
+  return e === !1 ? qi() ? de(1, o) : -de(1, o) : de(-o, o);
 }
 function rs(o) {
   let e = G.primes();
@@ -3066,7 +3072,7 @@ function os(o, e) {
   return e === void 0 && (e = 1), o.length <= 0 ? Object.values(o) : ki(o).slice(0, e);
 }
 function pi(o) {
-  return o.length === 0 ? null : o[fe(0, o.length - 1)];
+  return o.length === 0 ? null : o[de(0, o.length - 1)];
 }
 function ki(o) {
   const e = Object.values(o);
@@ -3102,7 +3108,7 @@ class B extends x {
     return e.array = this.copy(), e.asPoint = !0, e;
   }
 }
-var Ie, S, z, U, ne, Q, Pe, be;
+var Pe, S, z, U, ne, Q, Ie, be;
 const Ze = class Ze {
   /**
    * Value can be a mix of:
@@ -3110,14 +3116,14 @@ const Ze = class Ze {
    * @param values
    */
   constructor(...e) {
-    p(this, Ie);
+    p(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, Ie);
     p(this, be, "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));
@@ -3189,19 +3195,19 @@ 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, 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, Ie, s(this, Q).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, ne, e.clone()), h(this, Q, t.clone()), h(this, Ie, s(this, Q).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, U, s(this, U).clone()), h(this, Q, s(this, Q).clone()), h(this, ne, s(this, ne).clone()), h(this, Ie, s(this, Ie).clone()), this));
     // ------------------------------------------
     // Mathematical operations
     // ------------------------------------------
@@ -3239,7 +3245,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, U, new u().zero()), h(this, ne, new x()), h(this, Q, new x()), h(this, Ie, new x()), h(this, Pe, !0), e.length > 0 && this.parse(...e), this;
   }
   get a() {
     return s(this, S);
@@ -3274,12 +3280,12 @@ const Ze = class Ze {
     h(this, Q, e);
   }
   get n() {
-    return s(this, Pe);
+    return s(this, Ie);
   }
   // ------------------------------------------
   getEquation() {
-    const e = new H(new q().parse("xy", s(this, S), s(this, z), s(this, U)), new q("0"));
-    return s(this, Ie) ? e.simplify() : e;
+    const e = new H(new O().parse("xy", s(this, S), s(this, z), s(this, U)), new O("0"));
+    return s(this, Pe) ? e.simplify() : e;
   }
   // get system(): { x: Equation, y: Equation } {
   //     const e1 = new Equation(
@@ -3315,13 +3321,13 @@ const Ze = class Ze {
       case "equation":
         return this.getEquation().reorder().tex;
       case "mxh":
-        return this.slope.isInfinity() ? "x=" + this.OA.x.tex : "y=" + new q().parse("x", this.slope, this.height).tex;
+        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 q(s(this, ne).x).add(new M(s(this, Q).x).multiply(new M("k"))).reorder("k", !0).tex}\\\\ 
-            y &= ${new q(s(this, ne).y).add(new M(s(this, Q).y).multiply(new M("k"))).reorder("k", !0).tex}
+        return s(this, Pe) && 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}
             \\end{aligned}\\right.`;
       }
       default: {
@@ -3331,10 +3337,10 @@ 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);
@@ -3342,10 +3348,10 @@ const Ze = class Ze {
       case "equation":
         return this.getEquation().reorder().display;
       case "mxh":
-        return this.slope.isInfinity() ? "x=" + this.OA.x.display : "y=" + new q().parse("x", this.slope, this.height).display;
+        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}))`;
+        return s(this, Pe) && t.simplify(), `((x,y))=((${s(this, ne).x.display},${s(this, ne).y.display}))+k((${t.x.display},${t.y.display}))`;
       }
       default: {
         const t = this.getEquation();
@@ -3402,10 +3408,10 @@ const Ze = class Ze {
     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";
   }
 };
-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
+Pe = new WeakMap(), S = new WeakMap(), z = new WeakMap(), U = new WeakMap(), ne = new WeakMap(), Q = new WeakMap(), Ie = 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;
+var re, j, qe, jt, Ut, Gt, oe, Mi, Et, Ci, $i, Pi, ai;
 const Wt = class Wt {
   constructor(...e) {
     p(this, oe);
@@ -3452,7 +3458,7 @@ const Wt = class Wt {
       return [new R(e, t, we.Perpendicular)];
     });
     p(this, Ut, (e) => {
-      const t = this.center.x.clone().subtract(e.x), i = this.center.y.clone().subtract(e.y), n = new q("x"), r = new q("x^2+1");
+      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");
@@ -3506,10 +3512,10 @@ const Wt = class Wt {
     );
   }
   setRadius(e, t) {
-    return t ? h(this, j, new u(e)) : h(this, j, new u(e).pow(2)), O(this, oe, Et).call(this), this;
+    return t ? h(this, j, new u(e)) : h(this, j, new u(e).pow(2)), q(this, oe, Et).call(this), this;
   }
   parse(...e) {
-    return O(this, oe, Mi).call(this), typeof e[0] == "string" ? O(this, oe, ai).call(this, new H(e[0])) : e[0] instanceof H ? O(this, oe, ai).call(this, e[0]) : e[0] instanceof Wt ? O(this, oe, Ci).call(this, e[0]) : e[0] instanceof B && e.length > 1 && (e[1] instanceof B ? e[2] instanceof B || O(this, oe, Ii).call(this, e[0], e[1]) : (e[1] instanceof u || typeof e[1] == "number") && O(this, oe, $i).call(this, e[0], e[1], typeof e[2] == "boolean" ? e[2] : !1)), O(this, oe, Et).call(this), this;
+    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, Pi).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;
   }
   // private _parseThroughtThreePoints(A: Point, B: Point, C: Point): this {
   //     const T = new Triangle(A, B, C), mAB = T.remarquables.mediators.AB.clone(),
@@ -3522,14 +3528,14 @@ re = new WeakMap(), j = new WeakMap(), qe = new WeakMap(), jt = new WeakMap(), U
   return h(this, re, void 0), h(this, j, void 0), h(this, qe, void 0), this;
 }, Et = function() {
   h(this, qe, new H(
-    new q(`(x-(${this.center.x.display}))^2+(y-(${this.center.y.display}))^2`),
-    new q(this.squareRadius.display)
+    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()), O(this, oe, Et).call(this), this;
+  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) {
+}, Pi = function(e, t) {
   return h(this, re, e.clone()), h(this, j, new x(s(this, re), t).normSquare), this;
 }, ai = function(e) {
   if (e.moveLeft(), e.degree("x").value === 2 && e.degree("y").value === 2) {
@@ -3614,11 +3620,11 @@ const lt = class lt {
     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)}`,
       system: `\\left\\{\\begin{aligned}
-    x &= ${new q(s(this, L).x).add(new M(s(this, W).x).multiply(new M("k"))).reorder("k", !0).tex}\\\\ 
-    y &= ${new q(s(this, L).y).add(new M(s(this, W).y).multiply(new M("k"))).reorder("k", !0).tex}\\\\
-    z &= ${new q(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 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}
 \\end{aligned}\\right.`,
-      cartesian: `\\frac{ ${new q("x", 1, s(this, L).x.clone().opposite()).tex} }{ ${this.direction.x.tex} } = \\frac{ ${new q("y", 1, s(this, L).y.clone().opposite()).tex} }{ ${this.direction.y.tex} } = \\frac{ ${new q("z", 1, s(this, L).z.clone().opposite()).tex} }{ ${this.direction.z.tex} }`
+      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} }`
     };
   }
   get display() {
@@ -3708,8 +3714,8 @@ const gi = class gi {
   }
   get tex() {
     return new H(
-      new q("xyz", this.a, this.b, this.c, this.d),
-      new q(0)
+      new O("xyz", this.a, this.b, this.c, this.d),
+      new O(0)
     ).reduce().tex;
   }
   parse(e) {
@@ -3984,16 +3990,16 @@ function Ot(o) {
     },
     o
   ), t = new u();
-  if (e.negative ? t.numerator = Y(e.max, e.zero) : t.numerator = fe(e.zero ? 0 : 1, e.max), e.natural)
+  if (e.negative ? t.numerator = Y(e.max, e.zero) : t.numerator = de(e.zero ? 0 : 1, e.max), e.natural)
     t.denominator = 1;
   else {
     let i = 0;
     for (; t.isRelative() && i < 10; )
-      t.denominator = fe(1, e.max), i++;
+      t.denominator = de(1, e.max), i++;
   }
   return e.reduced ? t.reduce() : t;
 }
-function Pi(o) {
+function Ii(o) {
   const e = Object.assign(
     {
       letters: "x",
@@ -4033,10 +4039,10 @@ function Bi(o) {
   const e = Object.assign(
     as,
     o
-  ), t = new q().empty();
+  ), t = new O().empty();
   let i;
   for (let n = e.degree; n >= 0; n--)
-    i = Pi({
+    i = Ii({
       letters: e.letters,
       degree: n,
       fraction: e.fraction,
@@ -4044,7 +4050,7 @@ 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 = fe(1, t.length - 1);
+      const n = de(1, t.length - 1);
       t.monoms.splice(n, 1);
     }
   return t.reduce();
@@ -4069,7 +4075,7 @@ function ls(o) {
       }
     },
     o
-  ), t = new q().one();
+  ), t = new O().one();
   for (let i = 0; i < e.degree; i++) {
     const n = Bi({
       degree: 1,
@@ -4106,7 +4112,7 @@ function cs(o) {
     o
   ), t = ui(e.center);
   let i, n;
-  return e.pointsOnCircle === 8 ? (i = fe(1, 3), n = i ** 2 + (i + 1) ** 2) : n = fe(1, 20), new Tt(t, n, !0);
+  return e.pointsOnCircle === 8 ? (i = de(1, 3), n = i ** 2 + (i + 1) ** 2) : n = de(1, 20), new Tt(t, n, !0);
 }
 function us(o) {
   const e = Object.assign(
@@ -4146,9 +4152,9 @@ function fs(o) {
 const ds = {
   equation: (o) => ls(o),
   polynom: (o) => Bi(o),
-  monom: (o) => Pi(o),
+  monom: (o) => Ii(o),
   fraction: (o) => Ot(o),
-  number: (o, e, t) => fe(o, e, t),
+  number: (o, e, t) => de(o, e, t),
   numberSym: (o, e) => Y(o, e),
   prime: (o) => rs(o),
   bool: (o) => qi(o),
@@ -4168,11 +4174,11 @@ const ds = {
   Fraction: u,
   Root: vt,
   Monom: M,
-  Polynom: q,
+  Polynom: O,
   Equation: H,
   Matrix: hs,
   LinearSystem: ii,
-  Factor: ue,
+  Factor: le,
   PolyFactor: ni,
   LogicalSet: es,
   Random: ds,
@@ -4192,7 +4198,7 @@ export {
   H as Equation,
   Nt as EquationSolver,
   at as FACTOR_DISPLAY,
-  ue as Factor,
+  le as Factor,
   u as Fraction,
   R as Line,
   At as Line3,
@@ -4201,10 +4207,13 @@ export {
   hs as Matrix,
   M as Monom,
   vt as NthRoot,
+  _i as NumExp,
+  G as Numeric,
   li as Plane3,
   B as Point,
   ni as PolyFactor,
-  q as Polynom,
+  O as Polynom,
+  ds as Random,
   ci as Triangle,
   x as Vector,
   is as areVectorsColinears,