| 查看: 971 | 回复: 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 |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有74人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 注册化工工程师须知
已经有23人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
小丫鸭
2楼2015-11-04 15:50:01
非常感谢你的回复 问过老师 老师也不确定说是求逆 还是不太懂 3楼2015-11-04 16:56:51
4楼2015-11-04 16:57:49