actual final(?) update to polytrope.m