princeado

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 2729.5
紅花: 5
帖子: 874
在線: 868.6小時
蟲號: 1031978
注冊: 2010-05-30
性別: GG
專業(yè): 水文、水資源

[求助] 大家?guī)臀铱匆粋€fortran的程序，我總是計算不出正確的結(jié)果

最近剛學(xué)fortran，開始編小小程序試手，下面是我弄的一個程序，沒有錯誤，但總是運(yùn)行不出正確的結(jié)果，程序是想通過do循環(huán)對數(shù)組a進(jìn)行賦值，用的是語句a(i)=jiecheng(i)，jiecheng是我自定義的函數(shù)，求每個數(shù)對應(yīng)的階乘。所有程序如下，哪位高手幫我看看問題在哪里。
program test01
implicit none
integer :: i
real :: a(5)
real,external :: jiecheng
   do i=1,5
   a(i)=jiecheng(i)
   end do
      write(*,*) a
         stop
         end
   function jiecheng(x)
   implicit none
   integer x
         real jiecheng
         integer :: ii=1
         integer :: m=1
   do ii=1,x
         m=m*ii
         end do
   jiecheng=m
                                                return
   stop
   end

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

配置一臺1萬左右的計算機(jī)跑Fortran程序以及高斯的計算已經(jīng)有8人回復(fù)
求高人幫忙看一下我這個Fortran程序的錯誤出在哪里？已經(jīng)有12人回復(fù)
我編的Simpson積分法fortran程序給不出結(jié)果，大俠們看看哪里出了問題？已經(jīng)有4人回復(fù)
請教一個fortran小程序編譯出錯的問題，謝謝已經(jīng)有9人回復(fù)
請問，我想用fortran計算統(tǒng)計中的p值以及95%信度空間已經(jīng)有4人回復(fù)
寫了一個fortran90的小程序，編譯通不過，請大俠幫忙已經(jīng)有59人回復(fù)
【求助】如何寫FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離已經(jīng)有34人回復(fù)
【求助】同樣一個問題用matlab和fortran編程計算，結(jié)果相差比較大，不知道為什么已經(jīng)有4人回復(fù)
【求助】請教如何在超算上運(yùn)行fortran程序已經(jīng)有4人回復(fù)
【求助】用fortran怎么描述下面這樣的情況，我是新手，請大家?guī)兔?/a> 已經(jīng)有8人回復(fù)

【求助】Fortran的計算結(jié)果已經(jīng)有11人回復(fù)
【求助】FORTRAN中增加計算區(qū)域后出錯，已經(jīng)有5人回復(fù)
【求助】求助個Fortran計算，應(yīng)該不難，半天搞不定，請高手幫忙已經(jīng)有11人回復(fù)

1樓 2011-11-08 15:26:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

huruiusst

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 2687.3
帖子: 67
在線: 115.8小時
蟲號: 1103368
注冊: 2010-09-19
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

★
jjdg(金幣+1): 感謝參與 2011-11-09 02:06:18

今天看到這帖子才知道還有這種情況學(xué)習(xí)了聲明和賦值分開寫比較保險

贊一下(2人)

回復(fù)此樓

SLS,MD，HRSG

11樓2011-11-08 17:06:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★
princeado(金幣+20): 謝了，問題解決了，我再好好研究下原因 2011-11-08 16:10:31
余澤成(金幣+1): 謝謝參與應(yīng)助！ 2011-11-08 18:10:24

內(nèi)容已刪除

贊一下(1人)

回復(fù)此樓

2樓2011-11-08 15:50:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hhucy

實(shí)習(xí)版主 (著名寫手)

DOE鍋爐工

應(yīng)助: 53 (初中生)
貴賓: 0.264
金幣: 7775.6
散金: 1861
紅花: 13
沙發(fā): 7
帖子: 2588
在線: 1275.5小時
蟲號: 1079650
注冊: 2010-08-22
性別: GG
管轄: 能源

引用回帖:

2樓: Originally posted by snoopyzhao at 2011-11-08 15:50:49:
把 jiecheng 中的下列語句
[cdoe]
integer :: ii=1
integer :: m=1
[/code]
修改成下面就可以了：

CODE:

      integer :: ii

      integer :: m

      m=1

至于為什么，可以參考 ...

好可怕的變量

回復(fù)此樓

人生那么多不確定，你怕什么

3樓2011-11-08 16:06:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

princeado

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 2729.5
紅花: 5
帖子: 874
在線: 868.6小時
蟲號: 1031978
注冊: 2010-05-30
性別: GG
專業(yè): 水文、水資源

引用回帖:

2樓: Originally posted by snoopyzhao at 2011-11-08 15:50:49:
把 jiecheng 中的下列語句
[cdoe]
integer :: ii=1
integer :: m=1
[/code]
修改成下面就可以了：

