t = [0 2.5 5 10 15 20]';
P = [7.5 10.5 12.5 15.8 17.9 19.4]';
P0 = 7.5;