Dato un vettore v=(vx,vy,vz) e una funzione scalare a(x,y,z),si vuole dimostrare la seguente relazione:
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);  

Commenti
Posta un commento