| 7 | 1/1 | 返回列表 |
| 查看: 1107 | 回復(fù): 6 | |||
nll1989新蟲 (著名寫手)
|
[交流]
想用Fortran語言實(shí)現(xiàn)矩陣連乘,并且矩陣?yán)锩嬗幸粋(gè)變量 已有3人參與
|
|
想用Fortran實(shí)現(xiàn)對(duì)矩陣M的連乘,而且矩陣?yán)锩婧幸粋(gè)變量EM,當(dāng)EM<=15,EM=6.7. 當(dāng)EM>15,EM=8.7,但是用以下程序不能實(shí)現(xiàn)。希望大家?guī)兔鉀Q,不勝感激! PROGRAM MAIN USE IMSL COMPLEX:: M(2,2),P(2,2) REAL*8 E,TE,TM,EM,TD,H,N COMPLEX*8 M11,M12,M21,M22,A1,A2,A3 TD=0.8D0 TM=1.0D0 EM=7.75D0 N=10 C====================================================== OPEN(UNIT=30,FILE="DNA-DAN01.DAT" ![]() DO E=6.89,9.750,0.001D0 C================================================ cc================================================ P=M DO I=1,N IF (I .LT. 15) THEN EM=6.7 ELSE IF (I .GE. 15) THEN EM=8.7 END IF C==================================================== H=acos((E-EM)/(2*TM)) M11=CMPLX((E-EM)/TD,0) M12=CMPLX(-1,0) M21=CMPLX(1,0) M22=CMPLX(0,0) M=RESHAPE ((/M11,M21,M12,M22/),(/2,2/)) P=MATMUL(M,P) END DO C===================================================== A1=4*SIN(H)*SIN(H) A2=(P(1,2)-P(2,1)+(P(1,1)-P(2,2))*COS(H))**2 A3=(P(1,1)+P(2,2))**2*SIN(H)**2 TE=A1/(A2+A3) Write(30,*) E,TE END DO END |
新蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
工程和科學(xué)數(shù)值計(jì)算咨詢
新蟲 (著名寫手)
新蟲 (初入文壇)
新蟲 (著名寫手)
銅蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 生物學(xué)調(diào)劑 +4 | Surekei 2026-03-21 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 |  ̄^ ̄゜汗 2026-03-19 | 4/200 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|