| 查看: 1031 | 回复: 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 请教，这个的功能是怎样的，谢谢！ |
回复此楼» 猜你喜欢
为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有11人回复
-
什么是人一生最重要的？
已经有4人回复
-
版面费该交吗
已经有3人回复
-
今年春晚有几个节目很不错，点赞！
已经有12人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
找到一些相关的精华帖子，希望有用哦~
【求助】除去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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人