【数值分析】劈因子法（贝尔斯托法）求解-Python实现

本篇为数值分析课程代码实现-劈因子法的实现

仅供参考

配置环境

• Python3.6
• matplotlib
• numpy

• # -*- coding: utf-8 -*-
  """
  Created on Sun Nov 25 12:56:58 2018
  劈因子法（贝尔斯托法）
  @author: hhuaf
  """
  import numpy as np
  import matplotlib.pyplot as plt
  
  #input_u=input('请输入两个数，第一个')
  #iuput_v=input('请输入两个数，第二个')
  input_u=0
  iuput_v=0
  
  theta=1e-5
  num=1000
  
  #可以显示中文
  plt.rcParams["font.sans-serif"] = ["SimHei"]
  plt.rcParams['axes.unicode_minus'] = False
  
  # 设置风格
  plt.style.use('ggplot')
  
  
  # 因子式
  def f(x):
      return (x+1)*(x-4)*(x-5)*(x+3)*(x-2)
     
  # 一般式
  def f1(x):
      return x**5-7*x**4-3*x**3+79*x**2-46*x-120
  # 原函数系数
  
  coefficient=[1,-7,-3,79,-46,-120]
  root_list=[]
  
  def two(a):
      return -a[1]/a[0]
  
  def Bairstow(a = coefficient, u = 0, v = 0, num = 1000, theta = 1