nbmath 0.1.0__tar.gz

nbmath-0.1.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 tc0512
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

nbmath-0.1.0/PKG-INFO

@@ -0,0 +1,20 @@
+Metadata-Version: 2.4
+Name: nbmath
+Version: 0.1.0
+Summary: python零依赖数学库
+Author-email: tc0512 <tancheng_0812@qq.com>
+License: MIT
+Project-URL: Homepage, https://github.com/tc0512/nbmath
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Dynamic: license-file
+
+# nbmath
+
+一个实用的数学工具包，支持方程求解、几何计算、统计分析等功能。
+
+## 安装
+```bash
+pip install nbmath
+```

nbmath-0.1.0/README.md

@@ -0,0 +1,8 @@
+# nbmath
+
+一个实用的数学工具包，支持方程求解、几何计算、统计分析等功能。
+
+## 安装
+```bash
+pip install nbmath
+```

nbmath-0.1.0/nbmath/__init__.py

@@ -0,0 +1,16 @@
+from . import const
+from . import equation
+from . import geometry
+from . import plots
+from . import stats
+from . import utils
+
+__all__ = [
+    'const',
+    'equation', 
+    'geometry',
+    'plots',
+    'stats',
+    'utils'
+]
+__version__ = "0.1.0"

nbmath-0.1.0/nbmath/const.py

@@ -0,0 +1,20 @@
+def pi(): #π
+    return 3.1415926535897932
+def e(): #e
+    return 2.7182818284590452
+def tau(): #2π
+    return 6.2831853071795864
+def phi(): #φ
+    return 1.6180339887498948
+def G(): #万有引力常数(N·m^2/kg^2)
+    return 6.67430e-11
+def c(): #光速(m/s)
+    return 299792458
+def h(): #普朗克常数(J·s)
+    return 6.62607015e-34
+def g(): #重力加速度(N/kg)
+    return 9.80665
+def k(): #玻尔兹曼常数(J/K)
+    return 1.380649e-23
+def NA(): #阿伏伽得罗常数(mol^-1)
+    return 6.02215076e23

nbmath-0.1.0/nbmath/equation.py

