| 5 | 1/1 | 返回列表 |
| 查看: 1239 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
shangxj金蟲 (正式寫手)
|
[求助]
如何用mathematica8解這一方程的所有整數(shù)解?
|
|
|
如題 方程為0.91x+0.45y+0.75z=70 請(qǐng)大俠不吝相助,并另有獎(jiǎng)勵(lì)奉上。 |
鐵桿木蟲 (職業(yè)作家)
金蟲 (正式寫手)
木蟲 (著名寫手)
|
雖然你要求用mathematica8解這個(gè)方程的整數(shù)解。但是我不會(huì)mathematica8這個(gè)軟件。我用Fortran編寫了一段小程序,給你解了一下。為了減小計(jì)算量,首先解了一下非負(fù)整數(shù)解如下: x = 10 y = 2 z = 80 x = 10 y = 7 z = 77 x = 10 y = 12 z = 74 x = 10 y = 17 z = 71 x = 10 y = 22 z = 68 x = 10 y = 27 z = 65 x = 10 y = 32 z = 62 x = 10 y = 37 z = 59 x = 10 y = 42 z = 56 x = 10 y = 47 z = 53 x = 10 y = 52 z = 50 x = 10 y = 57 z = 47 x = 10 y = 62 z = 44 x = 10 y = 67 z = 41 x = 10 y = 72 z = 38 x = 10 y = 77 z = 35 x = 10 y = 82 z = 32 x = 10 y = 87 z = 29 x = 10 y = 92 z = 26 x = 10 y = 97 z = 23 x = 10 y = 102 z = 20 x = 10 y = 107 z = 17 x = 10 y = 112 z = 14 x = 10 y = 117 z = 11 x = 10 y = 122 z = 8 x = 10 y = 127 z = 5 x = 10 y = 132 z = 2 x = 25 y = 0 z = 63 x = 25 y = 5 z = 60 x = 25 y = 10 z = 57 x = 25 y = 15 z = 54 x = 25 y = 20 z = 51 x = 25 y = 25 z = 48 x = 25 y = 30 z = 45 x = 25 y = 35 z = 42 x = 25 y = 40 z = 39 x = 25 y = 45 z = 36 x = 25 y = 50 z = 33 x = 25 y = 55 z = 30 x = 25 y = 60 z = 27 x = 25 y = 65 z = 24 x = 25 y = 70 z = 21 x = 25 y = 75 z = 18 x = 25 y = 80 z = 15 x = 25 y = 85 z = 12 x = 25 y = 90 z = 9 x = 25 y = 95 z = 6 x = 25 y = 100 z = 3 x = 25 y = 105 z = 0 x = 40 y = 3 z = 43 x = 40 y = 8 z = 40 x = 40 y = 13 z = 37 x = 40 y = 18 z = 34 x = 40 y = 23 z = 31 x = 40 y = 28 z = 28 x = 40 y = 33 z = 25 x = 40 y = 38 z = 22 x = 40 y = 43 z = 19 x = 40 y = 48 z = 16 x = 40 y = 53 z = 13 x = 40 y = 58 z = 10 x = 40 y = 63 z = 7 x = 40 y = 68 z = 4 x = 40 y = 73 z = 1 x = 55 y = 1 z = 26 x = 55 y = 6 z = 23 x = 55 y = 11 z = 20 x = 55 y = 16 z = 17 x = 55 y = 21 z = 14 x = 55 y = 26 z = 11 x = 55 y = 31 z = 8 x = 55 y = 36 z = 5 x = 55 y = 41 z = 2 x = 70 y = 4 z = 6 x = 70 y = 9 z = 3 x = 70 y = 14 z = 0 在[-1000,1000]這個(gè)區(qū)間內(nèi)的整數(shù)就有33177組。估計(jì)你的這個(gè)問題也是有實(shí)際應(yīng)用的,所以我覺得你肯定是有區(qū)間界定的。因此我把我的代碼放在這兒,你可以通過調(diào)整這個(gè)區(qū)間的上下線來尋找你想要的結(jié)果: program sovl integer*4 i, j, k integer*4 max_x, max_y, max_z real*8 resu,EPS parameter (EPS=1.d-300) max_x=int(70d0/0.91d0)+1 max_y=int(70d0/0.45d0)+1 max_z=int(70d0/0.75d0)+1 open(unit=100,file="solv.dat" ![]() do i=-1000,1000,1 do j=-1000,1000,1 do k=-1000,1000,1 resu=(real(i)*0.91d0+real(j)*0.45d0+real(k)*0.75d0) if(abs(resu-70d0) write(100,'(a4,I5,2x,a4,I5,2x,a4,I5)')"x = ",i,"y = ",j,"z = ",k endif enddo enddo enddo end program sovl |

木蟲 (著名寫手)
|
考慮到你可能不能編譯fortran程序,我將fortran程序做成可執(zhí)行文件了。你只需要下載那個(gè)exe按照提示輸入你的求解范圍,就可以直接將結(jié)果輸出到文件中去了。 希望對(duì)你有用。 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎? +6 | redcom 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 283求調(diào)劑(080500) +13 | A child 2026-03-27 | 13/650 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +9 | Nadiums 2026-03-25 | 10/500 |
|
|
[考研] 317求調(diào)劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 291求調(diào)劑 +7 | Y-cap 2026-03-29 | 8/400 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|