Добавлю маленько. Если знать, как решается уравнение с двумя переменными, то увеличение числа переменных принципиальных трудностей не создаёт. Подробности можно поискать в известных учебниках по теории чисел. Вычленить из Бухштаба вычислительный алгоритм будет для начинающего несколько трудноват. Также не искал есть ли в инете нормальные сайты с примерами. Изложу одну из возможных схем, ограничиваясь этим примером.
Итак, требуется найти все целые решения уравнения
Временно считаем
![$z$ $z$](https://dxdy-04.korotkov.co.uk/f/f/9/3/f93ce33e511096ed626b4719d50f17d282.png)
ио параметра, то есть рассматриваем уравнение
![$15x+21y=8-35z, $ $15x+21y=8-35z, $](https://dxdy-02.korotkov.co.uk/f/d/f/2/df2383accf5aa181143e5fe45aea1d8882.png)
считая параметр
![$z$ $z$](https://dxdy-04.korotkov.co.uk/f/f/9/3/f93ce33e511096ed626b4719d50f17d282.png)
известным.
Так как левая часть делится на 3, то это уравнение явно неразрешимо, если
![$8-35z$ $8-35z$](https://dxdy-01.korotkov.co.uk/f/c/1/c/c1cfeb425779836d7873bb946bf28be282.png)
не делится на 3, иначе же
![$z\equiv 1 \pmod{3}$ $z\equiv 1 \pmod{3}$](https://dxdy-01.korotkov.co.uk/f/4/7/4/4749328b934415944d1b364f7c61e30b82.png)
, то есть
![$z=1+3u, \, u\in \mathbb Z$ $z=1+3u, \, u\in \mathbb Z$](https://dxdy-01.korotkov.co.uk/f/0/5/9/059f9f16938be5b5fc30f538f0df6d1e82.png)
.
Подставив в уравнение и сократив на 3, получим
![$5x+7y=-9-35u.$ $5x+7y=-9-35u.$](https://dxdy-02.korotkov.co.uk/f/9/9/b/99b9227df32f66c9b79b0cd98ca135f182.png)
Это уже стандартное уравнение общим решением:
![$$\begin{pmatrix}x\\ y\end{pmatrix}=\begin{pmatrix}-3\\ 2\end{pmatrix}(9+35u)+\begin{pmatrix}7\\ -5\end{pmatrix}v,$$ $$\begin{pmatrix}x\\ y\end{pmatrix}=\begin{pmatrix}-3\\ 2\end{pmatrix}(9+35u)+\begin{pmatrix}7\\ -5\end{pmatrix}v,$$](https://dxdy-01.korotkov.co.uk/f/4/9/e/49e8b880776ca2b7dc1541d1e4d2493a82.png)
Как итог переменные
![$x,y,z$ $x,y,z$](https://dxdy-03.korotkov.co.uk/f/2/4/4/244be3c7db382d3e1400c7c4caa1023a82.png)
выражены через произвольные целочисленные параметры
![$u$ $u$](https://dxdy-03.korotkov.co.uk/f/6/d/b/6dbb78540bd76da3f1625782d42d6d1682.png)
и
![$v$ $v$](https://dxdy-03.korotkov.co.uk/f/6/c/4/6c4adbc36120d62b98deef2a20d5d30382.png)
.
В общем случае, линейное уравнение с
![$n$ $n$](https://dxdy-02.korotkov.co.uk/f/5/5/a/55a049b8f161ae7cfeb0197d75aff96782.png)
неизвестными будет разрешимо только если правая часть делится на НОД коэффициентов левой части и общее решение в таком случае будет содержать
![$n-1$ $n-1$](https://dxdy-03.korotkov.co.uk/f/e/f/c/efcf8d472ecdd2ea56d727b5746100e382.png)
произвольных целочисленных параметров. Схема та же, как показано - одно из переменных ио параметра (чтобы заслужить такую честь, надо иметь соответствующий вид - что-то на что-то должно делиться), получаем уравнение с параметром от меньшего числа неизвестных, опять уменьшаем число неизвестных и т.д.