fuel_left/gauss_samps.m

   1 function s = gauss_samps(n,split,split2)\r
   2 \r
   3 s = zeros(5,n);\r
   4 for i = 1:n\r
   5     if rand < split\r
   6         %faster ROS\r
   7         s(1:4,i) = 1000*randn(4,1);\r
   8     else\r
   9         %slower ros\r
  10         s(1:4,i) = 5000*randn(4,1);\r
  11     end\r
  12 end\r
  13 %set fuel time constants\r
  14 for i = 1:n\r
  15     if rand < split2\r
  16         %slower burning fuels\r
  17         s(5,i) =  8 + 2000*rand;\r
  18     else\r
  19         %fast burning fuels\r
  20         s(5,i) =  8 + 10*rand;\r
  21     end\r
  22 end\r
  23 \r
  24 \r
  25 % figure,histogram(s(1,:))\r
  26 % title('t00')\r
  27 \r
  28 \r
  29 end % function\r