function [xave,hpave]=mpaveline(x,hp,MP,mptest)
i=find(MP<=mptest);
if isempty(i)
	disp('mptest is out of range')
end
imax1=i(1);
imax2=imax1-1;
if imax2<=0;
	imax2=1;
end
[itest1,xtest1,hptest1,MPtest1]=mpline(imax1,x,hp,MP);
[itest2,xtest2,hptest2,MPtest2]=mpline(imax2,x,hp,MP);

if imax2==imax1
	xave=xtest1;
	hpave=hptest1;
else
	f=(mptest-MP(imax2))/(MP(imax1)-MP(imax2));
	xave=f*xtest1+(1-f)*xtest2;
	hpave=f*hptest1+(1-f)*hptest2;
end
