pimath 0.0.91 → 0.0.92

package/src/maths/algebra/monom.ts

@@ -825,6 +825,31 @@ export class Monom {
         return r;
     };
 
+    evaluateAsNumeric = (values: { [Key: string]: number } | number): number => {
+        let r = this.coefficient.value
+
+        if (typeof values === 'number') {
+            let tmpValues: { [Key: string]: number } = {}
+            tmpValues[this.variables[0]] = values
+            return this.evaluateAsNumeric(tmpValues);
+        }
+
+        if (typeof values === 'object') {
+            if (this.variables.length === 0) {
+                return this.coefficient.value
+            }
+            for (let L in this._literal) {
+                if (values[L] === undefined) {
+                    return 0;
+                }
+
+                r *= values[L] ** (this._literal[L].value)
+            }
+        }
+
+        return r
+    }
+
     /**
      * Derivative the monom
      * @param letter

package/src/maths/algebra/polynom.ts

@@ -817,6 +817,15 @@ export class Polynom {
         return r;
     };
 
+    evaluateAsNumeric = (values: { [Key: string]: number } | number): number => {
+        let r = 0
+        this._monoms.forEach(monom => {
+            r += monom.evaluateAsNumeric(values)
+        })
+
+        return r
+    }
+
     derivative = (letter?: string): Polynom => {
         let dP = new Polynom();
 

package/src/maths/algebra/rational.ts

@@ -217,6 +217,10 @@ export class Rational {
         return N.divide(D)
     };
 
+    evaluateAsNumeric = (values: { [Key: string]: number } | number): number => {
+        return this._numerator.evaluateAsNumeric(values) / this._denominator.evaluateAsNumeric(values)
+    }
+
     study = (): RationalStudy => {
         return new RationalStudy(this)
     }

package/src/maths/algebra/study/rationalStudy.ts

@@ -15,6 +15,7 @@
  */
 import {
     ASYMPTOTE,
+    ASYMPTOTE_POSITION,
     FUNCTION_EXTREMA,
     IAsymptote,
     ITableOfSigns,
@@ -40,8 +41,7 @@ export class RationalStudy extends Study {
     };
 
     makeSigns(): ITableOfSigns {
-        let tos = this._getSigns(this.fx, this.zeroes)
-        return tos
+        return this._getSigns(this.fx, this.zeroes)
     };
 
     makeAsymptotes(): IAsymptote[] {
@@ -51,9 +51,10 @@ export class RationalStudy extends Study {
         let asymptotes: IAsymptote[] = []
         this.zeroes.filter(x => x.type === ZEROTYPE.DEFENCE).forEach(zero => {
             // Check if it's a hole or an asymptote
-            // TODO: Check for a hole ! Means calculate the limits !
             let Ztype = ASYMPTOTE.VERTICAL,
                 tex = `x=${zero.tex}`
+
+            // Check if it's a hole: the reduced polynom should not be null
             if (zero.exact instanceof Fraction) {
                 if (reduced.denominator.evaluate(zero.exact).isNotZero()) {
                     Ztype = ASYMPTOTE.HOLE
@@ -66,14 +67,48 @@ export class RationalStudy extends Study {
                 }
             }
 
+            // Get the position before and after the asymptote.
+            const delta = 0.000001
+            let before = this.fx.evaluateAsNumeric(zero.value - delta),
+                after = this.fx.evaluateAsNumeric(zero.value + delta),
+                position: ASYMPTOTE_POSITION[] = [],
+                pm = ""
+
+            if (after < -10000) {
+                position.push(ASYMPTOTE_POSITION.RB)
+                pm += "m"
+            } else if (after > 10000) {
+                position.push(ASYMPTOTE_POSITION.RT)
+                pm += "p"
+            }
+
+            if (before < -10000) {
+                position.push(ASYMPTOTE_POSITION.LB)
+                pm += "m"
+            } else if (before > 10000) {
+                position.push(ASYMPTOTE_POSITION.LT)
+                pm += "p"
+            }
+
+            // Left and right are to infinity
+            // TODO: handle the case were one side of the asymptote isn't infinity (not possible in rational study?!)
+            if (pm === "pp") {
+                pm = "+"
+            } else if (pm === "mm") {
+                pm = "-"
+            } else {
+                pm = `\\${pm}`
+            }
+
             asymptotes.push({
                 fx: null,
                 type: Ztype,
                 tex: tex,
                 zero: zero,
-                limits: `\\lim_{x\\to${zero.tex} }\\ f(x) = \\pm\\infty`,
+                limits: `\\lim_{x\\to${zero.tex} }\\ f(x) = ${pm}\\infty`,
                 deltaX: null,
-                tableOfSign: null
+                tableOfSign: null,
+                position
             })
         })
 
@@ -87,7 +122,7 @@ export class RationalStudy extends Study {
             let {reminder} = reduced.euclidian(),
                 deltaX = new Rational(reminder, reduced.denominator)
 
-
+            // Determine the position above or below on the left / right of the asymptote.
             asymptotes.push({
                 fx: new Polynom(H),
                 type: ASYMPTOTE.HORIZONTAL,
@@ -95,7 +130,8 @@ export class RationalStudy extends Study {
                 zero: null,
                 limits: `\\lim_{x\\to\\infty}\\ f(x) = ${Htex}`,
                 deltaX,
-                tableOfSign: this._getSigns(deltaX)
+                tableOfSign: this._getSigns(deltaX),
+                position: this._getHorizontalAsymptoteRelativePositon(deltaX)
             })
         } else if (DDegree.greater(NDegree)) {
             asymptotes.push({
@@ -105,7 +141,8 @@ export class RationalStudy extends Study {
                 zero: null,
                 limits: `\\lim_{x\\to\\infty}\\ f(x) = ${0}`,
                 deltaX: null,
-                tableOfSign: null
+                tableOfSign: null,
+                position: this._getHorizontalAsymptoteRelativePositon(this.fx)
             })
         } else if (NDegree.value - 1 === DDegree.value) {
             // Calculate the slope
@@ -119,13 +156,35 @@ export class RationalStudy extends Study {
                 zero: null,
                 limits: ``,
                 deltaX: new Rational(reminder, reduced.denominator),
-                tableOfSign: this._getSigns(deltaX)
+                tableOfSign: this._getSigns(deltaX),
+                position: this._getHorizontalAsymptoteRelativePositon(deltaX)
             })
         }
 
         return asymptotes
     };
 
+    _getHorizontalAsymptoteRelativePositon(deltaX: Rational, delta: number = 1000000): ASYMPTOTE_POSITION[] {
+
+        let position: ASYMPTOTE_POSITION[] = [],
+            before = deltaX.evaluateAsNumeric(-delta),
+            after = deltaX.evaluateAsNumeric(delta)
+
+        if (before >= 0) {
+            position.push(ASYMPTOTE_POSITION.LT)
+        } else {
+            position.push(ASYMPTOTE_POSITION.LB)
+        }
+
+        if (after >= 0) {
+            position.push(ASYMPTOTE_POSITION.RT)
+        } else {
+            position.push(ASYMPTOTE_POSITION.RB)
+        }
+
+        return position
+    }
+
     makeDerivative(): ITableOfSigns {
         let dx = this.fx.clone().derivative(),
             tos = this._getSigns(dx, this._getZeroes(dx), TABLE_OF_SIGNS.GROWS)

package/src/maths/algebra/study.ts

@@ -31,6 +31,12 @@ export enum ASYMPTOTE {
     HOLE = "hole"
 }
 
+export enum ASYMPTOTE_POSITION {
+    "LT" = "LT",
+    "RT" = "RT",
+    "LB" = "LB",
+    "RB" = "RB"
+}
 export interface IAsymptote {
     fx: Polynom,
     deltaX: StudyableFunction
@@ -38,7 +44,8 @@ export interface IAsymptote {
     tex: string,
     type: ASYMPTOTE,
     zero: IZero,
-    tableOfSign: ITableOfSigns
+    position: ASYMPTOTE_POSITION[]
+    tableOfSign: ITableOfSigns,
 }
 
 export enum FUNCTION_EXTREMA {
@@ -74,7 +81,7 @@ export interface ITableOfSigns {
 export enum TABLE_OF_SIGNS {
     SIGNS = "signs",
     GROWS = "grows",
-    VARIATIONS = "variatins"
+    VARIATIONS = "variations"
 }
 
 /**
@@ -150,7 +157,9 @@ export class Study {
         this._variations = this.makeVariation()
 
         this._signs.tex = this.texSigns
+
         this._derivative.tex = this.texGrows
+
         this._variations.tex = this.texVariations
 
     };

package/tests/algebra/study.test.ts

@@ -10,16 +10,6 @@ describe('Study tests', () => {
             // new Rational('x^2-4x-4', 'x+7')
             new Rational('(3x-2)(x-3)(x+4)', 'x^2-5x+6')
         )
-
-        // console.log(study.texSigns)
-        // console.log(study.asymptotes)
-        // console.log(study.derivative.fx.texFactors)
-        // console.log(study.texGrows)
-        //
-        //
-        // console.log('----------------')
-
-
         let AO = study.asymptotes.filter(x => x.type === ASYMPTOTE.SLOPE)[0]
 
         console.log(AO.tableOfSign.signs)
@@ -40,4 +30,12 @@ describe('Study tests', () => {
             "Z_9(1.3333333333333333,0)*")
     });
 
+    it('should get the before/after state of asymptotes', function () {
+        const study = new RationalStudy(
+            new Rational("x+5", "x-3")
+            // new Rational("x^2+5", "(x-3)^2")
+        )
+
+        expect(study.asymptotes[0].position).to.have.all.members(["LB", "RT"])
+    });
 })