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

//// phi1 : parallèle de référence de la projection.
//// lambda : longitude du méridien central.
//// exec('Projection_Bonne_inv.sci');

// Ayrens 	2.333333 44.983333
phi1=45;
lambda0=2.337229;

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


//////////
yB=1./tan(phi1)-y;
rhoB=x.^2+yB.^2;
rho=sqrt(rhoB);

phi=1./tan(phi1)+phi1-rho;
lambda=lambda0+rho./cos(phi).*atan(x./yB);

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

endfunction