【数值分析】试位法求解-Python实现

本篇为数值分析课程代码实现-试位法的实现

仅供参考

配置环境

• Python3.6
• matplotlib
• numpy

• # -*- coding: utf-8 -*-
  """
  Created on Fri Nov 16 21:33:10 2018
  试位法法
  @author: hhuaf
  """
  import numpy as np
  import matplotlib.pyplot as plt
  '''
  a :下限
  b :上限
  theta:阈值
  '''
  #a=input('输入下限\n:')
  #b=input('输入上限\n:')
  #a=float(a)
  #b=float(b)
  a=-1
  b=1
  theta=0.05
  
  #可以显示中文
  plt.rcParams["font.sans-serif"] = ["SimHei"]
  plt.rcParams['axes.unicode_minus'] = False
  
  # 设置风格
  plt.style.use('ggplot')
  
  # 定义函数，构造数值
  fun = lambda x: x**3-3*x+1
  #fun = lambda x: x**2-1
  x = np.arange(a,b,0.05)
  y = fun(x)
  
  # 画出图像
  
  fig = plt.figure(figsize = (15, 10))
  plt.subplot(2,3,1)
  plt.xlim(a, b)
  plt.xlabel('X')
  plt.ylabel('Y')
  plt.title('$f(x)=x^3-3x+1$ 图像')
  plt.plot(