function [ret] = fwd ( LAT, LONG, BEARING, DIST )
  LAT1r=LAT*pi/180;
  LONG1r=LONG*pi/180;
  BEARINGr=mod(BEARING,360)*pi/180;
  LAT2=asin(sin(LAT1r)*cos(DIST/3959)+cos(LAT1r)*sin(DIST/3959)*cos(BEARINGr))*180/pi;
  LAT2r=LAT2*pi/180; 
  LONG2= (LONG1r + atan2(sin(BEARINGr)*sin(DIST/3959)*cos(LAT1r), cos(DIST/3959)-sin(LAT1r)*sin(LAT2r)))*180/pi;
  x=[LAT2 LONG2];
  ret=x;
endfunction
