shangxj

金蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 1246.7
散金: 97
帖子: 412
在線: 37.8小時
蟲號: 250453
注冊: 2006-05-13
性別: GG
專業(yè): 水泥與耐火材料

[求助] 如何用mathematica8解這一方程的所有整數(shù)解？

如題
方程為0.91x+0.45y+0.75z=70

請大俠不吝相助，并另有獎勵奉上。

回復此樓

» 猜你喜歡

哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有6人回復
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有91人回復
求調(diào)劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復
光學工程學碩調(diào)劑信息已經(jīng)有30人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有65人回復
華東師范大學芯片設計徐瓏真實水平咋樣已經(jīng)有1人回復
調(diào)劑已經(jīng)有0人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

誰來幫我一下好不好？mathematica 方程組輸入一直錯誤，不知道怎么辦已經(jīng)有3人回復
這樣一個偏微分方程有沒有解析解？已經(jīng)有4人回復
【軟件】Mathmatica 軟件 7.0版已經(jīng)有68人回復
【求助】如何解這樣的方程組【有效期至2011年4月30日】已經(jīng)有9人回復
【求助】用mathematica 5.0求解一個非線性方程組失敗，特發(fā)帖求助！已經(jīng)有5人回復
【求助】求教多元指數(shù)方程如何解有軟件否？已經(jīng)有6人回復

1樓 2012-06-01 14:53:39

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

shangxj

金蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 1246.7
散金: 97
帖子: 412
在線: 37.8小時
蟲號: 250453
注冊: 2006-05-13
性別: GG
專業(yè): 水泥與耐火材料

還沒有達人看到嗎？
還是根本就無法解？

贊一下

回復此樓

2樓2012-06-04 12:58:02

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

★
感謝參與，應助指數(shù) +1
zhangguangping: 金幣+1, 謝謝指教！ 2012-06-04 20:37:11

有無窮多組解吧！
有無窮多組解吧！

贊一下

回復此樓

3樓2012-06-04 14:43:00

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangguangping

木蟲 (著名寫手)

計算強帖: 6
應助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時
蟲號: 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
shangxj: 金幣+10, ★★★很有幫助, 編程解決問題真是好，謝謝回應。 2012-06-06 17:50:45

雖然你要求用mathematica8解這個方程的整數(shù)解。但是我不會mathematica8這個軟件。我用Fortran編寫了一段小程序，給你解了一下。為了減小計算量，首先解了一下非負整數(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]這個區(qū)間內(nèi)的整數(shù)就有33177組。估計你的這個問題也是有實際應用的，所以我覺得你肯定是有區(qū)間界定的。因此我把我的代碼放在這兒，你可以通過調(diào)整這個區(qū)間的上下線來尋找你想要的結果：

   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(*,'(a4,I5,2x,a4,I5,2x,a4,I5)')"x = ",i,"y = ",j,"z = ",k
   write(100,'(a4,I5,2x,a4,I5,2x,a4,I5)')"x = ",i,"y = ",j,"z = ",k
   endif

      enddo
            enddo
         enddo

      end program sovl

贊一下

回復此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權問題，其責任在于用戶本人，如對版權有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : solv2.f

2012-06-04 21:48:57, 631 bytes

弘德明志博學篤行

4樓2012-06-04 21:50:20

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangguangping

木蟲 (著名寫手)

計算強帖: 6
應助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時
蟲號: 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
shangxj: 金幣+10, ★★★很有幫助, 寫程序解決問題對我來說有點難度，我還是想使用簡單的工具解決相關的一類問題。只是平時雜事頗多，靜不下心來深入去學。 2012-06-06 17:53:38
csgt0: 金幣+2, 多謝應助 2012-09-19 09:46:09

考慮到你可能不能編譯fortran程序，我將fortran程序做成可執(zhí)行文件了。你只需要下載那個exe按照提示輸入你的求解范圍，就可以直接將結果輸出到文件中去了。

希望對你有用。

贊一下(2人)

回復此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權問題，其責任在于用戶本人，如對版權有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : solv2.f

2012-06-04 22:10:38, 1.6 K

附件 2 : solv2.exe

2012-06-04 22:10:43, 344.03 K

弘德明志博學篤行

5樓2012-06-04 22:10:50

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

gotit

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3.5
帖子: 2
在線: 5分鐘
蟲號: 2011522
注冊: 2012-09-18

【答案】應助回帖

★
csgt0: 金幣+1, 多謝應助 2012-09-19 09:46:17

Solve[0.91 x + 0.45 y + 0.75 z == 70, {x, y, z}, Integers]

{{z -> ConditionalExpression[
93.3333 - 1.21333 x - 0.6 y, (x | y | z) \[Element] Integers]}}

贊一下

回復此樓

6樓2012-09-18 21:53:23

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

shangxj

金蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 1246.7
散金: 97
帖子: 412
在線: 37.8小時
蟲號: 250453
注冊: 2006-05-13
性別: GG
專業(yè): 水泥與耐火材料

引用回帖:

6樓: Originally posted by gotit at 2012-09-18 21:53:23
Solve

{{z -> ConditionalExpression}}

怎么沒有出來結果？
out顯示還是跟輸入的形式一樣，是我操作不對嗎？

贊一下

回復此樓

7樓2012-10-08 17:55:29

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 317求調(diào)劑 +4	申子申申 2026-03-19	9/450	2026-03-20 11:08 by 申子申申
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 321求調(diào)劑 +8	何潤采123 2026-03-18	10/500	2026-03-19 16:46 by 何潤采123
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 0703化學調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 本科鄭州大學物理學院，一志愿華科070200學碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 0703化學 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 一志愿華中科技大學，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 311求調(diào)劑 +6	26研0 2026-03-15	6/300	2026-03-18 14:43 by haxia
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 生物學071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 0703化學336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 293求調(diào)劑 +11	zjl的號 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu