Exactly solvable inhomogeneous fermion systems

15 exactly solvable inhomogeneous (spinless) fermion systems on one-dimensional lattices are constructed explicitly based on the discrete orthogonal polynomials of Askey scheme, e.g. the Krawtchouk, Hahn, Racah, Meixner, q-Racah polynomials. The Schrödinger and Heisenberg equations are solved explicitly, as the entire set of the eigenvalues and eigenstates are known explicitly. The ground state two point correlation functions are derived explicitly. The multi point correlation functions are obtained by Wick's Theorem. Corresponding 15 exactly solvable XX spin systems are also displayed. They all have nearest neighbour interactions. The exact solvability of Schrödinger equation means that of the corresponding Fokker-Planck equation. This leads to 15 exactly solvable Birth and Death fermions and 15 Birth and Death spin models. These provide plenty of materials for calculating interesting quantities, e.g. entanglement entropy, etc.