versuche = 10000;
treffer = 0;

xz = pi*rand(versuche,1);
yz = rand(versuche,1);

for i = 1:versuche
  if yz(i) < sin(xz(i)) 
    treffer = treffer + 1;
  end
end 

func = @(x)sin(x);
korrekt = quad(func,0,pi);
schaetzung = pi*treffer/versuche;
fprintf('Schaetzung: %f \n',schaetzung);
fprintf('Korrekter Wert: %f \n',korrekt);

figure
hold on
scatter(xz,yz);
x = 0:0.01:pi;
plot(x,sin(x),'Color','red','LineWidth',2);
legend(['\fontsize{18} Schaetzung: ' num2str(schaetzung)]);
hold off