CREATE DEFINER = `root`@`%` FUNCTION `NewProc`(`lon1` real,`lat1` real,`lon2` real,`lat2` real)
RETURNS float
BEGIN
#Routine body goes here...
DECLARE dist FLOAT;
DECLARE radLat1 FLOAT;
DECLARE radLat2 FLOAT;
DECLARE a FLOAT;
DECLARE b FLOAT;
DECLARE radius real;
SET radius = 6371000;
SET radLat1=lat1*PI()/180.0;
SET radLat2=lat2*PI()/180.0;
SET a=radLat1-radLat2;
SET b=lon1*PI()/180.0-lon2*PI()/180.0;
SET dist=2*ASIN(SQRT(POW(SIN(a/2),2)+COS(radLat1)*COS(radLat2)*POW(SIN(b/2),2)));
SET dist=dist*radius;
SET dist=ROUND(dist*10000)/10000;
RETURN dist;
END;