ideal gas limit, delta->0,

related to #11
thorade · Jul 30, 2013 · 85ddb57 · 85ddb57
1 parent 42cbfa7
commit 85ddb57
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/Interfaces/PartialHelmholtzMedium/EoS/Helmholtz_energy_derivatives.m b/Interfaces/PartialHelmholtzMedium/EoS/Helmholtz_energy_derivatives.m
@@ -46,7 +46,7 @@
 Phi = exp(-ni9*(delta-1)^2-ni10*(tau-1)^2)
 nonAnalyticalTerms = ni1*delta*Distance^ni6*Phi;
 % fr = +residualPolyTerms +BWRTerms +GaussTerms +nonAnalyticalTerms;
-fr = +nonAnalyticalTerms  
+fr = +GaussTerms  
 fr = simplify(fr)
 
 %% First derivatives
@@ -83,3 +83,6 @@
 frtdd = simplify(frtdd)
 frddd = diff(frdd,delta);
 frddd = simplify(frddd)
+
+%% ideal gas limit
+limit(frd,delta,0,'right')