Python计算开方、立方、圆周率,精确到小数点后任意位的方法


Posted in Python onJuly 17, 2018

Python计算的位数

在电脑上做了一个实验,看看python能计算到多少位,一下是结果。

x = math.sqrt((3))
print ("%.53f"%(x))
print ("%.63f"%(x))
print ("%.83f"%(x))
1.73205080756887719317660412343684583902359008789062500
1.732050807568877193176604123436845839023590087890625000000000000
1.73205080756887719317660412343684583902359008789062500000000000000000000000000000000
1.73205080756887719317660412343684583902359008789062500000000000000000000000000000000

计算立方根公式

设A = X^3,求X.称为开立方。 开立方有一个标准的公式:

Python计算开方、立方、圆周率,精确到小数点后任意位的方法

例如,A=5,,即求5介于1的3次方至2的3次方之间(1的3次方=1,2的3次方=8)

初始值X0可以取1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,都可以。例如我们取X0 = 1.9按照公式:

第一步:X1=1.9+(5/1.9²-1.9)1/3=1.7;

即5/1.9×1.9=1.3850416,1.3850416-1.9=-0.5149584,-0.5149584×1/3=-0.1716528,1.9+(-0.1716528)=1.7。即取2位数值,,即1.7。

第二步:X2=1.7+(5/1.7²-1.7)1/3=1.71;

即5/1.7×1.7=1.73010,1.73-1.7=0.03,0.03×1/3=0.01,1.7+0.01=1.71。取3位数,比前面多取一位数。

第三步:X3=1.71+(5/1.71²-1.71)1/3=1.709;

第四步:X4=1.709+(5/1.709²-1.709)1/3=1.7099;

这种方法可以自动调节,第一步与第三步取值偏大,但是计算出来以后输出值会自动转小;第二步,第四步输入值偏小,输出值自动转大。即5=1.7099³ 当然初始值X0也可以取1.1,1.2,1.3,。。。1.8,1.9中的任何一个,都是X1 = 1.7 。当然,我们在实际中初始值最好采用中间值,即1.5。 1.5+(5/1.5²-1.5)1/3=1.7。

增加计算步骤,可以得到任意精度的值

计算根号2任意位数

开平方公式

如果用这个公式开平方,只需将3改成2,2改成1。即

Python计算开方、立方、圆周率,精确到小数点后任意位的方法

import sys
n = 2
s = 0
while True:
 for si in range(9,-1,-1):
 nx = n - ((2*s*10+si)*si)
 if nx>=0:
  s = s*10+si
  n = nx*100
  sys.stdout.write(str(si))
  sys.stdout.flush()
  break

计算开5次方公式,精确至任意位数

Python计算开方、立方、圆周率,精确到小数点后任意位的方法

计算圆周率任意位数

因为根据马青公式π/4=4arctg1/5-arctg1/239

又因为arctgX=X-(1/3)X^3+(1/5)X^5-(1/7)X^7+……+[(-1)^(n-1)/((2n-1)]*X^(2n-1))

变形得π/4=(4/5-1/239)-1/3(4/5^3-1/239^3)+1/5(4/5^5-1/239^5)……

所以可以用python语言编写出求圆周率到任意位的程序如下:

n = int(raw_input('请键入想要计算到小数点后的位数n:')) #先键入字符串,再转化为整数

w = n+10 #多计算10位,防止尾数取舍的影响

b = 10**w #算到小数点后w位

x1 = b*4//5 #求含4/5的首项

x2 = b// -239 #求含1/239的首项

he = x1+x2 #求第一大项

n *= 2 #设置下面循环的终点,即共计算n项

for i in xrange(3,n,2): #循环初值=3,末值2n,步长=2

x1 //= -25 #求每个含1/5的项及符号

x2 //= -57121 #求每个含1/239的项及符号

x = (x1+x2) // i #求两项之和

he += x #求总和

pai = he*4 #求出π

pai //= 10**10 #舍掉后十位

print pai #输出圆周率π的值,

