Mmcc:aula 18

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)