semicolons suppress output

Extras/LaTeX/tables/RewriteDerivatives.m

@@ -13,16 +13,19 @@
 syms d2P_D2 d2P_T2 d2P_TD d2U_T2;
 
 % step 1: rewrite the first-order derivatives (wrt d and T)
-dU_D = -T/D^2*dP_T+P/D^2
-dS_T = 1/T*dU_T
-dS_D = -1/D^2*dP_T
-dH_T = dU_T+1/D*dP_T
-dH_D = 1/D*dP_D - T/D^2*dP_T
+dU_D = -T/D^2*dP_T+P/D^2 ;
+dS_T = 1/T*dU_T ;
+dS_D = -1/D^2*dP_T ;
+dH_T = dU_T+1/D*dP_T ;
+dH_D = 1/D*dP_D - T/D^2*dP_T ;
 
 % step 2: rewrite the second-order derivatives (wrt d and T)
-d2S_T2 = 1/T*d2U_T2 - 1/T^2*dU_T
-d2S_D2 = -1/D^2*d2P_TD + 2/D^3*dP_T
-d2S_TD = -1/D^2*d2P_T2
+d2S_T2 = 1/T*d2U_T2 - 1/T^2*dU_T ;
+d2S_D2 = -1/D^2*d2P_TD + 2/D^3*dP_T ;
+d2S_TD = -1/D^2*d2P_T2 ;
+d2H_T2 = d2U_T2 + 1/D*d2P_T2 ;
+d2H_D2 = -T/D^2*d2P_TD + 1/D*d2P_D2 - 1/D^2*dP_D + 2*T/D^3*dP_T ;
+d2H_TD = 1/D*d2P_TD - 1/D^2*dP_T - T/D^2*d2P_T2 ;
 
 
 %% engineering derivs
@@ -39,7 +42,7 @@
 dP_D_S = expand(dP_D_S)
 
 
-%% dynamic derivs (first)
+%% dynamic derivs (first order)
 dT_P_H = 1/(dP_T - dP_D*dH_T/dH_D)
 dT_P_H = simplify(dT_P_H)
 dT_P_H = expand(dT_P_H)
@@ -57,21 +60,33 @@
 dD_H_P = expand(dD_H_P)
 
 
+%% dynamic derivs (second order)
+a1 = (d2P_T2*dH_D - d2P_TD*dH_T) / (dP_T*dH_D - dP_D*dH_T) ;
+b1 = (d2P_TD*dH_D - d2P_D2*dH_T) / (dP_T*dH_D - dP_D*dH_T) ;
+c1 = (d2H_T2*dH_D - d2H_TD*dH_T) / (dP_T*dH_D - dP_D*dH_T) ;
+d1 = (d2H_TD*dH_D - d2H_D2*dH_T) / (dP_T*dH_D - dP_D*dH_T) ;
+d2T_P2_H = (-dH_D^2*a1 + dH_T*dH_D*b1 + dP_D*dH_D*c1 - dP_D*dH_T*d1) / (dP_D*dH_T - dP_T*dH_D)^2
+d2T_P2_H = simplify(d2T_P2_H)
+d2T_P2_H = simplifyFraction(d2T_P2_H)
+%d2T_P2_H = expand(d2T_P2_H)
+
+
 %% PVT second order derivatives
 d2T_D2_P = -(d2P_D2*dP_T - dP_D*d2P_TD)/dP_T^2 ...
-           +(d2P_TD*dP_T - dP_D*d2P_T2)*dP_D/dP_T^3
-d2T_D2_P = simplify(d2T_D2_P)
-d2T_D2_P = expand(d2T_D2_P)
+           +(d2P_TD*dP_T - dP_D*d2P_T2)*dP_D/dP_T^3 ;
+d2T_D2_P = simplify(d2T_D2_P);
+d2T_D2_P = expand(d2T_D2_P);
 
 d2D_T2_P = -(d2P_T2*dP_D - dP_T*d2P_TD)/dP_D^2 ...
-           +(d2P_TD*dP_D - dP_T*d2P_D2)*dP_T/dP_D^3
-d2D_T2_P = simplify(d2D_T2_P)
-d2D_T2_P = expand(d2D_T2_P)
+           +(d2P_TD*dP_D - dP_T*d2P_D2)*dP_T/dP_D^3 ;
+d2D_T2_P = simplify(d2D_T2_P);
+d2D_T2_P = expand(d2D_T2_P);
 
 
 %% derivative for fundamental derivative of gas dynamics
 d2P_D2_S = d2P_D2 -(dP_T*d2S_D2   + 2*d2P_TD*dS_D)/dS_T ...
                   +(d2P_T2*dS_D^2 + 2*dP_T*dS_D*d2S_TD)/dS_T^2 ...
-                  -(dP_T  *dS_D^2*d2S_T2)/dS_T^3
-d2P_D2_S = simplify(d2P_D2_S)
-d2P_D2_S = expand(d2P_D2_S)
+                  -(dP_T  *dS_D^2*d2S_T2)/dS_T^3 ;
+d2P_D2_S = simplify(d2P_D2_S) ;
+d2P_D2_S = simplifyFraction(d2P_D2_S) ;
+d2P_D2_S = expand(d2P_D2_S) ;