24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 989  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

physliebe

新虫 (小有名气)

[求助] 找不出毛病不知怎么改? 已有1人参与

计算系数矩阵为上三角的线性方程组Ax=b的解,N维
代码:
!##############################################
module tri_eq
contains
subroutine uptri(A,b,x,N)
implicit real*8(a-z) !a到z开头的变量都视为双精度浮点数
integer::i,j,k,N
real*8::A(N,N),b(N),x(N)
x(N)=b(N)/A(N,N)
!回代部分
do i=N-1,1,-1
do j=i+1,N
x(i)=b(i)-a(i,j)*x(j)
end do
x(i)=x(i)/A(i,i)
end do
end subroutine uptri
end module tri_eq
!############################################
module driver
contains
subroutine dri_main()
implicit real*8(a-z)
integer::ioerr !ioerr=0
character(20)::uptr
open(unit=11,file='fin.txt',iostat=ioerr)
open(unit=12,file='fout.txt')
do
read(11,*) uptr
!读输入文件
if(uptr(1:5)=='uptri') !提示这里有错误
call dri_up()
if(ioerr/=0) exit
end do
end subroutine dri_main
subroutine dri_up()
use tri_eq
implicit real*8(a-z)
integer,parameter::N=4
integer::i,j
real*8::A(N,N),b(N),x(N)
read(11,*)((A(i,j),j=1,N),i=1,N)
!读入b向量
read(11,*) b
call uptri(A,b,x,N)
write(12,101) x
101 format(T5,'上三角方程组的解',/,T4,'x='4(/f12.8))
end subroutine dri_up
end module driver
!################################################
program main
use driver
call dri_main()
end program main
!#################################################
其中fin.txt如下
uptri

      2  1  4  3
      0  3  2  5
      0  0  7  3
      0  0  0  2   

      50
      49
      53
      12
编译信息:
D:\fortran_programming\myfirst\main.f95:31.34:
          if(uptr(1:5)=='uptri')
                                           1
Error: Cannot assign to a named constant at (1)
D:\fortran_programming\myfirst\main.f95:51.8:
use driver
      1
Fatal Error: Can't open module file 'driver.mod' for reading at (1): No such file or directory
Process terminated with status 1 (0 minute(s), 0 second(s))
1 error(s), 0 warning(s) (0 minute(s), 0 second(s))
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

Whenever you feel like criticizing any one,just remember that all the people in this world haven’t had the advantages that you’ve had.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)

引用回帖:
3楼: Originally posted by physliebe at 2014-02-07 13:02:03
谢谢,是的。语法规定:if 如果不带 then 需要写成一行,如果写成两行,应该有 then 和 end if...

嗯。 if construct 和 if statement 的区别
4楼2014-02-07 20:40:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

btx97

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
physliebe: 金币+1, ★★★很有帮助 2014-02-07 20:01:04
if(uptr(1:5)=='uptri') !提示这里有错误
call dri_up()


后面没有 END IF?
2楼2014-02-07 19:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

physliebe

新虫 (小有名气)

引用回帖:
2楼: Originally posted by btx97 at 2014-02-07 19:27:13
if(uptr(1:5)=='uptri') !提示这里有错误
call dri_up()


后面没有 END IF?

谢谢,是的。语法规定:if 如果不带 then 需要写成一行,如果写成两行,应该有 then 和 end if
Whenever you feel like criticizing any one,just remember that all the people in this world haven’t had the advantages that you’ve had.
3楼2014-02-07 20:02:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 330分求调剂 +5 qzenlc 2026-03-29 5/250 2026-03-29 07:37 by 无际的草原
[考研] 289求调剂 +13 新时代材料 2026-03-27 13/650 2026-03-29 01:16 by 544594351
[考研] 085600,专业课化工原理,321分求调剂 +4 大馋小子 2026-03-28 4/200 2026-03-29 00:32 by 544594351
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 085701求调剂初试286分 +4 secret0328 2026-03-28 4/200 2026-03-28 21:09 by 15366876211
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 085404求调剂,总分309,本科经历较为丰富 +4 来财aa 2026-03-25 4/200 2026-03-28 07:41 by 棒棒球手
[考研] 340求调剂 +5 jhx777 2026-03-27 5/250 2026-03-28 04:18 by fmesaito
[考研] 322求调剂 +4 我真的很想学习 2026-03-23 4/200 2026-03-27 13:51 by 杨杨杨紫
[考研] 一志愿华东理工大学081700,初试分数271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by 惠州彭于晏
[考研] 调剂推荐 +5 清酒714 2026-03-26 6/300 2026-03-27 11:12 by 不吃魚的貓
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 一志愿北化求调剂 +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv风雷
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 8/400 2026-03-26 20:34 by Higraduate
[考研] 352求调剂 +4 大米饭! 2026-03-22 4/200 2026-03-26 16:40 by 不吃魚的貓
[考研] 086000生物与医药292求调剂 +6 小小陈小小 2026-03-22 9/450 2026-03-26 15:58 by dick_runner
[考研] 297求调剂 +6 田洪有 2026-03-26 6/300 2026-03-26 15:55 by 不吃魚的貓
[考研] 机械学硕总分317求调剂!!!! +4 Acaciad 2026-03-25 4/200 2026-03-25 19:59 by hanserlol
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
信息提示
请填处理意见