以上这篇Python计算开方、立方、圆周率,精确到小数点后任意位的方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Python进行数据提取的方法总结
Aug 22 Python
python机器学习实战之最近邻kNN分类器
Dec 20 Python
Django中STATIC_ROOT和STATIC_URL及STATICFILES_DIRS浅析
May 08 Python
Python 从相对路径下import的方法
Dec 04 Python
python实现画五角星和螺旋线的示例
Jan 20 Python
解决python3中的requests解析中文页面出现乱码问题
Apr 19 Python
python+openCV调用摄像头拍摄和处理图片的实现
Aug 06 Python
pytorch 在sequential中使用view来reshape的例子
Aug 20 Python
Pytorch Tensor的统计属性实例讲解
Dec 30 Python
python 使用多线程创建一个Buffer缓存器的实现思路
Jul 02 Python
python解决12306登录验证码的实现
Apr 18 Python
用Python爬取某乎手机APP数据
Jun 15 Python
python把数组中的数字每行打印3个并保存在文档中的方法
Jul 17 #Python
python3.x 将byte转成字符串的方法
Jul 17 #Python
Python 字符串转换为整形和浮点类型的方法
Jul 17 #Python
Python操作word常见方法示例【win32com与docx模块】
Jul 17 #Python
python查看列的唯一值方法
Jul 17 #Python
Python使用win32com模块实现数据库表结构自动生成word表格的方法
Jul 17 #Python
python 数字类型和字符串类型的相互转换实例
Jul 17 #Python
You might like
php分页查询mysql结果的base64处理方法示例
2017/05/18 PHP
PHP7下协程的实现方法详解
2017/12/17 PHP
php实现文件上传基本验证
2020/03/04 PHP
基于jQuery的js分页代码
2010/06/10 Javascript
return false;和e.preventDefault();的区别
2010/07/11 Javascript
javascript自然分类法算法实现代码
2013/10/11 Javascript
bootstrap网页框架的使用方法
2016/05/10 Javascript
JavaScript语言精粹经典实例(整理篇)
2016/06/07 Javascript
纯js实现图片匀速淡入淡出效果
2017/08/22 Javascript
Bootstrap 3多级下拉菜单实例
2017/11/23 Javascript
node vue项目开发之前后端分离实战记录
2017/12/13 Javascript
angular ng-model 无法获取值的处理方法
2018/10/02 Javascript
微信小程序文章详情页跳转案例详解
2019/07/09 Javascript
[01:13]2015国际邀请赛线下观战现场
2015/08/08 DOTA
python抓取网页内容示例分享
2014/02/24 Python
python 队列详解及实例代码
2016/10/18 Python
python模块之sys模块和序列化模块(实例讲解)
2017/09/13 Python
Python简单实现控制电脑的方法
2018/01/22 Python
python装饰器深入学习
2018/04/06 Python
flask中的wtforms使用方法
2018/07/21 Python
python进阶之多线程对同一个全局变量的处理方法
2018/11/09 Python
matlab中imadjust函数的作用及应用举例
2020/02/27 Python
Pandas将列表(List)转换为数据框(Dataframe)
2020/04/24 Python
使用HTML5中的contentEditable来将多行文本自动增高
2016/03/01 HTML / CSS
Html5调用手机摄像头并实现人脸识别的实现
2018/12/21 HTML / CSS
美国购买舞会礼服网站:Couture Candy
2019/12/29 全球购物
党委书记岗位职责
2013/11/24 职场文书
咖啡厅创业计划书范本
2014/01/22 职场文书
文秘应届生求职信
2014/07/05 职场文书
大学生找工作求职信
2014/07/09 职场文书
2014年教师思想工作总结
2014/12/03 职场文书
秋季运动会加油词
2015/07/18 职场文书
离婚民事起诉状
2015/08/03 职场文书
家长对孩子的寒假评语
2015/10/09 职场文书
2015年党务工作者个人工作总结
2015/10/22 职场文书
OpenCV-Python实现怀旧滤镜与连环画滤镜
2021/06/09 Python