// Example 5_7
clc;funcprot(0);
// Given data
D=0.100;// m
T_1=200;// °C
p_1=0.140;// MPa
h=3.50;// W/(m^2.K)
T_infinitive=15.0;// °C
c_v=3.123;// kJ/kg.K
R=2.077;// kJ/kg.K
t=5.00;// seconds

// Calculation
V=(%pi/6)*D^3;// m^3
A=%pi*D^2;// m^2
m=((p_1*10^3)*V)/(R*(T_1+273.15));// kg
hAbymc_v=(h*A)/(m*c_v*1000);// s^-1
T_2=((T_1-T_infinitive)*exp((-(h*A)/(m*c_v*1000))*t))+T_infinitive;// °C
delU=m*c_v*(T_2-T_1);// kJ
printf('\n(a)The final temperature of the helium,T_2=%2.1f°C \n(b)The change in total internal energy of the helium,U_2-U_1=%0.3f kJ',T_2,delU);