function [XH,XDH,XDDH]=project1(TP,TS,XP,XDP,BETA)
T=0.;
X=XP;
XD=XDP;
H=.001;
while T<=(TS-.0001)
	XDD=.0034*32.2*XD*XD*exp(-X/22000.)/(2.*BETA)-32.2;
	XD=XD+H*XDD;
	X=X+H*XD;
	T=T+H;
end
XH=X;
XDH=XD;
XDDH=XDD;