版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

下一頁(yè)

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線(xiàn): 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專(zhuān)業(yè): 物理電子學(xué)

[交流] 【求助】fortran出錯(cuò) 編譯通過(guò) 運(yùn)行溢出

subroutine  adams(t,h,m,y,n,z,f,d,eps,b,e,s,g,kp1,kp2)
external f,rkt
dimension y(n),d(n),z(n,m),b(4,n),e(n),s(n),g(n)
   real*8 kp1,kp2
a=t
do 10 i=1,n
10    z(i,1)=y(i)
   call f(t,y,n,d,kp1,kp2)
      do 20 i=1,n
20    b(1,i)=d(i)
      do 50 i=2,4
         if (i.le.m)then
      t=a+(i-1)*h
      call rkt(t,h,y,n,f,eps,d,e,s,g,kp1,kp2)
         do 30 k=1,n
30          z(k,i)=y(k)
         call f(t,y,n,d,kp1,kp2)
         do 40 k=1,n
40          b(i,k)=d(k)
         end if
50    continue
   do 100 i=5,m
      do 60 j=1,n
   q=55*b(4,j)-59*b(3,j)+37*b(2,j)-9*b(1,j)
   y(j)=z(j,i-1)+h*q/24.0
   b(1,j)=b(2,j)
   b(2,j)=b(3,j)
   b(3,j)=b(4,j)
60       continue
         t=a+(i-1)*h
      call f(t,y,n,d,kp1,kp2)
      do 70 k=1,n
70          b(4,k)=d(k)
               do 80 j=1,n
               q=9*b(4,j)+19*b(3,j)-5*b(2,j)+b(1,j)
               z(j,i)=z(j,i-1)+h*q/24.0
               y(j)=z(j,i)
80                continue
            do 90 k=1,n
90          b(4,k)=d(k)
100    continue
   t=a
   return
   end


   subroutine rkt(t,h,y,n,f,eps,d,b,c,g,kp1,kp2)
   dimension y(n),d(n),a(4),b(n),c(n),g(n)
         real*8 kp1,kp2
   h1=h
   m=1
   p=1+eps
   x=t
   do 10 i=1,n
10    c(i)=y(i)
20    if (p.ge.eps)then
         a(1)=h1/2.0
         a(2)=a(1)
         a(3)=h1
         a(4)=h1
         do 30 i=1,n
         g(i)=y(i)
            y(i)=c(i)
30          continue
            dt=h/m
            t=x
            do 80 j=1,m
            call f(t,y,n,d,kp1,kp2)
            do 40 i=1,n
40             b(i)=y(i)
               do 60 k=1,3
               do 50 i=1,n
                  y(i)=y(i)+a(k)*d(i)
                     b(i)=b(i)+a(k+1)*d(i)/3.0
50                continue
                  tt=t+a(k)
                  call f(tt,y,n,d,kp1,kp2)
60             continue
               do 70 i=1,n
70                   y(i)=b(i)+h1*d(i)/6.0
                  t=t+dt
80          continue
      p=0.0
      do 90 i=1,n
q=abs(y(i)-g(i))
if(q.gt.p) p=q
90    continue
      h1=h1/2.0
      m=m+m
      goto 20
   endif
t=x
return
end

program madams
c    -----------------------------------------------------------------------
c    this rpogram is used to solve a systwem of differential equagtion
c    by using adams method
c    -------------------------------------------------------------------
   external f
dimension  y(11),d(11),z(11,500000),b(4,11),e(11),s(11),g(11)
real*8 wucha,kp1,kp2,t,tao,variable,Parameters

   open(26,file='variable.dat')
open(28,file='Parameters.dat')
   open(10,file='wucha.dat')
h=0.001d0

y(1)=0.01d0
y(2)=0.02d0
y(3)=0.03d0

y(4)=0.4d0
y(5)=0.2d0
y(6)=0.8d0

