Iceburn Stars burn subplot(2,3,3) plot(r/ro,log(rho),'g','LineWidth',2); set(gca(),"grid",[2 2]) title('Density'); ('square')