
clear
valg=2; % 1: f=exp(x)     2: f=exp(x).*sqrt(1-x.^2);
%vi måler tiden!
t=cputime;
n=100;

% Gitteret
theta=linspace(0,pi,n);
x=cos(theta);

%funksjonsverdier

if valg==1
    f=exp(x);
end

if valg==2
    f=exp(x).*sqrt(1-x.^2);
end

%vekter
w=zeros(n,1);
w(1)=1/(n-1)^2;w(n)=w(1);
for i =2:n-1
    w(i)=2*(1-sum(2./(4*(1:(n-2)/2).^2-1).*cos(2*pi*(1:(n-2)/2)*(i-1)/(n-1))))/(n-1);
end

%kvadraturen
int=f*w;

%sammenlikning med eksakt verdi
if valg==1 
    disp(abs(int-exp(1)+exp(-1)))
end
if valg==2   
    disp(abs(int-1.7754996892121809469))
end

%cputime-t

