% See also kfinit, kffk, kfhk, kfplot. subplot(321), plot(t, [res(:,1:2),xkpk(:,1:2)]/glv.deg), xygo('phixy'); subplot(323), plot(t, [res(:,3),xkpk(:,3)]/glv.deg ...