Mudanças entre as edições de "Mmcc:aula 18"
De WikiLICC
m |
m |
||
| Linha 6: | Linha 6: | ||
[[Imagem:mmcc18_4.jpg]] | [[Imagem:mmcc18_4.jpg]] | ||
| + | |||
| + | 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) | ||
| + | |||
[[Imagem:mmcc18_5.jpg]] | [[Imagem:mmcc18_5.jpg]] | ||
Edição atual tal como às 13h39min de 6 de maio de 2010
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)
