pub fn solve_22(a: Mat22, b: Vec2) -> Vec2
Solve A * x = b, where b is a column vector. This is more efficient than computing the inverse in one-shot cases.