Witam, mam zadanie do wykonania a mianowicie napisanie prostego programu z pętlą for i instrukcją if. Wstęp programu już mam, prosiłbym o pomoc w dokończeniu i komentarzach przy każdej z linijek.
Program ma obliczać liczbę Pi metodą Monte Carlo. Coś oblicza, jednak zaraz po obliczeniu wyłącza się i nie jestem pewnien czy robi to tak jak należy.
function p =pi % oblicza Pi metodą MC
nkolo = 0;
rand('state',sum(100clock)); % zmienia "seed" przy każdym uruchomieniu
n = input('Podaj n = ') % prosi o podanie liczby losowanych punktów
for k = 1: n
x = rand;
y = rand;
r = sqrt(x2 + y2);
if (r <= 1) % jeżeli wylosowany punkt jest w kole o promieniu r =1
nkolo = nkolo + 1;
end%if
end%for
p = 4nkolo/n;
d = abs(pi - p); % błąd bezwzględny
dw = 100*d/pi; % błąd względny podany w procentach
disp("Pi = num2str(p)")
disp("delta = num2str(d)")
disp(["delta % = num2str(dw)")