24小时热门版块排行榜

返回列表

zzlovejyt

金虫 (著名写手)

应助: 7 (幼儿园)
金币: 1287.8
散金: 66
红花: 3
帖子: 1048
在线: 67.4小时
虫号: 1986630
注册: 2012-09-09
性别: GG
专业: 光学

[求助] 将C语言编写的蒙特卡罗方法求解圆周率改写为fortran语言

说明：本例中先取0到20000之间的随机数，然后再用这些点除以10000即可得到0到2之间的随机数，随机数精度为五位有效数。
以下是用c++语言模拟的求解圆周率的程序：

#include
#include //setw函数,setprecision函数
#include    //time函数
#include //srand,rand函数
using namespace std;
int main()
{
int dwCount;
cout<<"请输入打点的数目"< cin>>dwCount;          //自己输入所要取的点数。
double j,k;
double x,y;
double m=0;
for(int i=0;i {
j=(int)(20000.0*rand()/(RAND_MAX+0));//取从0到20000之间的随机数
k=(int)(20000.0*rand()/(RAND_MAX+0));//取从0到20000之间的随机数
      x=j/10000; /*用刚才所取得随机数除以10000就是0到2之间的随机数，5位有效数字，代表横坐标*/
y=k/10000; /*用刚才所取得随机数除以10000就是0到2之间的随机数，5位有效数字，代表纵坐标*/
if(((x-1)*(x-1)+(y-1)*(y-1))<=1) //判断所取得点是否在圆内
{
         m++; //如果这个所取得点在圆内，则计数器加1
}
cout< }
cout< double s=m/dwCount; //用圆内的点数除以所取得点数,得到一个比例
cout< double t=4*s; /*用刚才的点数乘以4就是圆的面积，因为这里做了特殊，所以这个值也就是圆周率*/
cout< return 0;
}

回复此楼

» 猜你喜欢

面上没中，邀请各位路过的虫友分析一下分数已经有13人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有8人回复
两块石头已经有8人回复
初秋的晨风已经有6人回复
学科评审组评审是指会评吗？已经有9人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有14人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有7人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有17人回复
Ei源刊怎么投已经有4人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复

人生苦短，何妨一试,每个读硕读博的人都曾在一段苦逼的岁月里像个二逼的坚持着!

1楼 2012-11-29 17:49:57

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zzlovejyt 的主题更新

返回列表