Function presents data wrong, Monte Carlo method
 

<p>I am trying to write a program that does integration with monte carlo method. One of its features is to place dots on the graph with different colours, blue or red depending on the if statement. The if statement is put in a "for" loop and i dont know why but it seems like the first option is ignored after first iteration. The whole thing looks like this :
<a href="https://i.stack.imgur.com/2rEwV.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/2rEwV.png" alt="enter image description here"></a></p>

<p>but it should look like this :</p>

<p><a href="https://i.stack.imgur.com/OvaPz.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/OvaPz.png" alt="enter image description here"></a></p>

<p>In addition i dont know why but it looks like the plot makes some additional empty space at the top</p>

<p>The whole code is not finished yet, its just a matter of a few lines but these dots are so annoying that I want to figure out whats wrong first. Heres the code.</p>

<pre><code>function p=montecarlo(f, a, b, n, t)
%f is a function provided by user
%a and b is a range
%n is the amount of random points
%t is a t=a:01:b vector to draw a plot
upper=max(f(t));
lower=min(f(t));
x=a+(b-a).*(rand(n,1)) %generates vector of random numbers from a to b
y=lower+(upper-lower).*(rand(n,1)) %generates vector of ranom numbers from min to max
hold on
for i=1:n
if y(i)>=f(i)
plot(x(i),y(i),'bo')
else
plot(x(i),y(i),'ro')
end
plot(t,f(t),'k')
end
end
</code></pre>

<p>Arguments provided to the function : f= x.^2+3*x+5, a= -4 , b= 2, n= 1000 . </p>



[url=https://stackoverflow.com/questions/59010160/function-presents-data-wrong-monte-carlo-method]More...[/url]