function Qw4=wpcn4(Pw4,well4)
x=well4(:,1);
Y=well4(:,2:4);
Qw4= interp1(x,Y,Pw4)';