function out=fVol(type,teta)
% Programme de calcul du volume et de la variation de volume par degré vilebrequin
% type: 'vol' et 'dvol'
% teta: angle ou vecteur ligne d'angle 
global Cu Vm lambd 
r=2*pi/360;

switch type
case 'vol' %Volum de la chambre
   out=Vm+Cu/2*(1+lambd-cos(teta*r)-(lambd^2-sin(teta*r).^2).^0.5);
case 'dvol' %variation de volume par degré vilebrequin 
   out=pi/180*Cu/2*sin(teta*r).*(1+cos(teta*r)./(lambd^2-sin(teta*r).^2).^0.5);
end