| 查看: 1199 | 回复: 3 | |||
[求助]
找不出毛病不知怎么改？ 已有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)) |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Materials Letters投稿一周后返回技术审查意见，不知道怎么修改，求高手指教
已经有8人回复
- 博士做了一年的实验感觉自己试验设计出问题，都不知道能不能毕业了
已经有18人回复
- castrp计算出错，不知什么原因，求大神解答
已经有12人回复
- lammps实验拉伸Al纳米线的in文件的错误提示好久找不出？
已经有8人回复
- 论文提交时图形不符合要求，但不知怎么修改
已经有8人回复
- 语言修改。请大家帮忙分析comments的意思。
已经有8人回复
- 导师不给改论文，怎么办？
已经有18人回复
- 自己写的一个小程序，出现问题不知道怎么解决
已经有9人回复
- 2012年5月31日基金评审体会 (转载)
已经有27人回复
- 唱歌找不着调也是病 被称为“失歌症”
已经有27人回复
- 投稿时审稿人给的关于XPS的修改意见，不知道该怎么处理了，求办法…………
已经有3人回复
- 采用matlab中nlinfit拟合，出现问题了，不知道是哪错了，该怎么解决，求高手指点一下
已经有9人回复
- chemcomm大修，第一次投文章，一些问题不知道怎么回答，虫子们帮帮忙啊！
已经有17人回复
- 英文写出来了，老师没提意见，论文不知如何修改，不敢投啊？
已经有6人回复
- 【求助/交流】关于重组质粒的构建，总是失败，希望高手指点一下啊～～
已经有16人回复
- 论文大修，问题一大堆，修改意见四页多，修后录用可能性大不？
已经有30人回复
- 【求助】有很多错误的程序，找不出来了
已经有18人回复
- 在企业里面研发做不出产品怎么办
已经有29人回复
2楼2014-02-07 19:27:13
3楼2014-02-07 20:02:03
4楼2014-02-07 20:40:04
