厄巧克力-CSDN博客

import sympyx=[]y=[]W=[]a=[]aa=[]beita=[]n=int(input("请输入最高项的次数:"))xxx=(input("请输入x的值，以空格分隔:"))yyy=(input("请输入y的值，以空格分隔:"))'''www=(input("请输入w的值，以空格分隔:"))'''x=xxx.split(" ")m=len(x)y=yyy.split(" ")'''W=www.split(" ")'''x=[float(x) for x in x

2021-05-06 20:27:36 206

原创 python龙贝格算法

龙贝格算法import numpy as npimport sympyx=sympy.symbols("x")yuan=input("请输入原函数：")quyu=input("请输入积分区间（用空格隔开）：")wucha=float(input("请输入误差："))Q=[]Q=quyu.split(" ")Q=[float(x) for x in Q]a=Q[0]b=Q[1]fa=sympy.sympify(yuan).subs(x,a)fb=sympy.sympify(yuan)

2021-05-06 20:24:39 790

原创利用python 解对称正定矩阵方程组的平方根法

大二在读，代码水平有限，轻喷import matha=[]L=[]U=[]ling=0jieshu=int(input("矩阵的阶数："))print("请输入方程的系数矩阵")for i in range(jieshu): onerow=[] Lonerow=[] Uonerow=[] for j in range(jieshu): x=float(input()) onerow.append(x) Loner

2021-04-15 23:45:01 1073

原创利用python解决三次样条插值

大二在读，代码水平有限只能保证第一第二边界条件正确import sympydef catchup(a,b,c,f,jieshu): beita=[] arf=[] yyy=[] m=[] for i in range(jieshu-1): beita.append(0) for i in range(jieshu): arf.append(0) for i in range(jieshu): y

2021-04-15 23:42:43 1730 2

原创利用python实现利用追赶法解决线性方程组

追赶法大二在读，能力有限，轻喷python实现追赶法解方程组，通篇都是低级函数，适合初学者学习jieshu=int(input("矩阵的阶数："))print("请输入a")a=[]for i in range(jieshu-1): x=float(input()) a.append(x)print("请输入b")b=[]for i in range(jieshu): x=float(input()) b.append

2021-04-09 17:30:23 2789

原创利用python实现杜利特尔分解法

杜利特尔分解法博主目前大二在读，代码水平有限轻喷使用的是python语言，适于数学者，通篇没有用到高级函数a=[]L=[]U=[]ling=0jieshu=int(input("矩阵的阶数："))print("请输入方程的系数矩阵")for i in range(jieshu): onerow=[] Lonerow=[] Uonerow=[] for j in range(jieshu): x=float(input()) o

2021-04-09 17:27:23 1770 5

weixin_49530111的博客

原创曲线拟合的最小二乘法：最小二乘法 python

原创曲线拟合的最小二乘法：正交二乘法 python

原创曲线拟合的最小二乘法：加权最小二乘法 python

原创非线性方程的数值解法：正割法 python

原创非线性方程的数值解法：牛顿下山法 python

原创非线性方程的数值解法：牛顿法 python

原创非线性方程的数值解法：二分法 python

原创非线性方程的数值解法：迭代法 python

原创非线性方程的数值解法：Aitken加速 python

原创 python 经典R-K(龙格库塔)

原创 python改进欧拉法

原创 python正交二乘法