x=linspace(0,0.45,20); a=VaporPressureLiCl(15,x); b=VaporPressureLiCl(30,x); c=VaporPressureLiCl(45,x); d=VaporPressureLiCl(60,x); plot(x,a,'-+',x,b,'-^',x,c,'-s',x,d,'-d') xlabel('Mass Fraction [kg_s_a_l_t/kg_w_a_t_e_r]','FontSize',24) ylabel('Equilibrium Vapor Pressure [kPa]','FontSize',24) set(gca,'FontSize',20) h=legend('T=15\circC','T=30\circC','T=45\circC','T=60\circC'); set(h,'FontSize',18) grid on