利用拉马努金公式和蒙特卡洛方法计算圆周率
拉马努金公式
#求阶乘函数
def myreduce(m):value = 1for i in range(1,m+1):value = value*ireturn value#求和项的函数
def sumfun(m):div = (396**(4*m))*((myreduce(m))**4)#公式的分母部分return (26390*m + 1103)*myreduce(4*m)/divCOUNT = 1000#求和项的个数
sumvalue = 0#累计的值
#进行求和
for m in range(0, COUNT+1):sumvalue += sumfun(m)
temp = 2*(2**0.5)*sumvalue/(99**2)#1/piprint('%.4f'%(1/temp))
蒙特卡洛方法
from random import randomCOUNT = 1000000#总实验次数
WIDTH = 1#正方形的边长
num_in_cir = 0#在1/4圆内的实验次数
for i in range(1, COUNT+1):x = random()y = random()if (x**2 + y**2) <= WIDTH:#判断实验点是否落入圆内num_in_cir += 1
pi = 4*(num_in_cir/COUNT)#计算pi
print('%.4f'%(pi))
利用拉马努金公式和蒙特卡洛方法计算圆周率相关推荐
- 使用Python通过拉马努金公式快速求π
使用Python通过拉马努金公式快速求π 一.前言 π是一个数学常数,定义为:圆的周长与直径的比值. π是一个无理数,也是一个超越数,它的小数部分无限不循环. π可以用来精确计算圆周长.圆面积.球体积 ...
- 利用蒙特卡洛方法计算圆周率并将撒点分布和π估计值可视化(Python大作业)
这个大作业做的比较简单,其他同学动辄两三百行,我只是出于纯个人兴趣,想了解一下第三方库matplotlib的用法以及更加直观的了解蒙特卡洛方法. 首先是蒙特卡洛方法求圆周率的部分: >>& ...
- 蒙特卡洛方法计算圆周率
import random import mathdef calpai():n = 10000000 #随机实验次数r = 1.0 #圆的半径a, b = (0.0, 0.0) #中心点x_neg, ...
- python使用蒙特卡洛方法计算圆周率的流程图怎么画_在python中用蒙特卡洛算法计算圆周率...
本文写给那些python初学者与对蒙特卡洛算法感兴趣,但却不知该如何理解或应用的人. (虽然我发现这个貌似有许多人做过了,但是程序都相对复杂,不便于理解,于是我就自己编写了一段程序,海龟的可视化请看下 ...
- 拉马努金公式计算圆周率 Python 3.7
极快的收敛速度,令人惊叹! 转载请注明出处,联系我: t39q@163.com 本人热衷于数据库技术及算法的研究,志同道合之士, 欢迎探讨 import math x=0 for k in range ...
- 从拉马努金到张益唐——数学是一个整体
分享一下我老师大神的人工智能教程!零基础,通俗易懂!http://blog.csdn.net/jiangjunshow 也欢迎大家转载本篇文章.分享知识,造福人民,实现我们中华民族伟大复兴! 分享到: ...
- 传奇数学家拉马努金留下的数学神谕,解开了多年悬而未决的神秘难题
◆ ◆ ◆ 序 印度数学奇才斯里尼瓦瑟· 拉马努金在短短的一生中写下了大量论文,虽然多数文章在生前未得以发表,但借助其中一些文章中的全新见解,数学家肯恩· 小野解开了长久以来悬而未决的神秘数学难题. ...
- 拉马努金:“与神对话”的数学天才
不管你 一个无法表达神的思想的方程, 对我而言毫无意义 -- 斯里尼瓦瑟·拉马努金 本文摘自人民邮电出版社图灵文化出版的图书<悠扬的素数:二百年数学绝唱黎曼假设> 当哈代和利特尔伍德步履维 ...
- 拉马努金的圆周率计算公式
拉马努金的几个神仙公式.据说他的论文里一共发表了14条圆周率的计算公式,但目前只收集到十个. 1π=18∑m=0∞(20m+3)(−1)m(4m)!(42)4m(m!)41π=123∑m=0∞(8m+ ...
最新文章
- Java调用 shell脚本阻塞
- Calendar日历小程序
- WSL安装xfce4
- 基础知识回顾——异常处理
- 亲试白天使:华硕家用级无线路由RT-N11+
- python3 正则 去除 html标签、提取正文内容_Python通过正则表达式去除(过滤)HTML标签,提取文字...
- 20个MySQL运维案例,请查收!
- bootstrap学习笔记(2)表单设计
- 关于反射中.getDeclaredContructor()返回构造方法顺序的问题(转)
- 重庆首个程序员社区成立
- sketchup边线设置_草图大师SketchUp改变模型边线颜色的方法
- 为什么学习计算机的时候学语文,小学语文“为什么要学拼音呢”?告诉你学拼音技巧!...
- Android开发-窗口跳转
- CodeForces 372A Counting Kangaroos is Fun 动物PK
- Openstack 使用Evacuate来实现HA
- Control.DataBinding数据绑定细解
- python-InputWhile基础知识
- Mininet系列实验(三):Mininet命令延伸实验扩展
- 代码解读一 文件名“ANO_Imu.c”
- 哈工大操作系统学习笔记十——信号量与死锁
热门文章
- 福昕阅读器无法注释、高亮、打字等情况
- php图片搜索引擎,Everypixel:基于人工智能的图片搜索引擎
- Centos7安装minio文件服务器
- 需要控制好情绪,调整好心态
- 如何用OpenSSL从https网站上导出SSL的CA证书?
- htc+m7+android系统+电源,HTC One M7手机怎么刷机 htc手机刷机教程【步骤】
- cass10.1+鸿业生成平纵横数据
- 机器人学——2.2-分段轨迹
- 计算机网络原理04741教学课件,【精编】04741计算机网络原理自考本科2018年4月课件.ppt...
- OkHttp基础概念解释