How to fit the curve to equality?

Question

I have a question on fitting to equality. Here is the code:

import numpy as np
from scipy.optimize import fsolve

phi_min_initial = 0.001
zc = 51
M = 25
x = n