python微积分近似值_《用 Python 学微积分》笔记 1

《用 Python 学微积分》原文见参考资料 1。

1、多项式

f(x)=x3-5x2+9

deff(x):return x**3 - 5*x**2 + 9

print f(3)print f(1)importnumpy as np

x= np.linspace(-5, 5, num = 100)

y=f(x)importmatplotlib.pyplot as plt

plt.plot(x,y)

plt.show()

View Code

2、指数函数

exp(x)=ex

importnumpy as npimportmatplotlib.pyplot as pltdefexp(x):return np.e**xprint exp(2)print np.exp(2)

x= np.linspace(-5, 5, num = 100)

y=exp(x)

plt.plot(x,y)

plt.show()

View Code

3、对数函数

importnumpy as npimportmatplotlib.pyplot as plt

x= np.linspace(0.1,10,99,endpoint =False)

y1=np.log2(x)

y2=np.log(x)

y3=np.log10(x)

plt.plot(x,y1,'red',x,y2,'yellow',x,y3,'blue')

plt.show()

View Code

4、三角函数

sin(x)

importnumpy as npimportmatplotlib.pyplot as plt

plt.plot(np.linspace(-2*np.pi,2*np.pi),np.sin(np.linspace(-2*np.pi,2*np.pi)))

plt.show()

View Code

5、函数的复合

h(x)=x2+1

importnumpy as npimportmatplotlib.pyplot as pltdef f(x): return x+1

def g(x): return x**2

def h(x): returnf(g(x))

x= np.array(range(-10,10))

y= np.array([h(i) for i inx])

plt.plot(x, y,'bo')#h2 = lambda x: f(g(x))#plt.plot(x,h2(x),'rs')

plt.show()

View Code

6、逆函数

w=x2

winv(x)=x1/2