oxygen_solve.m