| 查看: 1035 | 回复: 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 请教，这个的功能是怎样的，谢谢！ |
回复此楼» 猜你喜欢
307求调剂
已经有3人回复
-
085600材料与化工专硕329 求调剂
已经有14人回复
-
材料考研求调剂总分280
已经有17人回复
-
085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干
已经有4人回复
-
265求调剂
已经有15人回复
-
化工求调剂！
已经有20人回复
-
266求调剂
已经有22人回复
-
化工调剂303分，过四级
已经有29人回复
-
生物学363调剂求助
已经有4人回复
-
22408 一志愿双一流人工智能300分 四六级，数据分析国奖
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
【求助】除去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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
