Here you go:
phsr1 = [4.008 59.935; 7.086 27.616; 7.75 14.37];
phsr2 = [5.99 -30.07; 3.716 -62.32; 1.946 -75.92];
figure(1)
compass(phsr1(:,1).*cosd(phsr1(:,2)), phsr1(:,1).*sind(phsr1(:,2)))
figure(2)
compass(phsr2(:,1).*cosd(phsr2(:,2)), phsr2(:,1).*sind(phsr2(:,2)))
Since this is likely a homework assignment, I’ll let you figure out how it works. The important information are in the documentation for the various functions.
