hold on
grid
plot(xrr,4*phir,':')
%plot(xaa,phia,'--b')
%plot(xx98,phi98,'-+')
%plot(xtt,phit,'.-r')
plot(xb,phib)
hold off