AxiomX 0.1.2__tar.gz → 0.1.3__tar.gz

axiomx-0.1.3/AxiomX/constants.py

@@ -0,0 +1,35 @@
+from AxiomX.calculus import *
+from AxiomX.functions import *
+
+pi = 3.141592653589793
+e = 2.718281828459045
+tau = 2*pi
+lemniscate = 2.622057554292119
+gauss = lemniscate / pi
+euler_mascheroni = 0.577215664901533
+sqrt_2 = 2**0.5
+sqrt_3 = 3**0.5
+sqrt_5 = 5**0.5
+golden_ratio = (1 + sqrt_5) / 2
+silver_ratio = 1 + sqrt_2
+gelfond = e**pi
+gelfond_schneider = 2**sqrt_2
+infinity = float("inf")
+
+def metallic_ratio(n):
+    return (n + (n**2 + 4)**0.5) / 2
+    
+def harmonic_number(n):
+    return summation(1, n, lambda x:1/x)
+    
+def bernoulli(n):
+    B = [0] * (n + 1)
+    B[0] = 1
+
+    for m in range(1, n + 1):
+        B[m] = 0
+        for k in range(m):
+            B[m] -= (bc(m + 1, k) * B[k])
+        B[m] /= (m + 1)
+
+    return B[n]

{axiomx-0.1.2 → axiomx-0.1.3}/AxiomX/exp.py

@@ -1,3 +1,5 @@
+from AxiomX.constants import *
+
 def exp(n):
     return e**n
 

{axiomx-0.1.2 → axiomx-0.1.3}/AxiomX/functions.py

@@ -1,3 +1,5 @@
+from AxiomX.constants import *
+
 def absolute(x):
     if (x > 0):
         return x
@@ -75,7 +77,7 @@ def zeta(n):
     zetval = 0
     if n <= 1:
         raise ValueError("zeta(n) diverges for n <= 1")
-    for _ in range(1, 100001):
+    for _ in range(1, 1000001):
         zetval += (1 / _**n)
     return zetval
     
@@ -83,7 +85,7 @@ def beta(n):
     if n == 0:
         return 0.5
     total = 0.0
-    for i in range(100000):
+    for i in range(1000000):
         total += ((-1)**i) / ((2*i + 1)**n)
     return total
     
@@ -103,4 +105,12 @@ def pascal(row):
     pr = []
     for i in range(0, row+1):
         pr.append(int(bc(row, i)))
-    return pr
+    return pr
+
+def perm(n, r):
+    if n < r or n < 0 or r < 0:
+        return 0
+    result = 1
+    for i in range(r):
+        result *= (n - i)
+    return result

axiomx-0.1.3/AxiomX/graph.py

