| 查看: 1026 | 回复: 0 | ||
[求助]
求助：ludcmp&lubksb
|
|
subroutine ludcmp(coeff, numvar, np, indx, d) integer :: numvar, np, indx(numvar) doubleprecision :: d, coeff(np,np) integer, parameter :: nmax=500 doubleprecision, parameter :: tiny=1.0d-20 integer :: i, imax, j, k doubleprecision :: aamax, dum, sum, vv(nmax) d=1 do i = 1, numvar aamax = 0.d0 do j = 1, numvar if (abs(coeff(i,j)).gt.aamax) aamax = abs(coeff(i,j)) end do !if (aamax .eq. 0.d0) pause 'singular matrix in ludcmp' vv(i) = 1./aamax end do do j = 1, numvar do i= 1, j-1 sum = coeff(i,j) do k=1,i-1 sum = sum-coeff(i,k)*coeff(k,j) end do coeff(i,j) = sum end do aamax = 0.d0 do i = j,numvar sum = coeff(i,j) do k=1, j-1 sum = sum - coeff(i,k) * coeff(k,j) end do coeff(i,j) = sum dum = vv(i) * abs(sum) if (dum .ge. aamax) then imax = i aamax = dum end if end do if (j .ne. imax) then do k = 1, numvar dum = coeff(imax,k) coeff(imax,k) = coeff(j,k) coeff(j,k) = dum end do d = -d vv(imax) = vv(j) end if indx(j) = imax if (coeff(j,j) .eq. 0.d0) coeff(j,j) = tiny if (j .ne. numvar) then dum = 1./coeff(j,j) do i = j+1, numvar coeff(i,j) = coeff(i,j) * dum end do end if end do return end subroutine ludcmp 请教，这个的功能是怎样的，谢谢！ |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
【求助】除去2-溴丙酸乙酯
已经有3人回复
- 求助啊，SCI修稿后一月有余，询问编辑，给了一封回信，大家帮忙看一下是什么意思啊
已经有24人回复
- 求助：气相色谱测定醇醚平行性差如何解决
已经有8人回复
- 【求助】请教合成CdTe@SiO2
已经有6人回复
- 【求助/交流】NCBI给了我登录号但是用号查不到，这样算不算登录成功呢？
已经有8人回复
- 求助：Materials Chemistry and Physics的word 模板
已经有5人回复
- 【求助】分离黄酮和酚酸
已经有12人回复
- 【求助】羧酸的胺盐转化为酰胺
已经有5人回复
- 【求助】XRD几个小问题求助
已经有3人回复
- 【求助】请求帮助：测得一个晶体数据有这么多错误！怎么办？
已经有4人回复
- 【求助/交流】克隆Amp/Kan选择以及蓝白筛选
已经有10人回复
- 【求助】关于顺磁性物质和抗磁性物质
已经有12人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人