计算结果出现问题

dididada

如下图的源代码所示，我要用fortran计算一个有初值和边界条件限制的微分方程，但结果一直不如心意，希望能够在论坛得到解答。


主程序                 子程序

pasuka

数理方程的课程大作业吗？
建议不要纠结编程，先把梁昆淼老师那本教材的课后习题都做一遍或者答案认真抄三遍
很多时候，习题不做的滚瓜烂熟或者解题过程倒背如流，编程计算就是水中捞月

dididada

pasuka 发表于 2018-9-25 16:20
数理方程的课程大作业吗？
建议不要纠结编程，先把梁昆淼老师那本教材的课后习题都做一遍或者答案认真抄三 ...

不是作业，是导师让我入门做的东西。原理都很清楚，调试也是正确的，就是结果不是很满意

pasuka

为啥不去问导师呢？

vvt

没看见代码呀，楼主自己能看见？

dididada

不好意思，没把代码贴全。

[Fortran] 纯文本查看 复制代码

program main
  implicit none
  integer,parameter::N=11
  integer t
 real*8::A(N),B(N),C(N),D(N),H(N)
 integer,parameter::Dx=10,Dt=1
  !空间和时间的步长
 real,parameter:: k=0.1
 real,parameter:: lemad=real(k*Dt)/real(Dx**2)
  !含水层参数
 integer i,j
  !声明
    do i=2,N-1
     A(i)=-lemad
     B(i)=1+2*lemad
     C(i)=-lemad
     end do
    write(*,20)lemad
     20	format(/,5X,'lemad=',F9.6)

    do i =1,N
    H(i)=10
    end do
    write (*,*) H(2:N-1)
    H(1)=10
    H(N)=12
    write(*,15)
     15   format(5X,'TIME',5X,'HEAD')!赋初值和边界值
     t = 0
    do j=1,101
    t=t+Dt 
    do  i=2,N-1
       D(i)=H(i)
    D(2)=D(2)+H(1)*lemad
    D(N-1)=D(N-1)+H(N)*lemad
    end do
    call  chase (N,A,B,C,D,H)
     write(*,30) t,H(i)
    30	format(1X,I3,5X,11F9.6)
    end do
    read(*,*)t
    end program

    subroutine chase (N,A,B,C,D,H)
     !追赶法计算三对角线方程组
     implicit none
      integer:: N
     real*8::A(N),B(N),C(N),D(N),H(N)   
     real*8::L(N),U(N)
    integer::i,j
     L(2)=B(2)
     U(2)=D(2)/L(2)
     do  i = 3,N-1
      L(i) = B(i)-A(i)*C(i-1)/L(i-1)
      U(i) =(D(i)-A(i)*U(i-1))/L(i)
   end do
    !回代求H
    H(N-1)=U(N-1)
     do j= 2,N-3
     i=N-1-j
     H(i)=U(i)-C(i)*H(i+1)/L(i)
    end do
    return
    end subroutine chase

dididada

vvt 发表于 2018-9-26 06:44
没看见代码呀，楼主自己能看见？

不好意思，我现在把代码贴出来了，您能帮我看一下吗

dididada

pasuka 发表于 2018-9-25 21:53
为啥不去问导师呢？

老师最近有点忙，人还不在国内

pasuka

dididada 发表于 2018-9-26 08:54
老师最近有点忙，人还不在国内

老板不在，总有小老板在吧？
没有小老板，总有师兄师姐吧？
另外，lz就不能写个汇报总结把已经解决的和亟待解决的问题列举出来，再发送邮件给老板？

dididada

pasuka 发表于 2018-9-26 11:03
老板不在，总有小老板在吧？
没有小老板，总有师兄师姐吧？
另外，lz就不能写个汇报总结把已经解决的和亟 ...

哈哈，好的，你说的很有道理。