 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3372 | 回復: 9 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
為什么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 求大神指導，為啥這個程序運行后yingfeng.dat文件里面啥數(shù)據(jù)都沒有 |
回復此樓» 猜你喜歡
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有7人回復
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有53人回復
-
求調(diào)劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復
-
光學工程學碩調(diào)劑信息
已經(jīng)有31人回復
-
歡迎加入課題組
已經(jīng)有0人回復
-
散金幣，求好運，祝面上順利！
已經(jīng)有68人回復
-
華東師范大學芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復
-
調(diào)劑
已經(jīng)有0人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- Intel fortran與visual fortran
已經(jīng)有11人回復
- FORTRAN
已經(jīng)有0人回復
- Fortran
已經(jīng)有7人回復
- 已經(jīng)安裝studio，為什么還是不能安裝Fortran呢？
已經(jīng)有0人回復
- 開始學習FORTRAN，結(jié)果殺毒軟件總是認為生成的exe程序是木馬，要求刪掉，為什么呢？
已經(jīng)有7人回復
- fortran
已經(jīng)有4人回復
- Fortran77 代碼調(diào)試中提示有Syntax error，但卻不知道為什么出錯。望高人指點！
已經(jīng)有4人回復
- fortran循環(huán)求助
已經(jīng)有3人回復
- Fortran 求解復數(shù)矩陣SVD時，用到cgesvd，結(jié)果與matlab不同，為什么？
已經(jīng)有15人回復
- fortran里的intrisic function: precision()x結(jié)果表示什么意思
已經(jīng)有5人回復
- fortran 程序修改后改回來，計算結(jié)果不同？？？為什么會這樣
已經(jīng)有2人回復
- data 給fortran 的二維數(shù)組如何命名？為什么老出錯，解決送金幣
已經(jīng)有1人回復
- fortran輸出結(jié)果為0
已經(jīng)有3人回復
- fortran 程序有write和沒有write兩種情況計算的結(jié)果為什么會不一樣?
已經(jīng)有14人回復
- fortran輸出求助
已經(jīng)有6人回復
- 幫忙把公式轉(zhuǎn)化為FORTRAN子程序
已經(jīng)有6人回復
- fortran中計算數(shù)據(jù)寫入文本為 NAN的問題
已經(jīng)有4人回復
- 剛開始學Fortran 幫忙看一下這個程序 為什么算不出結(jié)果
已經(jīng)有3人回復
- fortran斷點
已經(jīng)有2人回復
- fortran 運行問題
已經(jīng)有2人回復
- 【求助】請問有什么軟件可以將C語言代碼轉(zhuǎn)換為fortran嗎？
已經(jīng)有7人回復
bennett_14
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學生)
- 金幣: 7204.8
- 紅花: 4
- 沙發(fā): 1
- 帖子: 932
- 在線: 421.5小時
- 蟲號: 1073877
- 注冊: 2010-08-13
- 性別: GG
- 專業(yè): 計算機科學
3樓2016-01-02 16:35:26
bennett_14
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學生)
- 金幣: 7204.8
- 紅花: 4
- 沙發(fā): 1
- 帖子: 932
- 在線: 421.5小時
- 蟲號: 1073877
- 注冊: 2010-08-13
- 性別: GG
- 專業(yè): 計算機科學
2樓2016-01-02 16:09:59
4樓2016-01-02 17:51:51
【答案】應(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的程度么？ 相當于500/3，得到的余數(shù)是多少？ 條件達不到，當然不會輸出。程序在邏輯上是沒錯的。 希望樓主能再審視一下相應(yīng)的條件，讓程序輸出你要的結(jié)果。 另外，希望樓主能注意以下幾個方面： 1 t1和t0都是雙精度，但是t卻沒有定義�。。�！按照默認的I-N規(guī)則，t是單精度型，可能會引起誤差，從而引起不必要的麻煩哦。 2 open之后沒有對應(yīng)的close語句，這點很不好��！很不好�。『懿缓茫。侵髂壳皩懙氖切〕绦虿挥X得，程序一大，很容易搞混的，到時候想寫這個文件卻寫到那個文件去了，或者想寫卻發(fā)現(xiàn)文件沒打開之類的。而且這種錯誤會很難調(diào)試。 3 stop語句一般只用在這樣兩個方面：調(diào)試程序 或 出錯后的強行停止。不建議用stop來結(jié)束程序，讓程序自己運行完畢后停止豈不更好？ 4 前面有program xxx的話，最后建議用end program xxx結(jié)束。這樣構(gòu)造一個封閉的block，今后也會很容易讀代碼，便于維護管理。 5 樓主的程序使用了默認的I-N規(guī)則，建議使用Implicit None，取消默認規(guī)則，使每個變量都要定義類型（和C一樣強類型）。使用默認規(guī)則的話，一旦有什么問題都是非常隱蔽極難調(diào)試的。 最后，樓主如果有程序方面的問題，歡迎移步“程序語言”版塊提問。 |
» 本帖已獲得的紅花（最新10朵）
夜月丨噩夢5樓2016-01-02 19:32:34
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學071000 329分求調(diào)劑 +4 | 我愛生物生物愛?/a> 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學獎覬?/a> 2026-03-15 | 8/400 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|