% Quick trim the BoilerExpress
disp('Start here')
Vt=28% ft/sec
Hp=1000
constant=BoilEX;
[xIC,uIC,CL,CD,CM,alphadeg]=QuickTrim(Vt,Hp,constant)
flag=1;
t=0;
[sys,x0,str,ts] = aircraft(t,xIC,uIC,flag,constant,xIC)