function [x,w] = ptomed_composta(a,b,N)

% Formula dei trapezi composta
% -----------------------------
% ---- input ----
% a,b : estremi di integrazione
% N : numero di subintervalli
%
% ---- output ----
% x : nodi di integrazione (vettore colonna)
% w : pesi di integrazione (vettore colonna)

% passo di integrazione
h = (b-a)/N;
% nodi di integrazione
x_ex = a:h:b; 
x = (x_ex(1:end-1)+x_ex(2:end))/2;
x = x';
% pesi di integrazione
w = h*ones(N,1);

end