y(7)=0.8d0
y(8)=2.5d0
y(9)=0.9d0
y(10)=6d0
y(11)=1d0

m=500000
n=11
eps=0.00001d0

   do kp1=1.50d0,1.60d0,0.1d0
   do kp2=0.00d0,0.50d0,0.1d0
call adams(a,h,m,y,n,z,f,d,eps,b,e,s,g,kp1,kp2)
do i=1,m
   t=(i-1)*h
   detae=((z(1,i)-z(4,i))**2.0+(z(2,i)+z(5,i))**2.0+
   $ (z(3,i)-z(6,i))**2.0)**0.5

!    write(26,'(2f15.6)')t,detae
!    write(28,'(6f15.6)')t,z(7,i),z(8,i),z(9,i),z(10,i),z(11,i)
if(detae>50)then
write(10,*)kp1,kp2,1
   write(*,*)kp1,kp2,1
goto 20
elseif(detae>1.0d-3)then
tao=t
endif
enddo
if(tao<=50)then
write(10,*)kp1,kp2,0
   write(*,*)kp1,kp2,0
else
write(10,*)kp1,kp2,1
   write(*,*)kp1,kp2,1
endif
20    continue
   enddo
enddo
   close(2)
   end

subroutine f(t,y,n,d,kp1,kp2)
c    ----------------------------------------------------------------------
! real*8 kp1,kp2
   dimension  y(n),d(n)
double precision ex,ey,ez,u1,kp1,kp2
cccccccccccccc
   ex=y(1)-y(4)
   ey=y(2)+y(5)
   ez=y(3)-y(6)

c       u1=ey-ez+y(8)*(y(1)+y(4))*ex-y(10)*(y(1)+y(4))*ey+4d0*y(11)*ez
c       u2=0d0

cc    u1=ey+y(8)*(y(1)+y(4))*ex-y(10)*(y(1)+y(4))*ey
cc    u2=-ex+4d0*y(11)*ex

      u1=ey+y(8)*(y(1)+y(4))*ex-y(10)*(y(1)+y(4))*ey+4d0*y(11)*ez
      u2=-ex
!    do kp1=0.00d0,10.0d0,0.01d0
!    do kp2=0.00d0,10.0d0,0.01d0
! kp1=2.0d0
!  kp2=0.5d0
cccccccccccccc
   d(1)=y(2)-1d0*y(1)**3+3d0*y(1)**2-y(3)+3d0

   d(2)=1d0-5d0*y(1)**2-y(2)

   d(3)=0.006d0*(4d0*(y(1)+1.6d0)-y(3))

cccccccccccccccccccccccc
   d(4)=y(5)-y(7)*y(4)**3+y(8)*y(4)**2-y(6)+3d0+u1

   d(5)=y(9)-y(10)*y(4)**2-y(5)

   d(6)=y(11)*(4d0*(y(4)+1.6d0)-y(6))+u2
cccccccccccccccccccccccccccccccccc

d(7)=-kp1*kp2*y(1)**3*ex

d(8)=kp1*kp2*y(1)**2*ex

d(9)=kp1*kp2*ey

d(10)=-kp1*kp2*y(1)**2*ey

   d(11)=kp1*kp2*(4d0*(y(1)+1.6d0)-y(3))*ez
  !    enddo
! enddo

return
end
高手看看這個(gè)程序出錯(cuò)了為什么  溢出了�。�

回復(fù)此樓

» 猜你喜歡

0856材料專(zhuān)碩353求調(diào)劑已經(jīng)有4人回復(fù)
山東省面上項(xiàng)目限額評(píng)審已經(jīng)有4人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有5人回復(fù)
生物學(xué)071000 329分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿華中科技大學(xué)071000，求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)

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

【求助】求fortra程序已經(jīng)有4人回復(fù)

1樓 2010-09-14 19:58:52

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

snoopyzhao

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

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

★
余澤成(金幣+1):謝謝參與應(yīng)助！ 2010-09-15 13:32:14

