| 查看: 497 | 回复: 2 | |||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | |||
[求助]
fortran编程
|
|||
| 求大神帮忙解答一下：fortran中如何将一个二维数组求行列式，然后将该行列式的值作为另一个数组的元素 |
回复此楼» 猜你喜欢
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有3人回复
-
需要合成515-64-0，50g，能接单的留言
已经有3人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有10人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
-
26申博自荐
已经有7人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有4人回复
-
带资进组求博导收留
已经有9人回复
-
求助大佬们，伤口沾上了乙腈
已经有6人回复
2楼2016-10-31 23:46:00
|
do i=1,nele,1 PP(1:3)=(/nodx(nod(i,1)), nody(nod(i,1)),nodz(nod(i,1)) /) QQ(1:3)=(/nodx(nod(i,2)), nody(nod(i,2)),nodz(nod(i,2)) /) RR(1:3)=(/nodx(nod(i,3)), nody(nod(i,3)),nodz(nod(i,3)) /) SS(1:3)=(/nodx(nod(i,4)), nody(nod(i,4)),nodz(nod(i,4)) /) a1(1:3,1:3)=(/QQ(1),QQ(2),QQ(3),RR(1),RR(2),RR(3),SS(1),SS(2),SS(3)/) b1(1:3,1:3)=-(/1.0,QQ(2),QQ(3),1.0,RR(2),RR(3),1.0,SS(2),SS(3)/) c1(1:3,1:3)=(/1.0,QQ(1),QQ(3),1.0,RR(1),RR(3),1.0,SS(1),SS(3)/) d1(1:3,1:3)=-(/1.0,QQ(1),QQ(2),1.0,RR(1),RR(2),1.0,SS(1),SS(2)/) aae1(i)=det(a1(1:3,1:3)) bbe1(i)=det(b1(1:3,1:3)) cce1(i)=det(c1(1:3,1:3)) dde1(i)=det(d1(1:3,1:3)) end do （这是我程序中的一部分，数组都以正确定义，但运行时a1,b1,c1,d1部分总会出现 Error: The shapes of the array expressions do not conform. 是什么原因啊？） |
3楼2016-11-01 09:15:53