function[x,y,z]=XYZ(Lat,Long,Alt)
dLat=Lat-Lat(1);
dLong=Long-Long(1);
x=dLat.*40075160.*cosd(Lat)./360;
y=dLong.*40008000./360;
z=Alt.*0.3048;
