Il codice richiede di inserire soltanto il numero dei valori in cui si desidera dividere l'intervallo e la funzione. Maggiori saranno gli intervalli del dominio e più preciso sarà il calcolo dell' integrale definito.E' comunque dimostrato che il metodo Montecarlo ha una convergenza molto lenta e richiede, in generale, un numero elevato di intervalli per ottenere una stima precisa. %calcolo di integrali di funzioni in X=[0,1] function[]=metodo_montecarlo() disp('inserire numero di elementi del dominio') t=input('t='); X=linspace(0,1,t); disp('definire la funzione f(X),(x maiuscola)') f=input('f(X)='); y=f'; s=length(X); E=sum(y(:))/s; disp('integrale definito della funzione in X=[0,1],E='); disp(E);