应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 新手请教最大值输出问题
91/1返回列表
查看: 1282  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

19870617

木虫 (著名写手)

[求助] 新手请教最大值输出问题

自己编的很小的一段程序

  
program main
    implicit none
    integer :: n,i
    real(8) ::pi,u1,u2,r   
   pi=3.1415926535
    call random_seed ()
    n=1000      
   do i=1,n
    call random_number(u1)
    call random_number(u2)
    r=sqrt(-log(u1))*cos(2*pi*u2)  
    print *,r
    enddo
    end
现在我想把r的值存储到一个数组当中，以后每次调用数组中的一个r值，并输出数组中随机生成的r的最大值，请问大侠该如何编写，刚刚学习FORTRAN，希望大家多多指教，不胜感激。
金币少了点，希望大家多多帮忙，谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-27 17:23:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢应助 2012-08-28 00:16:19
把 r 定义为一个数组就可以了……

找本入门的书看一下吧
赞一下 回复此楼
2楼2012-08-27 17:56:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

19870617

木虫 (著名写手)
引用回帖:
2楼: Originally posted by snoopyzhao at 2012-08-27 17:56:59
把 r 定义为一个数组就可以了……

找本入门的书看一下吧

r定义为数组？一维大小为1000嘛？试了一下输出的r有3列，感觉不对啊
赞一下 回复此楼
3楼2012-08-27 18:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 19870617 at 2012-08-27 18:30:01
r定义为数组？一维大小为1000嘛？试了一下输出的r有3列，感觉不对啊...

你是怎么输出的？

一个一维数组为啥会输出三列（二维）？
赞一下 回复此楼
4楼2012-08-27 19:01:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

19870617

木虫 (著名写手)
引用回帖:
4楼: Originally posted by snoopyzhao at 2012-08-27 19:01:09
你是怎么输出的？

一个一维数组为啥会输出三列（二维）？...

print *,r 和write(*,*) r这样都试过了，都是3列，我也觉得奇怪呀
赞一下 回复此楼
5楼2012-08-27 20:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by 19870617 at 2012-08-27 20:00:41
print *,r 和write(*,*) r这样都试过了，都是3列，我也觉得奇怪呀...

把完整的代码贴出来吧……
赞一下 回复此楼
6楼2012-08-27 20:07:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

19870617

木虫 (著名写手)
引用回帖:
6楼: Originally posted by snoopyzhao at 2012-08-27 20:07:45
把完整的代码贴出来吧……...

您好！源码是这样的
        program main
     implicit none
     integer :: n,i
     real(8) ::pi,u1,u2,maxr
         real(8) :: r(1)  
    pi=3.1415926535
     call random_seed ()
     n=1000      
    do i=1,n
     call random_number(u1)
     call random_number(u2)
     r=sqrt(-log(u1))*cos(2*pi*u2)  
     print *,r
     enddo
         maxr=maxval(r)
         print *,"the largest is",maxr
    end

输出现在是1列了，但是求最大值的时候输出的maxr是最后一个r值，而不是最大值，请问maxval这样用对嘛？
赞一下 回复此楼
7楼2012-08-27 21:54:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
19870617: 金币+5, ★★★★★最佳答案, 谢谢！ 2012-08-28 08:42:26
引用回帖:
7楼: Originally posted by 19870617 at 2012-08-27 21:54:47
您好！源码是这样的
        program main
     implicit none
     integer :: n,i
     real(8) ::pi,u1,u2,maxr
         real(8) :: r(1)  
    pi=3.1415926535
     call random_seed ()
     n=1000      
    d ...

你还是去找本书看一下吧……
CODE:
program main
   implicit none
   integer :: n,i
   real(8) ::pi,u1,u2,maxr
   real(8) :: r(1000)  
   pi=3.1415926535
   call random_seed ()
   n=1000      
   do i=1,n
      call random_number(u1)
      call random_number(u2)
      r(i)=sqrt(-log(u1))*cos(2*pi*u2)  
      print *, r(i)
   enddo
   maxr=maxval(r)
   print *,"the largest is",maxr
end

赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

19870617
8楼2012-08-27 22:10:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

19870617

木虫 (著名写手)
送鲜花一朵
引用回帖:
8楼: Originally posted by snoopyzhao at 2012-08-27 22:10:01
你还是去找本书看一下吧……


program main
   implicit none
   integer :: n,i
   real(8) ::pi,u1,u2,maxr
   real(8) :: r(1000)  
   pi=3.1415926535
   call random_seed ()
   n=1000      
  ...

嗯，谢谢！还是要静下心来慢慢学！
赞一下 回复此楼
9楼2012-08-27 23:56:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 19870617 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 4/200 2026-03-16 21:52 by Losir
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
信息提示
关闭
请填处理意见
关闭