075 <span class="keyword">end</span>;
0076 
0077 <span class="comment">%-------------------------------------------------</span>
0078 <span class="comment">% project observation point (lon,lat) =&gt; (x,y)</span>
0079 <span class="comment">%-------------------------------------------------</span>
0080 [xobs,yobs] = my_project(lon,lat,<span class="string">'forward'</span>); 
0081 
0082 
0083 <span class="comment">%-------------------------------------------------</span>
0084 <span class="comment">% compute e