از سری کد ها و مجموعه ی کاربردی زبان متلب در علم ژئوماتیک در ادامه کد ها را ببینید:

clc
clear 
format long
a=input('enter the a parameter Vaule=');
f=input('enter the Scale factor ='); %f=1/298.25722221
b=a*(1-f);
e=2*f-f^2;
X=input('Enter X=');
Y=input('Enter Y=');
Z=input('Enter Z=');
P=sqrt(X^2+Y^2);
lan=2*atend((-X+sqrt(X^2+Y^2))/Y);
N0=a;
h0=((X^2+Y^2+Z^2)^.5)-(a*b)^.5;
phi0=atand(Z/P*(1-(e^2*N0)/(N0+h0)))^(-1);
N=1;h=1;phi=1;
while(abs(h-h0)>.00001)&&(abs(phi-phi0)>0000001)
    
    N=a/(1-(e^2)*(sin(phi0))^2)^.5;
    h=P/cos(phi0)-N0;
    phi=atand(z/p*(1-e^2*N0/(N0+h0))^-1);
    N0=N;h0=h;phi0=phi;
end
    phi=abs([fix((fix(phi)-phi)*60)+round(((ix((ix(phi)-phi)*60))-((fix(phi)-phi)*60))*60)]);
    phi=abs([fix((fix(lan)-lan)*60)+round(((ix((ix(lan)-phi)*lan))-((fix(lan)-lan)*60))*60)]);
   disp(phi);
   disp(lan);