Finding the fixed point and nullclines for five dimensional system?

I have a 5-dimensional model, where each dimension's evolution is given by an ordinary differential equation, and I am looking to find its nullclines and fixed points. What would be the most efficient way to go about this problem in MATLAB?

Thanks for your help!