GPS Software Receivers - Fundamentals of Global Positioning System Receivers: A Software Approach - page 216

%p96.m SATPOSM.M Use ephemeris data to calculate satellite
position
% modified for data generated from OUR OWN COLLECTED DATA
receiver
% generat find PR from coarse PR
% JT 24 Sept 96
function[outp]=p(inpdat);
tc=inpdat(1);
toe=inpdat(2);
deln=inpdat(3);
asq=inpdat(4);
ra=inpdat(5);
i=inpdat(6);
w=inpdat(7);
delra=inpdat(8);
M=inpdat(9);
idot=inpdat(10);
e=inpdat(11);
crs=inpdat(12);
cuc=inpdat(13);
cus=inpdat(14);
cic=inpdat(15);
cis=inpdat(16);
crc=inpdat(17);
af0=inpdat(18);
af1=inpdat(19);
af2=inpdat(20);
toc=inpdat(21);
tgd=inpdat(22);
PRc=inpdat(23);
%**** define data
wer=7.2921151467e-5;
GM=3986005e8;
a=asq^2;
mu=GM;
n=(mu/(a^3))^.5+deln;
if tc-toe>302400;
tc=tc-604800;
elseif tc-toe<-302400;
tc=tc+604800;
end
M=M+n*(tc-toe);
Eold=M;