计算结果出现问题
如下图的源代码所示，我要用fortran计算一个有初值和边界条件限制的微分方程，但结果一直不如心意，希望能够在论坛得到解答。主程序 子程序
数理方程的课程大作业吗？
建议不要纠结编程，先把梁昆淼老师那本教材的课后习题都做一遍或者答案认真抄三遍
很多时候，习题不做的滚瓜烂熟或者解题过程倒背如流，编程计算就是水中捞月 pasuka 发表于 2018-9-25 16:20
数理方程的课程大作业吗？
建议不要纠结编程，先把梁昆淼老师那本教材的课后习题都做一遍或者答案认真抄三 ...
不是作业，是导师让我入门做的东西。原理都很清楚，调试也是正确的，就是结果不是很满意 为啥不去问导师呢？ 没看见代码呀，楼主自己能看见？ 不好意思，没把代码贴全。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
doi=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
callchase (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)
doi = 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
vvt 发表于 2018-9-26 06:44
没看见代码呀，楼主自己能看见？
不好意思，我现在把代码贴出来了，您能帮我看一下吗 pasuka 发表于 2018-9-25 21:53
为啥不去问导师呢？
老师最近有点忙，人还不在国内 dididada 发表于 2018-9-26 08:54
老师最近有点忙，人还不在国内
老板不在，总有小老板在吧？
没有小老板，总有师兄师姐吧？
另外，lz就不能写个汇报总结把已经解决的和亟待解决的问题列举出来，再发送邮件给老板？
pasuka 发表于 2018-9-26 11:03
老板不在，总有小老板在吧？
没有小老板，总有师兄师姐吧？
另外，lz就不能写个汇报总结把已经解决的和亟 ...
哈哈，好的，你说的很有道理。
页:
[1]