Mmcc:aula 18
De WikiLICC
Usando o código (em Scilab)
a=1; b=3; // Intervalo [a,b] N=5; // N=numero de pontos Deltax=(b-a)/(N-1) // Tamanho de um intervalo x=linspace(a, b, N); // Pontos no intervalo y=1 ./ x; // y=f(x)=1/x plot(x,y,'r.-') pesos=ones(N,1); // [0 1 ... 1 1] pesos(1)=0; Area=y*pesos*Deltax // Area dos retangulos x=[x(1) x]; // Truque para plotar corretamente. y=[y y(N)]; plot2d2(x,y)