 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3376 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
為什么Fortran調(diào)試無錯誤運行不出結(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)，為啥這個程序運行后yingfeng.dat文件里面啥數(shù)據(jù)都沒有 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有249人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 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ù)
- 開始學(xué)習(xí)FORTRAN，結(jié)果殺毒軟件總是認(rèn)為生成的exe程序是木馬，要求刪掉，為什么呢？
已經(jīng)有7人回復(fù)
- fortran
已經(jīng)有4人回復(fù)
- Fortran77 代碼調(diào)試中提示有Syntax error，但卻不知道為什么出錯。望高人指點！
已經(jīng)有4人回復(fù)
- fortran循環(huán)求助
已經(jīng)有3人回復(fù)
- Fortran 求解復(fù)數(shù)矩陣SVD時，用到cgesvd，結(jié)果與matlab不同，為什么？
已經(jīng)有15人回復(fù)
- fortran里的intrisic function: precision()x結(jié)果表示什么意思
已經(jīng)有5人回復(fù)
- fortran 程序修改后改回來，計算結(jié)果不同？？？為什么會這樣
已經(jīng)有2人回復(fù)
- data 給fortran 的二維數(shù)組如何命名？為什么老出錯，解決送金幣
已經(jīng)有1人回復(fù)
- fortran輸出結(jié)果為0
已經(jīng)有3人回復(fù)
- fortran 程序有write和沒有write兩種情況計算的結(jié)果為什么會不一樣?
已經(jīng)有14人回復(fù)
- fortran輸出求助
已經(jīng)有6人回復(fù)
- 幫忙把公式轉(zhuǎn)化為FORTRAN子程序
已經(jīng)有6人回復(fù)
- fortran中計算數(shù)據(jù)寫入文本為 NAN的問題
已經(jīng)有4人回復(fù)
- 剛開始學(xué)Fortran 幫忙看一下這個程序 為什么算不出結(jié)果
已經(jīng)有3人回復(fù)
- fortran斷點
已經(jīng)有2人回復(fù)
- fortran 運行問題
已經(jīng)有2人回復(fù)
- 【求助】請問有什么軟件可以將C語言代碼轉(zhuǎn)換為fortran嗎？
已經(jīng)有7人回復(fù)
4樓2016-01-02 17:51:51
bennett_14
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 7204.8
- 紅花: 4
- 沙發(fā): 1
- 帖子: 932
- 在線: 421.5小時
- 蟲號: 1073877
- 注冊: 2010-08-13
- 性別: GG
- 專業(yè): 計算機科學(xué)
2樓2016-01-02 16:09:59
bennett_14
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 7204.8
- 紅花: 4
- 沙發(fā): 1
- 帖子: 932
- 在線: 421.5小時
- 蟲號: 1073877
- 注冊: 2010-08-13
- 性別: GG
- 專業(yè): 計算機科學(xué)
3樓2016-01-02 16:35:26
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
夜月丨噩夢: 金幣+30, ★★★很有幫助 2016-01-05 11:42:26
夜月丨噩夢: 金幣+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 ... 可以看到，這個程序能否向yingfeng.dat中寫入數(shù)據(jù)的關(guān)鍵是：t的數(shù)值相距0.5的差值是否小到0.0001以內(nèi)，和其它問題完全無關(guān)。 如果分析我給出的這些語句的話，就會發(fā)現(xiàn)問題所在： t從0.003開始累加，要和0.5比？？？能達到相差小于0.0001的程度么？ 相當(dāng)于500/3，得到的余數(shù)是多少？ 條件達不到，當(dāng)然不會輸出。程序在邏輯上是沒錯的。 希望樓主能再審視一下相應(yīng)的條件，讓程序輸出你要的結(jié)果。 另外，希望樓主能注意以下幾個方面： 1 t1和t0都是雙精度，但是t卻沒有定義�。。。“凑漳J(rèn)的I-N規(guī)則，t是單精度型，可能會引起誤差，從而引起不必要的麻煩哦。 2 open之后沒有對應(yīng)的close語句，這點很不好��！很不好�。『懿缓茫�！樓主目前寫的是小程序不覺得，程序一大，很容易搞混的，到時候想寫這個文件卻寫到那個文件去了，或者想寫卻發(fā)現(xiàn)文件沒打開之類的。而且這種錯誤會很難調(diào)試。 3 stop語句一般只用在這樣兩個方面：調(diào)試程序 或 出錯后的強行停止。不建議用stop來結(jié)束程序，讓程序自己運行完畢后停止豈不更好？ 4 前面有program xxx的話，最后建議用end program xxx結(jié)束。這樣構(gòu)造一個封閉的block，今后也會很容易讀代碼，便于維護管理。 5 樓主的程序使用了默認(rèn)的I-N規(guī)則，建議使用Implicit None，取消默認(rèn)規(guī)則，使每個變量都要定義類型（和C一樣強類型）。使用默認(rèn)規(guī)則的話，一旦有什么問題都是非常隱蔽極難調(diào)試的。 最后，樓主如果有程序方面的問題，歡迎移步“程序語言”版塊提問。 |
» 本帖已獲得的紅花（最新10朵）
夜月丨噩夢5樓2016-01-02 19:32:34
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +10 | 申子申申 2026-03-19 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿福大288有機化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|