Υπολογισμός του π με προσομοίωση Monte Carlo στο Octave
Το κλασικό πρόβλημα υπολογισμού του π. Η λύση που δίνεται εδώ είναι ενδεικτική και σίγουρα όχι η πλέον βέλτιστη. Είναι χρήσιμη ωστόσο από εκπαιδευτική άποψη.
Κάτι που πρέπει να κάνετε είναι να πειραματιστείτε με την παράμετρο N και να δείτε την επίδραση των αλλαγών στο αποτέλεσμα. Μπορείτε να βγάλετε κάποιο συμπέρασμα;
clear; N = 5000; # Number of points x = unifrnd(0, 1, N, 1); y = unifrnd(0, 1, N, 1); r = sqrt(x .^ 2.0 + y .^2.0); N1 = sum((r<1.0)); mcpi = 4*N1/N; mpci
Συνδεθείτε για περισσότερες δυνατότητες αλληλεπίδρασης,
σχολιασμοί, εξωτερικοί σύνδεσμοι, βοήθεια, ψηφοφορίες, αρχεία, κτλ.
σχολιασμοί, εξωτερικοί σύνδεσμοι, βοήθεια, ψηφοφορίες, αρχεία, κτλ.
Εκπαιδευτικό υλικό από τον
Αθανάσιο Σταυρακούδη
σας παρέχετε κάτω από την άδεια
Creative Commons Attribution-NonCommercial-ShareAlike 4.0 License.
Σας παρακαλώ να ενημερωθείτε για κάποιους επιπλέον περιορισμούς
http://stavrakoudis.econ.uoi.gr/stavrakoudis/?iid=401.