hold on
grid
plot(xrr,phir2,':y')
plot(xaa,phia2,'--g')
plot(xx98,phi982,'-+')
plot(xtt,phit2,'.-r')
plot(xb,phib)
hold off