function data=trimalpha(alphamin,alphamax,Bhp,N_alpha,parameter)
% trim as a function of alpha for a given Bhp
% v and hdot will vary as a function of alpha
deltax=[1;1];
x0=[250;0];
u=zeros(2,1);
u(2)=Bhp;
inc=(alphamax-alphamin)/(N_alpha-1);
i=0;
for alpha=alphamin:inc:alphamax
	i=i+1;
	u(1)=alpha;
	[x,cost]=solve(x0,u,deltax,parameter);
	temp=[x;u;cost];
	data(i,:)=temp';
	x0=x;
end
%data=[v(ft/sec),hdot(ft/sec),alpha(rad),Bhp(hp),cost]
