gaussiana N=2

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function S=gaussiana(a,b,n)
// n numero de intervalos
//GAUSSIANA PARA DOIS NÓS
h=(b-a)/n
x=linspace(a,b,n+1)
w1=1; t1=-sqrt(3)/3;
w2=w1; t2=-t1;
S=0
for i=1:n
 alpha=(x(i+1)-x(i))/2
 bet =(x(i+1)+x(i))/2 
 x1=alpha*t1+bet;
 x2=alpha*t2+bet;
 
 
 A =(w1*f(x1)+w2*f(x2))* h/2
 S=S+A
end
endfunction
function y=f(x)
 y=sin(sin(sin(x)))
endfunction