Re: Re: Diophantic equations

*To*: mathgroup at smc.vnet.net*Subject*: [mg32788] Re: [mg32771] Re: [mg32746] Diophantic equations*From*: Andrzej Kozlowski <andrzej at tuins.ac.jp>*Date*: Sat, 9 Feb 2002 23:39:34 -0500 (EST)*Sender*: owner-wri-mathgroup at wolfram.com

In my original reply I failed to notice you only wanted non-negative solutions. One can of course obtain them from the solution below by trying various values for m, but if you prefer to do it more automatically here is one way: First load the InequalitySolve package: In[1]:= << "Algebra`InequalitySolve`" since a and b have to be >=0 we must have: In[2]:= InequalitySolve[{6300 + m*19 >= 0, -16200 - m*49 >= 0}, {m}] Out[2]= -(6300/19) <= m <= -(16200/49) In other words: In[3]:= N[%] Out[3]= -331.57894736842104 <= m <= -330.61224489795916 There is only one such integer In[4]:= m = -331; thus In[5]:= a = 6300 + m*19 Out[5]= 11 In[6]:= b = -16200 - m*49 Out[6]= 19 In[7]:= c = 100 - 5*a - 2*b Out[7]= 7 Andrzej Kozlowski Toyama International University JAPAN http://platon.c.u-tokyo.ac.jp/andrzej/ On Saturday, February 9, 2002, at 07:11 PM, Andrzej Kozlowski wrote: > Here is how to do it in your case. You can use this method to solve more > general linear cases. > > First eliminate one variable from the equations: > > In[1]:= > eqs={a+b+10 c==100,5 a+2 b+c==100} > > In[2]:= > Eliminate[eqs,{c}] > > Out[2]= > 900-19 b==49 a > > Clearly we only need to find all integer solutions of the equation 49 > a + 19 b == 900, since c will then be given by c=100-a-2b. > > The next step is to find just a single integer solution. The general way > to do this is to use the Euclidean algorithm. > > In[3]:= > ExtendedGCD[49,19] > > Out[3]= > {1,{7,-18}} > > This tells us that: > > In[4]:= > 7*49-18*19 > > Out[4]= > 1 > > From this we see at once that a solution will be given by: > > In[7]:= > a=900*7 > > Out[7]= > 6300 > > In[8]:= > b=900*-18 > > Out[8]= > -16200 > > indeed > > In[9]:= > 49 a+19 b > > Out[9]= > 900 > > Now that we have one solution (a,b) all the other solutions are given by > (a+m 19,b-m 49) = (6300+ m*19,-16200-m*49) where m is any integer. > Indeed: > > In[11]:= > {6300+m*19,-16200-m*49}.{49,19}//Simplify > > Out[11]= > 900 > > > > On Friday, February 8, 2002, at 05:49 PM, Max Ulbrich wrote: > >> Hi, >> >> I try to solve equations where only integers are allowed, like >> >> a + b + 10 c = 100 >> 5 a + 2 b + c = 100 >> >> a, b, c should be integers and positve. >> >> How can I do this? >> >> Max >> >> >> >> >> > > >