close all
clear all

c0 = 7897;
cinf = 245;

V = 1000;

T = 4:.1:20;

H = 2;  %m

Pa = cinf*H^2*T.^3;
Pb = c0*V*H./T;

plot(T,Pa,T,Pb);
