Step 1. Take a derivative of your function 1 (F1) and 2 (F2) and save them dF1 and dF2. You can use syms x and diff() here.
Step 2. Compute their values (F1(x), F2(x), dF1(x), dF2(x)) at x (e.g.: x = -pi... pi). You can use subs() or just plugin the values of x into the expressions of F1, F2, dF1, dF2.
Step 3. Plot the computed values: F1 vs. dF1 in Figure 1 and F2 vs. dF2 in Figure 2. This is a phase plot.
