/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created by wxMaxima version 0.7.4 ] */

/* [wxMaxima: input   start ] */
Cp_ice(T):=0.505+0.0018*T;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Cp_water(T):=4.216-0.0026*T;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Cp_steam(T):=0.991+T*1e-4;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
H_ice(T):=integrate(Cp_ice(T),T,173,T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
H_ice(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
T:273$ H_a:%o5$ remvalue(T)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
H_water(T):=H_a+334+integrate(Cp_water(T),T,273,T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
H_water(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
T:373$ H_b:%o10$ remvalue(T)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
H_steam(T):=H_b+2272+integrate(Cp_water(T),T,373,T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([H_ice(T),H_water(T),H_steam(T)], [T,173,373])$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
H_steam(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
H(T):=if T<273 then %o5 elseif T<373 then %o10 else %o16;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([H(T)],[T,173,473]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
plot2d([H(T)],[T,173,473],[y,-500,4000],[xlabel,"Temperature (K)"],[ylabel,"Specific Enthalpy (J/g)"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S_ice(T):=integrate(Cp_ice(T)/T,T,173,T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S_ice(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
T:273$ H_a:(2525*log(T)+9*T)/5000-(2525*log(173)+1557)/5000$ remvalue(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S_water(T):=H_a+334/273+integrate(Cp_water(T)/T,T,273,T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S_water(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
T:373$ H_b:(21080*log(T)-13*T)/5000-(21080*log(273)-3549)/5000+(2525*log(273)+2457)/5000-(2525*log(173)+1557)/5000+334/273$ remvalue(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S_steam(T):=H_b+2272/373+integrate(Cp_steam(T)/T,T,373,T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S_steam(T);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([S_ice(T),S_water(T),S_steam(T)], [T,173,373])$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S(T):=if T<273 then (2525*log(T)+9*T)/5000-(2525*log(173)+1557)/5000 elseif T<373 then (21080*log(T)-13*T)/5000-(21080*log(273)-3549)/5000+(2525*log(273)+2457)/5000-(2525*log(173)+1557)/5000+334/273 else (9910*log(T)+T)/10000+(21080*log(373)-4849)/5000-(9910*log(373)+373)/10000-(21080*log(273)-3549)/5000+(2525*log(273)+2457)/5000+(-2525*log(173)-1557)/5000+744838/101829;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([S(T)],[T,173,473]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
plot2d([S(T)],[T,173,473],[y,-2,12],[xlabel,"Temperature (K)"],[ylabel,"Specific Entropy
  (J/g K)"]);
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$