python积分计算高等数学_奇思妙想（1）python解积分高数题

从今天开始把一些电脑里积攒的当年奇思妙想拿出来，大家看个乐子

今天做高炉碳氧比监测模型的时候，有一个用到积分的数学模型，然后我就想到了当时专门写过一个python解积分的笔记。

首先这个数学模型是这样的

一个掺杂了很多化学知识的监测模型。实现也很简单。

from sympy import *

def compute_Wc(time,W_st,Q_gas,C_0):

C_de=integrate(t, (t, 0, time))

Wc=0.1*(C_0-C_de)/W_st

'''这个问题一氧化碳和二氧化碳摩尔量随时间变化这个信息没给

但是考虑到实际监测，我跟老师提出了这几种想法：

①使用递推函数，每个监测时间节点的函数值参考上一个节点做计算

②将每次监测数据结果以数组的形式储存，但是会耗费大量空间

③寻求CO和CO2量和其他参数的函数关系

'''

然后我就想到了scipy包，随后对当时的笔记和积分函数做了学习。

import scipy.integrate as integrate

import scipy.special as special

result = integrate.quad(lambda x: 2*x, 9, 10)

print(result)

#这是最基本的quad函数，quad函数第一个参数应当是一个python对象或者方法，lambda匿名函数是最合适的

#当然也可以用单个参数的def函数定义，返回值必须是和参数有关的表达式&#