lmax=5.8;a=.804;
NDSolve[{PH'[l]==X[l]^2/a^3,PH[0]==0.,AR'[l]==Y[l] Cos[PH[l]],Y'[l]==Sin[PH[l]],Y[0]==0.,X'[l]==Cos[PH[l]],X[0]==0.,AR[0]==0},{PH,X,Y,AR},{l,0,lmax}];
{ph[u_],x[u_],y[u_],ar[u_]}={PH[u],X[u],Y[u],AR[u]}/.First[%];
ParametricPlot[{x[l],y[l]},{l,0,lmax}];
Plot[Tooltip[{ph[l],ar[l], x [l],y[l] }],{l,0,lmax}]
I prefer the Tootip to read labels ph,ar, x,y as named but not again echo the same "Interpolating Function" [...]