uiopen('e:\Math 447\polarGPS.m',1) % Modify expression to add input arguments. % Example: % a = [1 2 3; 4 5 6]; % foo(a); polarGPS a a = 15600 11760 11610 15170 b b = 7540 2750 14630 610 c c = 10380 16190 7680 13320 rho rho = 26570 phi phi = 0.401355155089376 0.655220171546875 0.293232029093986 0.525120502950172 theta theta = 0.313380971199111 0.130904999001313 0.612814916506477 0.026536329364190 R R = 1.0e+004 * 1.778459164557905 1.556574765310038 1.872284967626456 1.669740997879611 t t = 0.059423012207262 0.052021745319892 0.062552704117942 0.055796564517297 TruePos = GPSCalc([0;0;6370;0],zeros(4,1)); err(:,1) = GPSCalc([0;0;6370;0],[1;0;1;0]); err(:,2) = GPSCalc([0;0;6370;0],[-1;0;1;0]); err(:,2) = GPSCalc([0;0;6370;0],[-1;-1;1;1]); err(:,2) = GPSCalc([0;0;6370;0],[-1;0;1;0]); err(:,3) = GPSCalc([0;0;6370;0],[-1;-1;1;1]); err(:,3) = GPSCalc([0;0;6370;0],[1;-1;1;-1]); err(:,4) = GPSCalc([0;0;6370;0],[0;-1;1;-1]); TruePos- err(:,1) ans = 0.012359814694069 0.009886534855839 0.002341249935853 0.000000041894346 TruePos- err(:,2) ans = -0.039907392616364 -0.018858631608364 -0.016082120543615 -0.000000145528319 TruePos- err(:,3) ans = 0.024719650043469 0.019773085658853 0.004682502171818 0.000000093788779 TruePos- err(:,4) ans = -0.001414040682569 0.005400450037504 -0.004529206234110 0.000000000077099 infNorm = norm(TruePos- err(:,1),inf) infNorm = 0.012359814694069 for k=1:4 infNorm(k) = norm(TruePos- err(:,k),inf) end infNorm = 0.012359814694069 infNorm = 0.012359814694069 0.039907392616364 infNorm = 0.012359814694069 0.039907392616364 0.024719650043469 infNorm = Columns 1 through 3 0.012359814694069 0.039907392616364 0.024719650043469 Column 4 0.005400450037504 max(infNorm) ans = 0.039907392616364 errMagFact = ans/.003 errMagFact = 13.302464205454593