function [x,y]=Projection_Bonne(lambda,phi)

// S.Ayrinhac (02/2012)
// lambda : longitude en dégrés décimaux
// phi : latitude en degrés décimaux

// phi1 : parallèle de référence de la projection.
phi1=45;
// lambda0 : longitude du méridien central.
lambda0=2.337229;

//// conversion des degrés décimaux en radians
phi=phi*%pi/180;
phi1=phi1*%pi/180;
lambda=lambda*%pi/180;
lambda0=lambda0*%pi/180;


beta0=phi;
beta1=phi1;

R=1./tan(beta1)+beta1-beta0;
epsilon=(lambda-lambda0).*cos(beta0)./R;

x=R.*sin(epsilon);
y=1./tan(beta1)-R.*cos(epsilon);

endfunction