To study this type of system, you need to build an equivalent small-signal (linear) version of it. The simplest one you can think of is below:

The operating point looks good but I used a perfect op amp to drive a largely-capacitive load such as a big MOSFET as you selected. Also, I used a simplified hybrid model for the MOSFET as you can see. So this sketch is mainly for illustration purposes. You can improve the analysis and use a more realistic model for the op amp and the MOSFET. As I type, I did not have the exact model your circuit uses but you see the spirit in the below picture:

The loop is physically open via the LoL/CoL elements inserted in series with the sense resistance. When SPICE computes its dc operating point, LoL is shorted and it can determine the bias point. Then, as the ac analysis starts, CoL allows the ac to pass while LoL blocks the return: the loop is closed in dc but open in ac. If you plot \$V_{out}\$ over \$V_{in}\$, you will have \$T(s)\$ the open-loop gain and see what compensation is needed on the op amp to stabilize the whole thing. If you want to run a symbolic analysis, replace the op amp by its macro model and the MOSFET by a more realistic small-signal model. I believe however that a SPICE analysis should already provide good insight. Good luck!