0001 function [emf, W] = varyWind() 0002 W = 106; 0003 for i = 1:3000, 0004 theta = tacoma([0 500],[0 0 0.001 0], 0.04, 3, W ); 0005 emf = theta / 0.001; 0006 W = W - 0.005; 0007 if( emf < 100.0005 ) break; 0008 end 0009 end