function [Ipm]=ptomedioc(a,b,M,f)
%PTOMEDIOC Metodo composito del punto medio 
%   [IPM]=PTOMEDIOC(A,B,M,F) calcola con il
%   metodo del punto medio composito l'integrale
%   approssimato IPM di una funzione F in [A,B]
%   usando M sotto-intervalli.
%
%   Si vedano anche TRAPC, SIMPADPT.

%   Saleri Fausto 23-01-01
 
H=(b-a)/M; 
x=[a+H/2:H:b-H/2]; 
fpm=eval(f); 
Ipm=H*sum(fpm);

return