@@ -0,0 +1,112 @@
+import turtle
+
+def plot(funcs, start, end, step=0.1, scale=20, limit=1000):
+    if not isinstance(funcs, list):
+        funcs = [funcs]
+    screen = turtle.Screen()
+    screen.title("AxiomX Graph")
+
+    t = turtle.Turtle()
+    t.speed(0)
+    t.hideturtle()
+
+# =========================
+# AXES + NUMBER LINES
+# =========================
+
+    # Axes
+    t.penup()
+    t.goto(-400, 0)
+    t.pendown()
+    t.goto(400, 0)
+
+    t.penup()
+    t.goto(0, -300)
+    t.pendown()
+    t.goto(0, 300)
+
+    # ----- X-axis ticks -----
+    # ----- X-axis ticks (with proper negative handling) -----
+
+    tick_step = max(1, int((end - start) / 10))
+
+    x_tick = int(start // tick_step) * tick_step
+
+    while x_tick <= end:
+        sx = x_tick * scale
+
+        # Tick
+        t.penup()
+        t.goto(sx, -5)
+        t.pendown()
+        t.goto(sx, 5)
+
+        # Label
+        t.penup()
+        t.goto(sx, -20)
+
+        label = "0" if abs(x_tick) < 1e-9 else str(x_tick)
+        t.write(label, align="center", font=("Arial", 8, "normal"))
+
+        x_tick += tick_step
+
+    # ----- Y-axis ticks -----
+    # Estimate visible y-range based on screen height
+    y_range = 300 / scale
+    y_tick_step = max(1, int((2 * y_range) / 10))
+
+    y_tick = int(-y_range)
+    while y_tick <= y_range:
+        sy = y_tick * scale
+
+        # Skip origin label overlap (optional)
+        if y_tick == 0:
+            y_tick += y_tick_step
+            continue
+
+        # Tick
+        t.penup()
+        t.goto(-5, sy)
+        t.pendown()
+        t.goto(5, sy)
+
+        # Label
+        t.penup()
+        t.goto(-25, sy - 5)
+        t.write(f"{y_tick}", align="right", font=("Arial", 8, "normal"))
+
+        y_tick += y_tick_step
+    colors = ["red", "blue", "green", "purple"]
+
+    for i, f in enumerate(funcs):
+        t.color(colors[i % len(colors)])
+        t.penup()
+
+        x = start
+        drawing = False
+
+        while x <= end:
+            try:
+                y = f(x)
+
+                if y is None or abs(y) > limit:
+                    raise ValueError
+
+                sx = x * scale
+                sy = y * scale
+
+                if not drawing:
+                    t.penup()
+                    t.goto(sx, sy)
+                    t.pendown()
+                    drawing = True
+                else:
+                    t.goto(sx, sy)
+
+            except turtle.Terminator:
+                return
+            except:
+                t.penup()
+                drawing = False
+
+            x += step

{axiomx-0.1.2 → axiomx-0.1.3}/AxiomX/hyperbolic.py

@@ -1,3 +1,6 @@
+from AxiomX.exp import ln
+from AxiomX.constants import e
+
 def sinh(x):
     return (e**x - e**(-x))/2
     
@@ -17,13 +20,13 @@ def cosech(x):
     return 1 / sinh(x)
     
 def arcsinh(x):
-    return ln(x + sqrt(x**2 + 1))
+    return ln(x + (x**2 + 1)**0.5)
     
 def arccosh(x):
-    return abs(arcsinh(sqrt(x**2 - 1)))
+    return abs(arcsinh((x**2 - 1)**0.5))
     
 def arccoth(x):
-    return 0.5 * ((x+1)/(x-1))
+    return 0.5 * ln((x+1)/(x-1))
     
 def arctanh(x):
     return arccoth(1/x)

{axiomx-0.1.2 → axiomx-0.1.3}/AxiomX/matrix.py

@@ -1,11 +1,17 @@
 class Matrix:
     def __str__(self):
+        """
+        Makes it suitable for printing.
+        """
         result = ""
         for row in self.matrix:
             result += " ".join(str(x) for x in row) + "\n"
         return result
     
     def __init__(self, rows, columns):
+        """
+        Initializes a matrix.
+        """
         self.rows = rows
         self.columns = columns
         self.order = (rows,columns)
@@ -19,6 +25,9 @@ class Matrix:
         
     @classmethod
     def from_list(cls, data):
+        """
+        This class method helps in loading a matrix from its list form.
+        """
         rows = len(data)
         columns = len(data[0]) if rows > 0 else 0
 
@@ -27,16 +36,26 @@ class Matrix:
         return obj
         
     def transpose(self):
+        """
+        This method helps in transposing a matrix.
+        """
         data = [[self.matrix[i][j] for i in range(self.rows)] for j in range(self.columns)]
         return Matrix.from_list(data)
         
     def zero(self):
+        """
+        This method makes a matrix a zero matrix, whether the matrix was already having values or not.
+        """
         self.matrix = []
         self.matrix = [[0 for _ in range(self.columns)] for _ in range(self.rows)]
 
     def identity(self):
+        """
+        This method makes a matrix an identity matrix, whether the matrix was already having values or not.
+        """
         self.matrix = []  # reset matrix
-
+        if (self.rows != self.columns):
+            raise ValueError("Identity matrix is a square matrix")
         for i in range(self.rows):
             row = []
             for j in range(self.columns):
@@ -47,12 +66,18 @@ class Matrix:
             self.matrix.append(row)
             
     def __getitem__(self, key):
+        """
+        This method helps in indexing.
+        """
         if isinstance(key, tuple):
             i, j = key
             return self.matrix[i][j]
         return self.matrix[key]
 
     def __setitem__(self, key, value):
+        """
+        This method helps in setting values with a particular index.
+        """
         if isinstance(key, tuple):
             i, j = key
             self.matrix[i][j] = value
@@ -60,6 +85,9 @@ class Matrix:
             self.matrix[key] = value
             
     def __add__(self, other):
+        """
+        This method overloads the + operator for adding two matrices.
+        """
         if self.rows != other.rows or self.columns != other.columns:
             raise ValueError("matrices cannot be added")
 
@@ -71,6 +99,9 @@ class Matrix:
         return Matrix.from_list(data)
         
     def __sub__(self, other):
+        """
+        This method overloads the - operator to subtract two matrices.
+        """
         if self.rows != other.rows or self.columns != other.columns:
             raise ValueError("matrices cannot be added")
 
@@ -82,6 +113,9 @@ class Matrix:
         return Matrix.from_list(data)
         
     def __mul__(self, other):
+        """
+        This method overloads the * operator for multiplying two matrices.
+        """
         if isinstance(other, (int, float)):
             return Matrix.from_list([
                 [self[i, j] * other for j in range(self.columns)]
@@ -102,6 +136,9 @@ class Matrix:
         return Matrix.from_list(data)
         
     def det(self):
+        """
+        This method calculates the determinant of any square matrix.
+        """
         if self.rows != self.columns:
             raise ValueError("determinant is only defined for square matrices")
 
@@ -131,12 +168,18 @@ class Matrix:
         return determinant
     
     def _minor(self, i, j):
+        """
+        This method returns the minor of any element in the matrix.
+        """
         return Matrix.from_list([
             [self[r, c] for c in range(self.columns) if c != j]
             for r in range(self.rows) if r != i
         ])
         
     def _cofactor_matrix(self):
+        """
+        This method returns the cofactor matrix of a matrix.
+        """
         data = []
 
         for i in range(self.rows):
@@ -150,6 +193,9 @@ class Matrix:
         return Matrix.from_list(data)
         
     def inverse(self):
+        """
+        This method returns the inverse of a matrix.
+        """
         if self.rows != self.columns:
             raise ValueError("Inverse only defined for square matrices")
 
@@ -169,9 +215,15 @@ class Matrix:
         return adjoint*(1 / det)
         
     def __truediv__(self, other):
+        """
+        This method overloads the / operator for division of two matrices.
+        """
         return self*other.inverse()
         
     def __pow__(self, other):
+        """
+        This method overloads the ** operator to multiply the matrix n times (works only for square matrices).
+        """
         if not isinstance(other, int):
             raise ValueError("exponent is not a number")
         else:
@@ -179,4 +231,38 @@ class Matrix:
             p.identity()
             for i in range(other):
                 p *= self
-            return p
+            return p
+    
+    def rank(self):
+        A = [row[:] for row in self]  # copy
+        n = len(A)
+        m = len(A[0])
+        
+        rank = 0
+        row = 0
+
+        for col in range(m):
+            pivot = None
+            for r in range(row, n):
+                if A[r][col] != 0:
+                    pivot = r
+                    break
+
+            if pivot is not None:
+                A[row], A[pivot] = A[pivot], A[row]
+
+                pivot_val = A[row][col]
+                for j in range(col, m):
+                    A[row][j] /= pivot_val
+
+                for r in range(n):
+                    if r != row:
+                        factor = A[r][col]
+                        for j in range(col, m):
+                            A[r][j] -= factor * A[row][j]
+
+                row += 1
+                rank += 1
+
+        return rank
+    

{axiomx-0.1.2 → axiomx-0.1.3}/AxiomX.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: AxiomX
-Version: 0.1.2
+Version: 0.1.3
 Summary: A dynamic Python math library containing numerous mathematical functions implemented from scratch.
 Author: Eric Joseph
 License: MIT
@@ -20,18 +20,18 @@ exploration.
 
 ## 📢 Community
 
-Join Discord for updates and feature requests. Join [here](https://discord.gg/CrEf8mDB).
+Join Discord for updates and feature requests. Join [here](https://discord.gg/yhN5eYwT).
 
 ------------------------------------------------------------------------
 
 ## 🚀 Features
 
-### 📐 Core Mathematics
+### 📐 Core Mathematics (AxiomX.constants)
 
 -   Constants (π, e, φ, √2, γ, ...)
 -   Custom implementations
 
-### 📈 Exponential & Logarithmic
+### 📈 Exponential & Logarithmic (AxiomX.exp)
 
 -   exp(x)
 -   ln(x)
@@ -39,22 +39,25 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   log2(x)
 -   log(arg, base)
 
-### 📐 Trigonometry
+### 📐 Trigonometry (AxiomX.trig)
 
 -   sin, cos, tan, cot, sec, cosec
 -   inverse functions
 
-### 🔢 Hyperbolic Functions
+### 🔢 Hyperbolic Functions (AxiomX.hyperbolic)
 
 -   sinh, cosh, tanh, coth
 
-### 🧮 General Functions
+### 🧮 General Functions (AxiomX.functions)
 
 -   absolute(x)
 -   sqrt(x)
 -   gamma(x)
 -   cbrt(x)
 -   zeta(n) and beta(n)
+-   bc(n, k) - **NEW**
+-   perm(n, K) - **NEW**
+-   pascal(row) - **NEW**
 
 ### ∫ Calculus
 
@@ -63,9 +66,7 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   convergence analysis
 -   limit evaluation
 
-------------------------------------------------------------------------
-
-## 🧱 Matrix Engine (NEW)
+### 🧱 Matrix Engine
 
 -   Matrix creation
 -   A\[i, j\] access
@@ -74,6 +75,7 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   Determinant
 -   Inverse
 -   Division (via inverse) by using / operator
+-   Rank of a Matrix by Gaussian elimination (**NEW**)
 
 If you want to create a zero matrix or identity of any order, you can do it this way:
 
@@ -94,6 +96,16 @@ a = Matrix.from_list([]) # enter your matrix in it.
 ```
 ------------------------------------------------------------------------
 
+## 📈 GRAPH ENGINE - **NEW** (included in AxiomX.graph)
+
+AxiomX has introduced a killer feature - graphing. Just like how `matplotlib` is for numpy, AxiomX has introduced a graph module. It has only one function:
+
+- plot(functions, start, end)
+
+The functions argument can a function or list of functions. It also contains a scale and you can provide a starting and ending limit. We will add zooming later.
+
+------------------------------------------------------------------------
+
 ## 📦 Installation
 
 `pip install axiomx`
@@ -104,7 +116,8 @@ a = Matrix.from_list([]) # enter your matrix in it.
 ```
 from AxiomX.matrix import Matrix
 
-A = Matrix.from_list(\[\[1, 2\], \[3, 4\]\]) print(A.det())
+A = Matrix.from_list([1, 2], [3, 4]])
+print(A.det())
 ```
 
 ------------------------------------------------------------------------

{axiomx-0.1.2 → axiomx-0.1.3}/AxiomX.egg-info/SOURCES.txt

@@ -7,6 +7,7 @@ AxiomX/calculus.py
 AxiomX/constants.py
 AxiomX/exp.py
 AxiomX/functions.py
+AxiomX/graph.py
 AxiomX/hyperbolic.py
 AxiomX/matrix.py
 AxiomX/trig.py

{axiomx-0.1.2 → axiomx-0.1.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: AxiomX
-Version: 0.1.2
+Version: 0.1.3
 Summary: A dynamic Python math library containing numerous mathematical functions implemented from scratch.
 Author: Eric Joseph
 License: MIT
@@ -20,18 +20,18 @@ exploration.
 
 ## 📢 Community
 
-Join Discord for updates and feature requests. Join [here](https://discord.gg/CrEf8mDB).
+Join Discord for updates and feature requests. Join [here](https://discord.gg/yhN5eYwT).
 
 ------------------------------------------------------------------------
 
 ## 🚀 Features
 
-### 📐 Core Mathematics
+### 📐 Core Mathematics (AxiomX.constants)
 
 -   Constants (π, e, φ, √2, γ, ...)
 -   Custom implementations
 
-### 📈 Exponential & Logarithmic
+### 📈 Exponential & Logarithmic (AxiomX.exp)
 
 -   exp(x)
 -   ln(x)
@@ -39,22 +39,25 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   log2(x)
 -   log(arg, base)
 
-### 📐 Trigonometry
+### 📐 Trigonometry (AxiomX.trig)
 
 -   sin, cos, tan, cot, sec, cosec
 -   inverse functions
 
-### 🔢 Hyperbolic Functions
+### 🔢 Hyperbolic Functions (AxiomX.hyperbolic)
 
 -   sinh, cosh, tanh, coth
 
-### 🧮 General Functions
+### 🧮 General Functions (AxiomX.functions)
 
 -   absolute(x)
 -   sqrt(x)
 -   gamma(x)
 -   cbrt(x)
 -   zeta(n) and beta(n)
+-   bc(n, k) - **NEW**
+-   perm(n, K) - **NEW**
+-   pascal(row) - **NEW**
 
 ### ∫ Calculus
 
@@ -63,9 +66,7 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   convergence analysis
 -   limit evaluation
 
-------------------------------------------------------------------------
-
-## 🧱 Matrix Engine (NEW)
+### 🧱 Matrix Engine
 
 -   Matrix creation
 -   A\[i, j\] access
@@ -74,6 +75,7 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   Determinant
 -   Inverse
 -   Division (via inverse) by using / operator
+-   Rank of a Matrix by Gaussian elimination (**NEW**)
 
 If you want to create a zero matrix or identity of any order, you can do it this way:
 
@@ -94,6 +96,16 @@ a = Matrix.from_list([]) # enter your matrix in it.
 ```
 ------------------------------------------------------------------------
 
+## 📈 GRAPH ENGINE - **NEW** (included in AxiomX.graph)
+
+AxiomX has introduced a killer feature - graphing. Just like how `matplotlib` is for numpy, AxiomX has introduced a graph module. It has only one function:
+
+- plot(functions, start, end)
+
+The functions argument can a function or list of functions. It also contains a scale and you can provide a starting and ending limit. We will add zooming later.
+
+------------------------------------------------------------------------
+
 ## 📦 Installation
 
 `pip install axiomx`
@@ -104,7 +116,8 @@ a = Matrix.from_list([]) # enter your matrix in it.
 ```
 from AxiomX.matrix import Matrix
 
-A = Matrix.from_list(\[\[1, 2\], \[3, 4\]\]) print(A.det())
+A = Matrix.from_list([1, 2], [3, 4]])
+print(A.det())
 ```
 
 ------------------------------------------------------------------------

{axiomx-0.1.2 → axiomx-0.1.3}/README.md

@@ -8,18 +8,18 @@ exploration.
 
 ## 📢 Community
 
-Join Discord for updates and feature requests. Join [here](https://discord.gg/CrEf8mDB).
+Join Discord for updates and feature requests. Join [here](https://discord.gg/yhN5eYwT).
 
 ------------------------------------------------------------------------
 
 ## 🚀 Features
 
-### 📐 Core Mathematics
+### 📐 Core Mathematics (AxiomX.constants)
 
 -   Constants (π, e, φ, √2, γ, ...)
 -   Custom implementations
 
-### 📈 Exponential & Logarithmic
+### 📈 Exponential & Logarithmic (AxiomX.exp)
 
 -   exp(x)
 -   ln(x)
@@ -27,22 +27,25 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   log2(x)
 -   log(arg, base)
 
-### 📐 Trigonometry
+### 📐 Trigonometry (AxiomX.trig)
 
 -   sin, cos, tan, cot, sec, cosec
 -   inverse functions
 
-### 🔢 Hyperbolic Functions
+### 🔢 Hyperbolic Functions (AxiomX.hyperbolic)
 
 -   sinh, cosh, tanh, coth
 
-### 🧮 General Functions
+### 🧮 General Functions (AxiomX.functions)
 
 -   absolute(x)
 -   sqrt(x)
 -   gamma(x)
 -   cbrt(x)
 -   zeta(n) and beta(n)
+-   bc(n, k) - **NEW**
+-   perm(n, K) - **NEW**
+-   pascal(row) - **NEW**
 
 ### ∫ Calculus
 
@@ -51,9 +54,7 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   convergence analysis
 -   limit evaluation
 
-------------------------------------------------------------------------
-
-## 🧱 Matrix Engine (NEW)
+### 🧱 Matrix Engine
 
 -   Matrix creation
 -   A\[i, j\] access
@@ -62,6 +63,7 @@ Join Discord for updates and feature requests. Join [here](https://discord.gg/Cr
 -   Determinant
 -   Inverse
 -   Division (via inverse) by using / operator
+-   Rank of a Matrix by Gaussian elimination (**NEW**)
 
 If you want to create a zero matrix or identity of any order, you can do it this way:
 
@@ -82,6 +84,16 @@ a = Matrix.from_list([]) # enter your matrix in it.
 ```
 ------------------------------------------------------------------------
 
+## 📈 GRAPH ENGINE - **NEW** (included in AxiomX.graph)
+
+AxiomX has introduced a killer feature - graphing. Just like how `matplotlib` is for numpy, AxiomX has introduced a graph module. It has only one function:
+
+- plot(functions, start, end)
+
+The functions argument can a function or list of functions. It also contains a scale and you can provide a starting and ending limit. We will add zooming later.
+
+------------------------------------------------------------------------
+
 ## 📦 Installation
 
 `pip install axiomx`
@@ -92,7 +104,8 @@ a = Matrix.from_list([]) # enter your matrix in it.
 ```
 from AxiomX.matrix import Matrix
 
-A = Matrix.from_list(\[\[1, 2\], \[3, 4\]\]) print(A.det())
+A = Matrix.from_list([1, 2], [3, 4]])
+print(A.det())
 ```
 
 ------------------------------------------------------------------------

{axiomx-0.1.2 → axiomx-0.1.3}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "AxiomX"
-version = "0.1.2"
+version = "0.1.3"
 authors = [
     { name = "Eric Joseph" }
 ]

{axiomx-0.1.2 → axiomx-0.1.3}/setup.py

@@ -2,7 +2,7 @@ from setuptools import setup, find_packages
 
 setup(
     name="AxiomX",  # rename to your project name
-    version="0.1.2",
+    version="0.1.3",
     description="A custom Python math library with special functions.",
     author="Eric Joseph",  # change if needed
     packages=find_packages(),

axiomx-0.1.2/AxiomX/constants.py

@@ -1,17 +0,0 @@
-pi = 3.141592653589793
-e = 2.718281828459045
-tau = 2*pi
-lemniscate = 2.622057554292119
-gauss = lemniscate / pi
-euler_mascheroni = 0.577215664901533
-sqrt_2 = 2**0.5
-sqrt_3 = 3**0.5
-sqrt_5 = 5**0.5
-golden_ratio = (1 + sqrt_5) / 2
-silver_ratio = 1 + sqrt_2
-gelfond = e**pi
-gelfond_schneider = 2**sqrt_2
-infinity = float("inf")
-
-def metallic_ratio(n):
-    return (n + (n**2 + 4)**0.5) / 2