function start
global bond usclose gold ct 
global dm 
global usbp

K=menu('which data set', '30 year us bond','US treasury bill', 'gold price', 'cotton price', 'US$ Dmark exchange rate','US$British pd exchange','simulated (1800 daily values)');
par=[3.2, -.4, .01 2.5];
st=str2mat('z=bond(:,2);' ,'z=usclose;' ,'z=gold;' ,'z=ct;','z=dm;','z=usbp;', '[z,x,parstart]=milstei2(1800,.004);');
st2=str2mat('parstart=[.0052 2.5910];', 'parstart=[5.1616 .1252];','parstart=[.004 1.62];', 'parstart=[.3029 .9391];','parstart=[.0778 .4979];','parstart=[.0988 1.2158];','parstart=parstart;');
s=st(K,:);
eval(s);
eval(st2(K,:));
pause(4);
bonddemo(z,2,.004,parstart)