计算距已知坐标一定距离内的经纬度取值范围

import math
d = 1000
R=6370996.81
pi=3.1415926

lng1=96.486654
lat1=39.490295

dx= math.acos((math.cos(d/R)-math.sin(lat1*pi/180)*math.sin(lat1*pi/180))/(math.cos(lat1*pi/180)*math.cos(lat1*pi/180)))*180/pi
dy=180*d/R/pi

#lng2=lng1
lng2=lng1-dx
#lat2=lat1+dy
lat2=lat1

distance=R*math.acos(math.cos(lat1*pi/180 )*math.cos(lat2*pi/180)*math.cos(lng1*pi/180 -lng2*pi/180)+math.sin(lat1*pi/180 )*math.sin(lat2*pi/180))
print distance

print dx
print dy