function [geo]=utm2geo(E,N)

k0=0.9996; %scalefactor
a=6378137.00000; %semi-major axis
b=6356752.31414; %semi-minor axis
f= 0.003352810681225; %flattening
e= sqrt(f*(2-f));%eccentricity
FN=0;%false northing
FE=500000.00;%false easting

landa0=-3*pi/180;%??

M0=0;

M1=M0+(N-FN)/k0;
mu1=M1/(a*(1-(e*e)/4-(3*e*e*e*e)/64-(5*e*e*e*e*e*e)/256));
e1=(1-sqrt(1-e*e))/(1+sqrt(1-e*e));
phi1=mu1+(3*e1/2-27*e1*e1*e1/32)*sin(2*mu1)+(21*e1*e1/16-55*e1*e1*e1*e1/32)*sin(4*mu1)+(151*e1*e1*e1/96)*sin(6*mu1)+(1097*e1*e1*e1*e1/512)*sin(8*mu1);
ro1=(a*(1-e*e))/((1-e*e*sin(phi1)*sin(phi1))^(3/2));
u1=a/sqrt(1-e*e*sin(phi1)*sin(phi1));

D=(E-FE)/(u1*k0);
ep2=(e*e)/(1-e*e);
C1=ep2*cos(phi1)*cos(phi1);
T1=tan(phi1)*tan(phi1);

landa=landa0+(D-(1+2*T1+C1)*D*D*D/6+(5-2*C1+28*T1-3*C1*C1+8*ep2+24*T1*T1)*D*D*D*D*D/120)/cos(phi1);
phi=phi1-(u1*tan(phi1)/ro1)*(D*D/2-(5+3*T1+10*C1-4*C1*C1-9*ep2)*D*D*D*D/24+(61+90*T1+298*C1+45*T1*T1-252*ep2-3*C1*C1)*D*D*D*D*D*D/720);
geo(1)=landa;
geo(2)=phi;