我這里很順利地結(jié)束了啊，呵呵……

我用 gfortran 編譯的……

贊一下(1人)

回復(fù)此樓

2樓2010-09-14 20:22:34

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

maomao1210

金蟲(chóng) (正式寫(xiě)手)

程序強(qiáng)帖: 5
應(yīng)助: 2 (幼兒園)
金幣: 1432.8
散金: 242
紅花: 16
沙發(fā): 1
帖子: 991
在線(xiàn): 442.1小時(shí)
蟲(chóng)號(hào): 253215
注冊(cè): 2006-05-20
性別: MM
專(zhuān)業(yè): 考古理論

★
余澤成(金幣+1):謝謝參與應(yīng)助！ 2010-09-15 13:32:24

do kp1=1.50d0,1.60d0,0.1d0
      do kp2=0.00d0,0.50d0,0.1d0
      call adams(a,h,m,y,n,z,f,d,eps,b,e,s,g,kp1,kp2)
      do i=1,m
         t=(i-1)*h
   detae=((z(1,i)-z(4,i))**2.0+(z(2,i)+z(5,i))**2.0+
   $       (z(3,i)-z(6,i))**2.0)**0.5

!    write(26,'(2f15.6)')t,detae
!    write(28,'(6f15.6)')t,z(7,i),z(8,i),z(9,i),z(10,i),z(11,i)
      if(detae>50)then
      write(10,*)kp1,kp2,1
   write(*,*)kp1,kp2,1
      goto 20
      elseif(detae>1.0d-3)then
      tao=t
      endif
      enddo
      if(tao<=50)then
      write(10,*)kp1,kp2,0
   write(*,*)kp1,kp2,0
      else
      write(10,*)kp1,kp2,1
   write(*,*)kp1,kp2,1
      endif
20    continue
   enddo
      enddo

