24小時熱門版塊排行榜

>論壇更新日志 (6452)
>導師招生 (2520)
>考研 (2385)
>蟲友互識 (516)
>文獻求助 (453)
>考博 (348)
>碩博家園 (265)
>招聘信息布告欄 (133)
>論文投稿 (116)
>休閑灌水 (116)
>博后之家 (88)
>基金申請 (84)
>公派出國 (50)
>綠色求助(高懸賞) (43)
>找工作 (42)
>教師之家 (41)

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

shangxj

金蟲 (正式寫手)

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

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

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

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

回復此樓

» 猜你喜歡

面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有268人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！已經(jīng)有1人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8 已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有2人回復
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院已經(jīng)有9人回復

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

誰來幫我一下好不好？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的回帖

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ū)間的上下線來尋找你想要的結(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(*,'(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

贊一下

回復此樓