當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

夜月丨噩夢(mèng)

應(yīng)助: 0 (幼兒園)
金幣: 694.8
散金: 50
紅花: 1
帖子: 20
在線: 17.4小時(shí)
蟲(chóng)號(hào): 4040736
注冊(cè): 2015-08-28
性別: GG
專業(yè): 半導(dǎo)體晶體與薄膜材料

[求助] 為什么Fortran調(diào)試無(wú)錯(cuò)誤運(yùn)行不出結(jié)果已有2人參與

program main
parameter  n=600
double precision u(0:n),u1(0:n),x0,x1,x,t0,t1,h,dett,lambda,y(0:n)
integer i,j
x0=0.0d0
x1=6.0d0
t0=0.0d0
lambda=0.3d0
t1=0.5d0
h=0.01d0
dett=0.003d0
x=x0
!do i=0,n,1
   !if(i<=250)then
   !y(i)=1.0
   !else
! y(i)=0
   !end if
! end do
do i=0,n,1
   if(x>=1 .and. x<=3) then
   u(i)=1.0d0
   else
   u(i)=0.0d0
   end if
   x=x+h
   end do
   open(1,file="yingfeng.dat"

   t=t0+dett
   x=x0
   do while(t<=t1)
      x=x0+h
   do j=1,n,1
   u1(j)=u(j)-3*lambda*(u(j)-u(j-1))

   if(abs(t-0.5)<=0.0001)then
   write(1,*)x,u1(j),y(j)
   end if
   x=x+h
      end do
      do j=1,n-1
   u(j)=u1(j)
   u(0)=1
   u(1)=0
   end do

      t=t+dett
   end do
stop
end
求大神指導(dǎo)，為啥這個(gè)程序運(yùn)行后yingfeng.dat文件里面啥數(shù)據(jù)都沒(méi)有

回復(fù)此樓

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有299人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有9人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

Intel fortran與visual fortran 已經(jīng)有11人回復(fù)
FORTRAN 已經(jīng)有0人回復(fù)
Fortran 已經(jīng)有7人回復(fù)
已經(jīng)安裝studio，為什么還是不能安裝Fortran呢？已經(jīng)有0人回復(fù)
開(kāi)始學(xué)習(xí)FORTRAN，結(jié)果殺毒軟件總是認(rèn)為生成的exe程序是木馬，要求刪掉，為什么呢？已經(jīng)有7人回復(fù)
fortran 已經(jīng)有4人回復(fù)
Fortran77 代碼調(diào)試中提示有Syntax error，但卻不知道為什么出錯(cuò)。望高人指點(diǎn)！已經(jīng)有4人回復(fù)
fortran循環(huán)求助已經(jīng)有3人回復(fù)
Fortran 求解復(fù)數(shù)矩陣SVD時(shí)，用到cgesvd，結(jié)果與matlab不同，為什么？已經(jīng)有15人回復(fù)
fortran里的intrisic function: precision()x結(jié)果表示什么意思已經(jīng)有5人回復(fù)
fortran 程序修改后改回來(lái)，計(jì)算結(jié)果不同？？？為什么會(huì)這樣已經(jīng)有2人回復(fù)
data 給fortran 的二維數(shù)組如何命名？為什么老出錯(cuò)，解決送金幣已經(jīng)有1人回復(fù)
fortran輸出結(jié)果為0 已經(jīng)有3人回復(fù)
fortran 程序有write和沒(méi)有write兩種情況計(jì)算的結(jié)果為什么會(huì)不一樣? 已經(jīng)有14人回復(fù)
fortran輸出求助已經(jīng)有6人回復(fù)
幫忙把公式轉(zhuǎn)化為FORTRAN子程序已經(jīng)有6人回復(fù)
fortran中計(jì)算數(shù)據(jù)寫入文本為 NAN的問(wèn)題已經(jīng)有4人回復(fù)
剛開(kāi)始學(xué)Fortran 幫忙看一下這個(gè)程序為什么算不出結(jié)果已經(jīng)有3人回復(fù)
fortran斷點(diǎn) 已經(jīng)有2人回復(fù)
fortran 運(yùn)行問(wèn)題已經(jīng)有2人回復(fù)
【求助】請(qǐng)問(wèn)有什么軟件可以將C語(yǔ)言代碼轉(zhuǎn)換為fortran嗎？已經(jīng)有7人回復(fù)

1樓 2015-12-29 21:04:08

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

fmying

銀蟲(chóng) (小有名氣)

應(yīng)助: 29 (小學(xué)生)
金幣: 618
散金: 50
紅花: 6
帖子: 127
在線: 113.6小時(shí)
蟲(chóng)號(hào): 1081069
注冊(cè): 2010-08-24
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
夜月丨噩夢(mèng): 金幣+30, ★★★很有幫助 2016-01-05 11:42:26

你的程序中最重要的部分如下：
double precision ...t0,t1,...,dett...
t0=0.0d0
t1=0.5d0
dett=0.003d0
...
t=t0+dett
do while(t<=t1)
...
if(abs(t-0.5)<=0.0001) then
write(1,*)x,u1(j),y(j)
endif
...

可以看到，這個(gè)程序能否向yingfeng.dat中寫入數(shù)據(jù)的關(guān)鍵是：t的數(shù)值相距0.5的差值是否小到0.0001以內(nèi)，和其它問(wèn)題完全無(wú)關(guān)。

如果分析我給出的這些語(yǔ)句的話，就會(huì)發(fā)現(xiàn)問(wèn)題所在：
t從0.003開(kāi)始累加，要和0.5比？？？能達(dá)到相差小于0.0001的程度么？
相當(dāng)于500/3，得到的余數(shù)是多少？
條件達(dá)不到，當(dāng)然不會(huì)輸出。程序在邏輯上是沒(méi)錯(cuò)的。
希望樓主能再審視一下相應(yīng)的條件，讓程序輸出你要的結(jié)果。

另外，希望樓主能注意以下幾個(gè)方面：
1 t1和t0都是雙精度，但是t卻沒(méi)有定義�。。�！按照默認(rèn)的I-N規(guī)則，t是單精度型，可能會(huì)引起誤差，從而引起不必要的麻煩哦。
2 open之后沒(méi)有對(duì)應(yīng)的close語(yǔ)句，這點(diǎn)很不好��！很不好��！很不好��！樓主目前寫的是小程序不覺(jué)得，程序一大，很容易搞混的，到時(shí)候想寫這個(gè)文件卻寫到那個(gè)文件去了，或者想寫卻發(fā)現(xiàn)文件沒(méi)打開(kāi)之類的。而且這種錯(cuò)誤會(huì)很難調(diào)試。
3 stop語(yǔ)句一般只用在這樣兩個(gè)方面：調(diào)試程序或出錯(cuò)后的強(qiáng)行停止。不建議用stop來(lái)結(jié)束程序，讓程序自己運(yùn)行完畢后停止豈不更好？
4 前面有program xxx的話，最后建議用end program xxx結(jié)束。這樣構(gòu)造一個(gè)封閉的block，今后也會(huì)很容易讀代碼，便于維護(hù)管理。
5 樓主的程序使用了默認(rèn)的I-N規(guī)則，建議使用Implicit None，取消默認(rèn)規(guī)則，使每個(gè)變量都要定義類型（和C一樣強(qiáng)類型）。使用默認(rèn)規(guī)則的話，一旦有什么問(wèn)題都是非常隱蔽極難調(diào)試的。

最后，樓主如果有程序方面的問(wèn)題，歡迎移步“程序語(yǔ)言”版塊提問(wèn)。

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

夜月丨噩夢(mèng)

上善若水

5樓2016-01-02 19:32:34

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

查看全部 10 個(gè)回答

bennett_14

木蟲(chóng) (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 7204.8
紅花: 4
沙發(fā): 1
帖子: 932
在線: 421.5小時(shí)
蟲(chóng)號(hào): 1073877
注冊(cè): 2010-08-13
性別: GG
專業(yè): 計(jì)算機(jī)科學(xué)

x0=0.0d0中的d0是啥？？？？

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

君子當(dāng)自強(qiáng)不息，死而后已！

2樓2016-01-02 16:09:59

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

bennett_14

木蟲(chóng) (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 7204.8
紅花: 4
沙發(fā): 1
帖子: 932
在線: 421.5小時(shí)
蟲(chóng)號(hào): 1073877
注冊(cè): 2010-08-13
性別: GG
專業(yè): 計(jì)算機(jī)科學(xué)

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

經(jīng)過(guò)測(cè)試，可以生成yingfeng.dat文件的，文件和程序在相同的目錄下，估計(jì)你沒(méi)有運(yùn)行，或運(yùn)行了但是找錯(cuò)了文件夾

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

君子當(dāng)自強(qiáng)不息，死而后已！

3樓2016-01-02 16:35:26

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

夜月丨噩夢(mèng)

金蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 694.8
散金: 50
紅花: 1
帖子: 20
在線: 17.4小時(shí)
蟲(chóng)號(hào): 4040736
注冊(cè): 2015-08-28
性別: GG
專業(yè): 半導(dǎo)體晶體與薄膜材料

引用回帖:

3樓: Originally posted by bennett_14 at 2016-01-02 16:35:26
經(jīng)過(guò)測(cè)試，可以生成yingfeng.dat文件的，文件和程序在相同的目錄下，估計(jì)你沒(méi)有運(yùn)行，或運(yùn)行了但是找錯(cuò)了文件夾

是可以運(yùn)行，但是運(yùn)行結(jié)果是0KB，里面啥都沒(méi)有

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

4樓2016-01-02 17:51:51

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

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

[考研] 286分調(diào)劑 +10	Faune 2026-03-30	11/550	2026-03-30 20:35 by 啊李999
[考研] 調(diào)劑 +4	GK72 2026-03-30	4/200	2026-03-30 20:32 by dick_runner
[考研] 0703化學(xué)321分求調(diào)劑 +10	三dd. 2026-03-30	11/550	2026-03-30 19:24 by markhwc
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[考研] 化工專碩求調(diào)劑 +5	question挽風(fēng) 2026-03-24	5/250	2026-03-30 18:34 by tianyyysss
[考研] 一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	7/350	2026-03-30 18:29 by 兔兔兔111223314
[考研] 337求調(diào)劑 +6	《樹(shù)》 2026-03-29	6/300	2026-03-30 10:15 by herarysara
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8	柳溪邊 2026-03-26	8/400	2026-03-29 20:42 by 唐沐兒
[考研] 299求調(diào)劑 +10	15188958825 2026-03-25	10/500	2026-03-29 17:51 by 王亮_大連醫(yī)科大
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12	c297914 2026-03-28	12/600	2026-03-29 16:21 by Serene1974
[考研] 0856材料化工調(diào)劑總分330 +14	zhubinhao 2026-03-27	14/700	2026-03-29 10:01 by Sjndkwm
[考研] 352分-085602-一志愿985 +5	海納百川Ly 2026-03-29	5/250	2026-03-29 09:57 by Sjndkwm
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 085602 307分求調(diào)劑 +7	不知道叫什么！ 2026-03-26	7/350	2026-03-28 09:57 by 神馬都不懂
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3	熠然好運(yùn)氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito
[考研] 274求調(diào)劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3	zzffylgg 2026-03-24	3/150	2026-03-27 10:42 by shangxh
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325