1. 程式人生 > >python實現一個簡單泰勒的計算

python實現一個簡單泰勒的計算

也沒有從泰勒公式開始,就根據具體函式寫了一個簡單計算

#coding=utf-8
from sympy import *
import math

#定義變數為x
x=Symbol("x")
#函式為
f = -0.1*x**4-0.15*x**3-0.5*x**2-0.25*x+1.2
#求出一到四階導數分別為
f1 = diff(f,x,1)
f2 = diff(f,x,2)
f3 = diff(f,x,3)
f4 = diff(f,x,4)
print f1
print f2
print f3
print f4

#傳入x=0求出各階導函式的具體數值
g0 = f.subs(x,0)
g1 = f1.subs(x,0
) g2 = f2.subs(x,0) g3 = f3.subs(x,0) g4 = f4.subs(x,0) #print g0,g1,g2,g3,g4 #當x改變時傳入x的值 x=int(raw_input("請輸入改變後x的值:")) #將x與導函式的值帶入泰勒公式中求出結果 g = g0/math.factorial(0) * x**0 + g1/math.factorial(1) * x**1 + g2/math.factorial(2) * x**2 + g3/math.factorial(3) * x**3 + g4/math.factorial(4) * x**4 print g