function OK = convertCoordinates_test_uk(varargin)
%CONVERTCOORDINATES_TEST_UK  test for UK points
%
%    OK = convertCoordinates_test_uk(<EPSG>)
%
%See also: convert_coordinates

   if nargin==1
   EPSG = varargin{1};
   else
   EPSG = load('EPSG');
   end

% test for UK system (fails currently)
% http://osgeo-org.1560.n6.nabble.com/DHDN-Gauss-Krueger-to-WGS84-td3847546.html

geo_code = 31466; % 4326 = WGS 84

lon = [11.446588];
lat = [49.169494];

x   = [4459750];
y   = [5448182];

[x1,y1] = convertCoordinates(lon ,lat , EPSG, 'CS1.code', geo_code, 'CS2.code', 4326);

OK = abs(x1-x) < 1 && abs(y1-y) < 1; % require 1 m accurate