@@ -0,0 +1,73 @@
+from cmath import sqrt
+I = 1j
+def solve_linear_equation_in_one_unknown(a, b): #一元一次方程
+    return float('nan') if a==0 else [-b/a]
+def solve_quadratic_equation(a, b, c): #一元二次方程
+    if a==0:
+        return float('nan')
+    x1 = (-b+sqrt(b**2-4*a*c))/(2*a)
+    x2 = (-b-sqrt(b**2-4*a*c))/(2*a)
+    return [x1, x2]
+def solve_cubic_equation_in_one_unknown(a, b, c, d):#一元三次方程
+    if a==0:
+        return float('nan')
+    return [-(-3*c/a + b**2/a**2)/(3*(sqrt(-4*(-3*c/a + b**2/a**2)**3 + (27*d/a - 9*b*c/a**2 + 2*b**3/a**3)**2)/2 + 27*d/(2*a) - 9*b*c/(2*a**2) + b**3/a**3)**(1/3)) - (sqrt(-4*(-3*c/a + b**2/a**2)**3 + (27*d/a - 9*b*c/a**2 + 2*b**3/a**3)**2)/2 + 27*d/(2*a) - 9*b*c/(2*a**2) + b**3/a**3)**(1/3)/3 - b/(3*a), -(-3*c/a + b**2/a**2)/(3*(-1/2 - sqrt(3)*I/2)*(sqrt(-4*(-3*c/a + b**2/a**2)**3 + (27*d/a - 9*b*c/a**2 + 2*b**3/a**3)**2)/2 + 27*d/(2*a) - 9*b*c/(2*a**2) + b**3/a**3)**(1/3)) - (-1/2 - sqrt(3)*I/2)*(sqrt(-4*(-3*c/a + b**2/a**2)**3 + (27*d/a - 9*b*c/a**2 + 2*b**3/a**3)**2)/2 + 27*d/(2*a) - 9*b*c/(2*a**2) + b**3/a**3)**(1/3)/3 - b/(3*a), -(-3*c/a + b**2/a**2)/(3*(-1/2 + sqrt(3)*I/2)*(sqrt(-4*(-3*c/a + b**2/a**2)**3 + (27*d/a - 9*b*c/a**2 + 2*b**3/a**3)**2)/2 + 27*d/(2*a) - 9*b*c/(2*a**2) + b**3/a**3)**(1/3)) - (-1/2 + sqrt(3)*I/2)*(sqrt(-4*(-3*c/a + b**2/a**2)**3 + (27*d/a - 9*b*c/a**2 + 2*b**3/a**3)**2)/2 + 27*d/(2*a) - 9*b*c/(2*a**2) + b**3/a**3)**(1/3)/3 - b/(3*a)]
+def solve_quartic_equation(a, b, c, d, e):
+    if a==0:
+        return float('nan')
+    x1 = -sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2))/2 - sqrt((2*d/a - b*c/a**2 + b**3/(4*a**3))/sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2)) + 2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) - 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 4*c/(3*a) + b**2/(2*a**2))/2 - b/(4*a)
+    x2 = -sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2))/2 + sqrt((2*d/a - b*c/a**2 + b**3/(4*a**3))/sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2)) + 2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) - 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 4*c/(3*a) + b**2/(2*a**2))/2 - b/(4*a)
+    x3 = sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2))/2 - sqrt(-(2*d/a - b*c/a**2 + b**3/(4*a**3))/sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2)) + 2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) - 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 4*c/(3*a) + b**2/(2*a**2))/2 - b/(4*a)
+    x4 = sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2))/2 + sqrt(-(2*d/a - b*c/a**2 + b**3/(4*a**3))/sqrt(-2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) + 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 2*c/(3*a) + b**2/(4*a**2)) + 2*(-e/a + b*d/(4*a**2) - c**2/(12*a**2))/(3*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3)) - 2*((c/a - 3*b**2/(8*a**2))**3/216 - (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/6 + sqrt((-e/a + b*d/(4*a**2) - c**2/(12*a**2))**3/27 + (-(c/a - 3*b**2/(8*a**2))**3/108 + (c/a - 3*b**2/(8*a**2))*(e/a - b*d/(4*a**2) + b**2*c/(16*a**3) - 3*b**4/(256*a**4))/3 - (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/8)**2/4) + (d/a - b*c/(2*a**2) + b**3/(8*a**3))**2/16)**(1/3) - 4*c/(3*a) + b**2/(2*a**2))/2 - b/(4*a)
+    return [x1, x2, x3, x4]
+def diff(fx: list): #求导
+    if len(fx)==0:
+        return float('nan')
+    elif len(fx)==1:
+        return 0
+    n = len(fx)
+    fx.reverse()
+    fpx = []
+    for i in range(n-1, 0, -1):
+        fpx.append(i*fx[i])
+    return fpx
+def polyval(f: list, x): #代入求值
+    if len(f)==0:
+        return float('nan')
+    elif len(f)==1:
+        return f[0]
+    fx = 0
+    power = len(f)-1
+    for i in f:
+        fx+=i*x**power
+        power-=1
+    return fx
+def newton_solver(fx: list, x0, depth):
+    x = x0
+    for i in range(depth):
+        fpx = diff(fx)
+        val_fpx = polyval(fpx, x)
+        if val_fpx==0:
+            raise ValueError("the derivative is zero so cannot iterate")
+        val_fx = polyval(fx, x)
+        x = x-val_fx/val_fpx
+    return x
+def solve(*args): #统一求解函数接口
+    if len(args)==2: #双参数->一元一次
+        a, b = args
+        return solve_linear_equation_in_one_unknown(a, b)
+    elif len(args)==3: #三参数->一元二次/牛顿法
+        if isinstance(args[0], list): #第一参数为列表->牛顿法
+            fx, x0, depth = args
+            return newton_solver(fx, x0, depth)
+        else: #否则(第一参数为数字)->一元二次
+            a, b, c = args
+            return solve_quadratic_equation(a, b, c)
+    elif len(args)==4: #四参数->一元三次
+        a, b, c, d = args
+        return solve_cubic_equation_in_one_unknown(a, b, c, d)
+    elif len(args)==5: #五参数->一元四次
+        a, b, c, d, e = args
+        return solve_quartic_equation(a, b, c, d, e)
+    else: #都不是->不支持,报错
+        raise NotImplementedError("don't support this type of parameter")

nbmath-0.1.0/nbmath/geometry.py

@@ -0,0 +1,92 @@
+import math
+class Point: #点
+    def __init__(self, x, y):
+        self.x = x
+        self.y = y
+    def expression(self): #表达式
+        return f"({self.x}, {self.y})"
+def dist(p1, p2): #两点距离
+    return math.sqrt((p2.x-p1.x)**2+(p2.y-p1.y)**2)
+class Line: #线段
+    def __init__(self, p1, p2): #p1,p2均为Point几何对象
+        self.p1 = p1
+        self.p2 = p2
+    def length(self): #长度
+        return dist(self.p1, self.p2)
+    def midpoint(self): #中点
+        return Point((self.p1.x+self.p2.x)/2, (self.p1.y+self.p2.y)/2)
+    def expression(self): #表达式
+        k = (self.p2.y-self.p1.y)/(self.p2.x-self.p1.y)
+        b = self.p1.y-k*self.p1.x
+        x_min = min(self.p1.x, self.p2.x)
+        x_max = max(self.p1.x, self.p2.x)
+        if k!=0:
+            if b>0:
+                return f"y={k}x+{b}{{{x_min}≤x≤{x_max}}}"
+            elif b==0:
+                return f"y={k}x{{{x_min}≤x≤{x_max}}}"
+            elif b<0:
+                return f"y={k}x-{-b}{{{x_min}≤x≤{x_max}}}"
+        else:
+            return f"y={b}{{{x_min}≤x≤{x_max}}}"
+class Circle: #圆
+    def __init__(self, O, r): #圆心是Point几何对象
+        self.O = O
+        self.r = r
+    def area(self): #面积
+        return math.pi*self.r**2
+    def circumference(self): #周长
+        return 2*math.pi*self.r
+    def diameter(self): #直径
+        return 2*self.r
+    def contains(self, point): #是否在圆内
+        return dist(self.O, point)<=self.r
+    def expression(self): #表达式
+        a = self.O.x
+        b = self.O.y
+        if a>0:
+            if b>0:
+                return f"(x-{a})^2+(y-{b})^2={self.r**2}"
+            elif b==0:
+                return f"(x-{a})^2+y^2={self.r**2}"
+            elif b<0:
+                return f"(x-{a})^2+(y+{-b})^2={self.r**2}"
+        elif a==0:
+            if b>0:
+                return f"x^2+(y-{b})^2={self.r**2}"
+            elif b==0:
+                return f"x^2+y^2={self.r**2}"
+            elif b<0:
+                return f"x^2+(y+{-b})^2={self.r**2}"
+        elif a<0:
+            if b>0:
+                return f"(x+{-a})^2+(y-{b})^2={self.r**2}"
+            elif b==0:
+                return f"(x+{-a})^2+y^2={self.r**2}"
+            elif b<0:
+                return f"(x+{-a})^2+(y+{-b})^2={self.r**2}"
+class Polygon:
+    def __init__(self, points): #points是Point几何对象
+        self.points = points
+    def area(self): #面积
+        n = len(self.points) 
+        if n<3:
+            return 0
+        s = 0
+        for i in range(n):
+            x1, y1 = self.points[i].x, self.points[i].y
+            x2, y2 = self.points[(i+1)%n].x, self.points[(i+1)%n].y
+            s+=x1*y2-x2*y1
+        return abs(s)/2
+    def circumference(self):
+        n = len(self.points)
+        if n<2:
+            return 0
+        elif n==2:
+            return dist(self.points[0], self.points[1])
+        C = 0
+        for i in range(n):
+            p1 = self.points[i]
+            p2 = self.points[(i+1)%n]
+            C+=dist(p1, p2)
+        return C

nbmath-0.1.0/nbmath/plots.py

@@ -0,0 +1,72 @@
+import turtle as t
+_length = _height = None
+_llx = _lly = _urx = _ury = None
+def window(length, height): #窗口(t.setup)
+    global _length, _height
+    _length, _height = length, height
+    t.setup(length, height)
+    t.hideturtle()
+    t.speed(0)
+    t.tracer(0)
+    t.penup()
+def setax(llx, lly, urx, ury): #设置坐标系范围(t.setworldcoordinates)
+    global _llx, _lly, _urx , _ury
+    _llx, _lly, _urx , _ury = llx, lly, urx, ury
+    t.setworldcoordinates(llx, lly, urx, ury)
+def drawaxhline(): #绘制坐标轴
+    global _length, _height, _llx, _lly, _urx , _ury
+    if None in (_length, _height):
+        raise ValueError("please use 'window' first to create a window")
+    elif None in (_llx, _lly, _urx , _ury):
+        x_min, x_max = -_length/2, _length/2
+        y_min, y_max = -_height/2, _height/2
+    else:
+        x_min, x_max = _llx, _urx
+        y_min, y_max = _lly, _ury
+    t.penup()
+    t.goto(x_min, 0)
+    t.pendown()
+    t.goto(x_max, 0)
+    t.penup()
+    t.goto(0, y_min)
+    t.pendown()
+    t.goto(0, y_max)
+    t.penup()
+    t.update()
+def point(x, y, color: str, size: int, label: str): #描点
+    t.penup()
+    t.goto(x, y)
+    t.pencolor(color)
+    t.dot(size)
+    t.write(label)
+    t.penup()
+    t.update()
+def line(x: list, y: list, color: str, linewidth: int): #绘制线段
+    t.penup()
+    t.goto(x[0], y[0])
+    t.pendown()
+    t.pencolor(color)
+    t.pensize(linewidth)
+    t.goto(x[1], y[1])
+    t.pensize(1)
+    t.update()
+def plot_function(f, x_min, x_max, color: str, linewidth: int, steps: int): #函数图像y=f(x)
+    t.penup()
+    dx = (x_max-x_min)/steps
+    for i in range(steps+1):
+        x = x_min+i*dx
+        y = f(x)
+        t.goto(x, y)
+        if i==0:
+            t.pendown()
+            t.pencolor(color)
+            t.pensize(linewidth)
+    t.penup()
+    t.update()
+def scatter(points: list, color: str): #散点图
+    t.penup()
+    t.pencolor(color)
+    for i in points:
+        t.goto(i)
+        t.dot(5)
+    t.update()

nbmath-0.1.0/nbmath/stats.py

@@ -0,0 +1,50 @@
+import statistics
+import math
+def mean(data: list):
+    if len(data)==0:
+        raise statistics.StatisticsError("mean requires at least one data point")
+    return sum(data)/len(data)
+def SS(data: list):
+    if len(data)==0:
+        raise statistics.StatisticsError("sum of deviation square requires at least one data point")
+    average = mean(data)
+    ss = 0
+    for i in data:
+        ss+=(i-average)**2
+    return ss
+def var(data: list):
+    if len(data)==0:
+        raise statistics.StatisticsError("variance requires at least one data point")
+    elif len(data)==1:
+        return 0
+    return SS(data)/len(data)
+def std(data: list):
+    if len(data)==0:
+        raise statistics.StatisticsError("standard deviation requires at least one data point")
+    elif len(data)==1:
+        return 0
+    return math.sqrt(var(data))
+def data_range(data: list):
+    if len(data)<2:
+        raise statistics.StatisticsError("range requires at least two data points")
+    return max(data)-min(data)
+def percentile(data: list, p: int):
+    if len(data)==0:
+        raise statistics.StatisticsError("percentile requires at least one data point")
+    elif len(data)==1:
+        return data[0]
+    data.sort()
+    pos = len(data)*p/100
+    if int(pos)==pos:
+        index1 = int(pos-1)
+        index2 = int(pos)
+        return (data[index1]+data[index2])/2
+    else:
+        index = int(pos)
+        return data[index]
+def mode(data):
+    freq = {}
+    for x in data:
+        freq[x] = freq.get(x, 0) + 1
+    max_freq = max(freq.values())
+    return [k for k, v in freq.items() if v == max_freq]

nbmath-0.1.0/nbmath/utils.py

@@ -0,0 +1,68 @@
+import time
+def diff(fx: list): #求导
+    if len(fx)==0:
+        return float('nan')
+    elif len(fx)==1:
+        return 0
+    n = len(fx)
+    fx.reverse()
+    fpx = []
+    for i in range(n-1, 0, -1):
+        fpx.append(i*fx[i])
+    return fpx
+def polyval(f: list, x): #代入求值
+    if len(f)==0:
+        return float('nan')
+    elif len(f)==1:
+        return f[0]
+    fx = 0
+    power = len(f)-1
+    for i in f:
+        fx+=i*x**power
+        power-=1
+    return fx
+def linspace(start, end, steps): #np.linspace的纯python实现
+    if steps==0:
+        return []
+    elif steps==1:
+        return [start]
+    dx = (end-start)/steps
+    result = []
+    for i in range(steps+1):
+        x = start+i*dx
+        result.append(x)
+    return result
+def timer(func): #计时器
+    def wrapper(*args, **kwargs):
+        start = time.time()
+        result = func(*args, **kwargs)
+        print(f"{func.__name__} 耗时: {time.time()-start:.4f}秒")
+        return result
+    return wrapper
+def gcd(a: int, b: int): #最大公约数
+    while b:
+        a, b = b, a%b
+    return abs(a)
+def lcm(a: int, b: int): #最小公倍数
+    return a*b/gcd(a, b)
+def floor(x): #向下取整
+    if x>=0:
+        return int(x)
+    return int(x)-1
+def trunc(x): #截断取整
+    return int(x)
+def frac(x): #小数部分
+    if x>=0:
+        return x-int(x)
+    return x-int(x)+1
+def fac(x): #阶乘
+    if x>0 and int(x)==x:
+        result = 1
+        for i in range(1, x+1):
+            result*=i
+        return result
+    elif x==0:
+        return 1
+    elif x<0 and int(x)==x:
+        return float('nan')
+    return math.gamma(x+1)

nbmath-0.1.0/nbmath.egg-info/PKG-INFO

@@ -0,0 +1,20 @@
+Metadata-Version: 2.4
+Name: nbmath
+Version: 0.1.0
+Summary: python零依赖数学库
+Author-email: tc0512 <tancheng_0812@qq.com>
+License: MIT
+Project-URL: Homepage, https://github.com/tc0512/nbmath
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Dynamic: license-file
+
+# nbmath
+
+一个实用的数学工具包，支持方程求解、几何计算、统计分析等功能。
+
+## 安装
+```bash
+pip install nbmath
+```

nbmath-0.1.0/nbmath.egg-info/SOURCES.txt

@@ -0,0 +1,14 @@
+LICENSE
+README.md
+pyproject.toml
+nbmath/__init__.py
+nbmath/const.py
+nbmath/equation.py
+nbmath/geometry.py
+nbmath/plots.py
+nbmath/stats.py
+nbmath/utils.py
+nbmath.egg-info/PKG-INFO
+nbmath.egg-info/SOURCES.txt
+nbmath.egg-info/dependency_links.txt
+nbmath.egg-info/top_level.txt

nbmath-0.1.0/nbmath.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

nbmath-0.1.0/nbmath.egg-info/top_level.txt

@@ -0,0 +1 @@
+nbmath

nbmath-0.1.0/pyproject.toml

@@ -0,0 +1,18 @@
+[build-system]
+requires = ["setuptools>=61.0"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "nbmath"
+version = "0.1.0"
+description = "python零依赖数学库"
+authors = [
+    {name = "tc0512", email = "tancheng_0812@qq.com"}
+]
+license = {text = "MIT"}
+readme = "README.md"
+requires-python = ">=3.8"
+dependencies = []
+
+[project.urls]
+Homepage = "https://github.com/tc0512/nbmath"

nbmath-0.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+