Added isentropic.m