Oefening 1

In de bovenstaande afbeelding zie je een mier in zijn natuurlijke leefomgeving: een rij hokjes. De mier zal gaan wandelen, op zoek naar voedsel of andere eerste levensbehoeften. We modelleren het wandelgedrag van de mier als volgt: iedere tijdseenheid (stap) loopt de mier met een kans van 50% 1 vakje naar links en zo niet, dan loopt hij 1 vakje naar rechts. Het doel voor de mier is om de suiker te bereiken in een van de uiterste vakjes.

Beantwoord de volgende vraag met een simulatie (experiment 10000 keer uitvoeren): hoeveel tijdseenheden (stapjes) duurt het gemiddeld tot de mier in het meest linkse vakje of het meest rechtse vakje aankomt?

Oefening 2

Het model in oefening 1 beschrijft een random walk in één dimensie. In deze oefening doe je een random walk in 2D.

Een man stapt om middernacht stomdronken uit een café midden in de stad. Hij wil naar huis, maar zijn richtingsgevoel en coördinatie laten hem enigszins in de steek: iedere stap die hij zet is in een willekeurige richting. Hij loopt steeds exact 1 stap in een willekeurige richting, noord, zuid, oost of west. Hoe groot is de kans dat de man in hooguit 10000 stappen terug komt bij het café?

Oefening 3

Wanneer we de lengte van de zijden $a$, $b$ en $c$ van een (mogelijke) driehoek willekeurig kiezen uit de uniforme verdeling $U(0, 1)$, wat is de kans dat er een daadwerkelijk een driehoek van te maken is?

Oefening 4

Bij de specialisatie Engineering wordt hard gewerkt aan het broodnodige integreren. Maar sommige functies zijn 'wat lastig' te integreren.

De integraal $\int_{0}^{1}\sqrt{1 - x^4} dx$ is een voorbeeld van zo'n functie:



In [18]:

    
%pylab
%matplotlib inline

x = linspace(0, 1, 400)
y = sqrt(1 - x**4)

plt.fill_between(x, y)
plt.axis('equal')
plt.ylim([0, 1.1])
plt.show()









    



Using matplotlib backend: Qt4Agg
Populating the interactive namespace from numpy and matplotlib

Met behulp van Monte Carlo simulatie kun je de integraal wel vrij eenvoudig benaderen. Je kunt een rechthoek kiezen waar je integraal binnen valt en dan voor willekeurige punten in de rechthoek bepalen of ze onder of boven de grafiek van de functie vallen.

Benader de gegeven integraal met een Monte Carlo simulatie.



In [ ]: