t=0:0.01:350;
p1=1.0+(0.00277*50-1.71*10^(-7)*50^3)*(1.0-1.15)^2+3.49*10^(-4)*50./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p2=1.0+(0.00277*25-1.71*10^(-7)*25^3)*(1.0-1.15)^2+3.49*10^(-4)*25./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p3=1.0+(0.00277*0.1-1.71*10^(-7)*0.1^3)*(1.0-1.15)^2+3.49*10^(-4)*0.1./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p4=0.88+(0.00277*50-1.71*10^(-7)*50^3)*(0.88-1.15)^2+3.49*10^(-4)*50./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p5=0.88+(0.00277*25-1.71*10^(-7)*25^3)*(0.88-1.15)^2+3.49*10^(-4)*25./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p6=0.88+(0.00277*0.1-1.71*10^(-7)*0.1^3)*(0.88-1.15)^2+3.49*10^(-4)*0.1./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p7=0.78+(0.00277*50-1.71*10^(-7)*50^3)*(0.78-1.15)^2+3.49*10^(-4)*50./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p8=0.78+(0.00277*25-1.71*10^(-7)*25^3)*(0.78-1.15)^2+3.49*10^(-4)*25./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
p9=0.78+(0.00277*0.1-1.71*10^(-7)*0.1^3)*(0.78-1.15)^2+3.49*10^(-4)*0.1./[0.927+3.18*10^(-4)*(t+17.78).^1.175];
plot(t,p1,'r',t,p2,'g',t,p3,'b',t,p4,'k',t,p5,'c',t,p6,'y',t,p7,'k--',t,p8,'k:',t,p9,'r-.');