CODE:

      integer :: ii

      integer :: m

      m=1

至于為什么，可以參考 ...

原因不是很清楚，您能再幫我解釋一下嗎？

贊一下

回復(fù)此樓

4樓2011-11-08 16:18:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
余澤成(金幣+1): 謝謝參與應(yīng)助！ 2011-11-08 18:10:46

fortran 對聲明階段賦值的變量自動賦予 SAVE 屬性，這樣，在這個子程的連續(xù)調(diào)過程，在下一次執(zhí)行時，其初始值為一次調(diào)用結(jié)束時的值……

另外，子程序的 return 之后，不要用 stop，stop 只在應(yīng)該結(jié)束程序時才使用……

贊一下(1人)

回復(fù)此樓

5樓2011-11-08 16:18:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

princeado

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 2729.5
紅花: 5
帖子: 874
在線: 868.6小時
蟲號: 1031978
注冊: 2010-05-30
性別: GG
專業(yè): 水文、水資源

引用回帖:

5樓: Originally posted by snoopyzhao at 2011-11-08 16:18:53:
fortran 對聲明階段賦值的變量自動賦予 SAVE 屬性，這樣，在這個子程的連續(xù)調(diào)過程，在下一次執(zhí)行時，其初始值為一次調(diào)用結(jié)束時的值……

另外，子程序的 return 之后，不要用 stop，stop 只在應(yīng)該結(jié)束程序時才使 ...

您的意思是是不是我原來的程序每次調(diào)用完之后，m的值都是上一次結(jié)束時的m值，所以才會出錯？

贊一下

回復(fù)此樓

6樓2011-11-08 16:28:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

princeado

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 2729.5
紅花: 5
帖子: 874
在線: 868.6小時
蟲號: 1031978
注冊: 2010-05-30
性別: GG
專業(yè): 水文、水資源

引用回帖:

還有就是，integer :: m
m=1
與 integer :: m=1 兩種聲明變量的方式的區(qū)別是什么��？

贊一下

回復(fù)此樓

7樓2011-11-08 16:37:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

princeado

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 2729.5
紅花: 5
帖子: 874
在線: 868.6小時
蟲號: 1031978
注冊: 2010-05-30
性別: GG
專業(yè): 水文、水資源

引用回帖:

6樓: Originally posted by princeado at 2011-11-08 16:28:48:
您的意思是是不是我原來的程序每次調(diào)用完之后，m的值都是上一次結(jié)束時的m值，所以才會出錯？

我測試過了，我原來的程序是每次m的初始值是上次計算完的值。但是還是不明白兩種聲明方式的區(qū)別。

贊一下

回復(fù)此樓

8樓2011-11-08 16:39:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★ ★
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-11-08 18:10:59
余澤成: 專家考核存檔！ 2011-11-08 18:11:11

引用回帖:

7樓: Originally posted by princeado at 2011-11-08 16:37:35:
還有就是，integer :: m
m=1
與 integer :: m=1 兩種聲明變量的方式的區(qū)別是什么啊？

區(qū)別在于
integer :: m
m=1
中的 m是一個局變變量，不具有 SAVE 屬性，每次調(diào)用該子程時，m 都會被重新賦值。

而 integer :: m=1 中的 m 雖然是一個局部變量，但具有 SAVE 屬性，這里的 1 只是主程序初始化時的取值，在隨后的執(zhí)行中，它的值會是上一次該子程序執(zhí)行結(jié)束時的值。

至于什么是 SAVE 屬性，去查一下標(biāo)準(zhǔn)吧，我還真講不清楚。大概的意思是，該變量的最值取值會被保留下來，直到主程序結(jié)束。

贊一下(2人)

回復(fù)此樓

9樓2011-11-08 16:53:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

princeado

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 2729.5
紅花: 5
帖子: 874
在線: 868.6小時
蟲號: 1031978
注冊: 2010-05-30
性別: GG
專業(yè): 水文、水資源

引用回帖:

9樓: Originally posted by snoopyzhao at 2011-11-08 16:53:36:
區(qū)別在于
integer :: m
m=1
中的 m是一個局變變量，不具有 SAVE 屬性，每次調(diào)用該子程時，m 都會被重新賦值。

而 integer :: m=1 中的 m 雖然是一個局部變量，但具有 SAVE 屬性，這里的 1 只是主程序初始化 ...

好的，非常感謝了，我再琢磨琢磨啊！

贊一下

回復(fù)此樓

10樓2011-11-08 16:56:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 0703化學(xué)調(diào)劑 +6	妮妮ninicgb 2026-03-15	9/450	2026-03-16 16:40 by houyaoxu
[教師之家] 焦慮 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 求調(diào)劑（材料與化工327） +4	愛吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 314求調(diào)劑 +7	無懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 化工學(xué)碩306求調(diào)劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant