%  TEST PROGRAM FOR SEQUENCE OF COMMANDS
%  ***************************************
H=100;
T=18;
d=150;
g=32.2;
HTdg=[H T d g]
say='[dol0,hohb,L0,Hb]=dimwv(H,T,d,g)'
[dol0,hohb,L0,Hb]=dimwv(H,T,d,g)
say='[hol0,hod,lol0,psi]=wvcon(dol0,hohb)'
[hol0,hod,lol0,psi]=wvcon(dol0,hohb)
say='yc=ycol0(dol0,hohb,hol0)'
yc=ycol0(dol0,hohb,hol0)

say='coef=stfn(dol0,hohb)'
coef=stfn(dol0,hohb);
say='Computing theta and y vectors'

theta=[0:1:180]';
y=[0:yc/100:yc]';
say='eta=profile(theta,dol0,hohb)'
eta=profile(theta,dol0,hohb);
say='[eta,vx,vy,ax,ay]=wvprop(theta,y,dol0,hohb,hol0,lol0,coef)'
[eta,vx,vy,ax,ay]=wvprop(theta,y,dol0,hohb,hol0,lol0,coef);