| 查看: 752 | 回复: 3 | ||
[求助]
Fortran中的一个高斯程序已有1人参与
|
|
各个大神门 有没有见过下面这个程序的 可以给解答下这是求什么的吗 ？编写的什么公式？？？ subroutine gauss(n,ns,BB) implicit none integer i,j,k,n1,nt,n,ns,ip,i1 double precision BB(2,3),x,AA(3,4) n1=n+1 nt=n+ns do i=1,2 do j=1,3 aa(i+1,j+1)=bb(i,j) enddo enddo if (n.gt.1) then do 555 i=2,n ip=i i1=ip x=abs(aa(i1,i1)) do 444 j=i+1,n+1 if (abs(aa(j,i1)).ge.x) then x=abs(aa(j,i1)) ip=j endif 444 continue if (ip.ne.I1) then do 333 j=i1,nt+1 x=aa(i1,j) aa(i1,j)=aa(ip,j) aa(ip,j)=x 333 continue endif do 222 j=i+1,n+1 if (abs(aa(i1,i1)).lt.1.d-10) then aa(i1,i1)=1.0 endif x=aa(j,i1)/aa(i1,i1) do 111 k=i+1,nt+1 aa(j,k)=aa(j,k)-x*aa(i1,k) 111 continue 222 continue 555 continue endif do 50 ip=1,n i=n1-ip+1 do 40 k=n1+1,nt+1 if (abs(aa(i,i)).lt.1.d-10) then aa(i,i)=1.0 endif aa(i,k)=aa(i,k)/aa(i,i) if (i.ne.2) then i1=i-1 do 30 j=2,i1 aa(j,k)=aa(j,k)-aa(i,k)*aa(j,i) 30 continue endif 40 continue 50 continue do i=1,n do j=1,n+1 bb(i,j)=aa(i+1,j+1) enddo enddo end |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有148人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 注册化工工程师须知
已经有23人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
小丫鸭
2楼2015-11-04 15:50:01
非常感谢你的回复 问过老师 老师也不确定说是求逆 还是不太懂3楼2015-11-04 16:56:51
4楼2015-11-04 16:57:49