=1</span>
0140 <span class="keyword">if</span> BFI&lt;0 || BFI&gt;1   
0141     BFI = NaN; 
0142     error_flag = 1;
0143     error_str = [<span class="string">'Warning: Estimated BFI outside allowed range (0 to 1).'</span>, error_str];
0144 <span class="keyword">end</span>
0145 
0146 <span class="comment">% optional plotting</span>
0147 <span class="keyword">if</span> plot_results
0148     fig = figure(<span class="string">'Position'</span>,[100 100 700 250]); hold on
0149     plot(t,Q)
0150     plot(t,Q_b