Buenas tardes Marianela, 

Dejo resolución adjunta en drive: 

https://drive.google.com/file/d/1kmWYvmuY5R8yiXLtUKY1MH4VJub8IX7h/view?usp=sharing

Al utilizar la ecuación integral, se debe trabajar con la densidad promedio entre la condición de entrada y la condición de salida, según deducción a partir de la ecuación de los gases ideales (adjunta en el archivo). En caso de resolver con la integral resuelta previamente, no se requiere el valor de densidad.

Como bien mencionas, si bien se pre-define con un valor aleatorio la variable que se desea resolver, siempre es conveniente hacerlo con un valor próximo a lo esperable, sino Mathcad puede arrojar error o, en el caso del uso de otras funciones de resolución como "root" encuentra la raíz mas próxima al valor pre-definido. 

Gracias por tu observación :)