Please recreate Figure 9.7 from your text. Make the figures look tidy, and include a vector field and nullclines. Make up your own parameter values but also look for some guidance within the book if you don’t know where to start.

If you want to create a multi-panel plot, then use the function par() and set the argument mfrow = c(n, m) with n being the number of rows and m being the number of columns (e.g., par(mfrow = c(2, 3)) creates a plot with 2 panels for rows and 3 columns).
There is often a lot of white space as a default around plot(). You can control that with mar = c(a, b, c, d)) as an argument in par(). ad respectfully adjust the whitespace at the bottom, left, top, and right of each panel.

It might be useful to make the graphs that go in the first column first (they use linear interaction terms) then the second column (they have the nonlinear interaction terms).