(val-.8)/1.1)^2, 1, 1000);
y = ( (val - 1.1*(t.^(1/2)))/.8 ).^2;
plot(t,y,<span class="string">'--'</span>, [.5 .5], [0 1], <span class="string">':'</span>, [0 1], [.5 .5], <span class="string">':'</span>);
axis([0 1 0 1]);
hold <span class="string">on</span>

<span class="keyword">for</span> k=1:m
   <span class="keyword">if</span> bounds(k,2) &lt; 1e-5,  <span class="comment">% preferred over (.5,.5)</span>
      dot = plot(data(k,1),data(k,2),<span class="string">'o'</span>);
      <span class="comment">