Why doesn't SymPy calculate the system of equations?

Question

Python

import sympy
a,b = sympy.symbols('a,b')
equations = [sympy.Eq((9*a+b)%32,14), sympy.Eq((15*a+b)%32, 5)]
sympy.solve(equations)
Traceback (most recent call last):
File "", line 1, in
File "/usr/local/lib/python3.4/site-packages/sympy/solvers/solvers.py" , line 911, in solve
solution = _solve_system(f, symbols, **flags)
File "/usr/local/lib/python3.4/site-packages/sympy/solvers/solvers.py", line 1646, in _solve_system
raise NotImplementedError('could not solve %s' % eq2)
NotImplementedError: could not solve Mod(9*a + b, 32) - 14

What am I doing wrong?

Answer 1

Well, these are Diophantine equations, as far as I understand. There is a separate module for them .