在循環(huán)的時(shí)候，=((z(1,i)-z(4,i))**2.0+(z(2,i)+z(5,i))**2.0+
   $       (z(3,i)-z(6,i))**2.0
是一個(gè)無(wú)窮大或者無(wú)窮小的數(shù)字，cvf不識(shí)別，不予計(jì)算，就提示溢出了。

gfortran 編譯計(jì)算的時(shí)候，即使是NaN 也會(huì)繼續(xù)計(jì)算不提示任何錯(cuò)誤。

贊一下(1人)

回復(fù)此樓

3樓2010-09-14 22:22:30

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

snoopyzhao

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

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

引用回帖:

在循環(huán)的時(shí)候，=((z(1,i)-z(4,i))**2.0+(z(2,i)+z(5,i))**2.0+
$ (z(3,i)-z(6,i))**2.0
是一個(gè)無(wú)窮大或者無(wú)窮小的數(shù)字，cvf不識(shí)別，不予計(jì)算，就提示溢出了。

gfortran 編譯計(jì)算的時(shí)候，即使是NaN 也會(huì)繼續(xù)計(jì)算不提示任何錯(cuò)誤。

多謝指教……

回復(fù)此樓

4樓2010-09-14 22:43:47

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

maomao1210

金蟲(chóng) (正式寫(xiě)手)

程序強(qiáng)帖: 5
應(yīng)助: 2 (幼兒園)
金幣: 1432.8
散金: 242
紅花: 16
沙發(fā): 1
帖子: 991
在線(xiàn): 442.1小時(shí)
蟲(chóng)號(hào): 253215
注冊(cè): 2006-05-20
性別: MM
專(zhuān)業(yè): 考古理論

引用回帖:

Originally posted by snoopyzhao at 2010-09-14 22:43:47:

多謝指教……

snoopyzhao ，嘿嘿，

我猜的，我沒(méi)試

，哈哈，在宿舍不方便調(diào)試。哈哈。

贊一下

回復(fù)此樓

5樓2010-09-14 22:57:57

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

zyj8119

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線(xiàn): 1329.4小時(shí)
蟲(chóng)號(hào): 664177
注冊(cè): 2008-11-29
性別: GG
專(zhuān)業(yè): 理論和計(jì)算化學(xué)

SQRT范圍錯(cuò)誤？就是SQRT的數(shù)值里面出現(xiàn)了負(fù)數(shù)？

贊一下

回復(fù)此樓

好好學(xué)習(xí)，天天向上。

6樓2010-09-15 08:54:13

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

maomao1210

金蟲(chóng) (正式寫(xiě)手)

程序強(qiáng)帖: 5
應(yīng)助: 2 (幼兒園)
金幣: 1432.8
散金: 242
紅花: 16
沙發(fā): 1
帖子: 991
在線(xiàn): 442.1小時(shí)
蟲(chóng)號(hào): 253215
注冊(cè): 2006-05-20
性別: MM
專(zhuān)業(yè): 考古理論

建議樓上的看看程序再發(fā)言

，
實(shí)數(shù)的平方和怎么可能是負(fù)數(shù)。

贊一下

回復(fù)此樓

7樓2010-09-15 09:07:29

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線(xiàn): 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專(zhuān)業(yè): 物理電子學(xué)

引用回帖:

Originally posted by snoopyzhao at 2010-09-14 20:22:34:
我這里很順利地結(jié)束了啊，呵呵……

我用 gfortran 編譯的……

運(yùn)行結(jié)果出現(xiàn)NAN啊你那沒(méi)有？

贊一下

回復(fù)此樓

8樓2010-09-15 09:22:15

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線(xiàn): 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專(zhuān)業(yè): 物理電子學(xué)

引用回帖:

Originally posted by maomao1210 at 2010-09-14 22:22:30:
do kp1=1.50d0,1.60d0,0.1d0
      do kp2=0.00d0,0.50d0,0.1d0
      call adams(a,h,m,y,n,z,f,d,eps,b,e,s,g,kp1,kp2)
      do i=1,m
         t=(i-1)*h
   detae=((z(1,i)-z(4,i))**2.0+( ...

如何解決這個(gè)問(wèn)題呢

贊一下

回復(fù)此樓

9樓2010-09-15 09:23:21

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

maomao1210

金蟲(chóng) (正式寫(xiě)手)

程序強(qiáng)帖: 5
應(yīng)助: 2 (幼兒園)
金幣: 1432.8
散金: 242
紅花: 16
沙發(fā): 1
帖子: 991
在線(xiàn): 442.1小時(shí)
蟲(chóng)號(hào): 253215
注冊(cè): 2006-05-20
性別: MM
專(zhuān)業(yè): 考古理論

★
余澤成(金幣+1):辛苦！ 2010-09-15 13:33:10

引用回帖:

Originally posted by brqhl_ing at 2010-09-15 09:23:21:

如何解決這個(gè)問(wèn)題呢

微分方程組的那一塊有點(diǎn)問(wèn)題，
或者給定解微分方程組的參數(shù)有問(wèn)題。
或者這塊缺少一些句子之類(lèi)。
具體需要你自己仔細(xì)排查。。。。

贊一下(1人)

回復(fù)此樓

10樓2010-09-15 09:48:19

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 brqhl_ing 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4	zxts12 2026-03-21	4/200	2026-03-22 07:20 by ilovexiaobin
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 326求調(diào)劑 +5	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	8/400	2026-03-21 19:33 by ColorlessPI
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 材料工程（專(zhuān)）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 北科281學(xué)碩材料求調(diào)劑 +5	tcxiaoxx 2026-03-20	5/250	2026-03-20 21:35 by laoshidan
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無(wú)懈可擊111
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽(yáng)w
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

24小時(shí)熱門(mén)版塊排行榜

[交流] 【求助】fortran出錯(cuò) 編譯通過(guò) 運(yùn)行溢出

» 猜你喜歡

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

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