C++随机生成圆内的一个坐标（double） 非均匀

极坐标原理，设圆心为x,y;半径为r。
新坐标（a,b）。
则
a= cos ⁡ α cosalpha cosα*r+x
b= sin ⁡ α sinalpha sinα*r+y

double PI=3.1415926;
//从[0，2*pi]中随机一个double值
double radian=((double)rand()/RAND_MAX)*2*PI; 
//从[0，r]中随机一个double值
double randR=((double)rand()/RAND_MAX)*r;

return {randR*cos(radian)+x,randR*sin(radian)+y};//随机坐标

C++11有个随机函数是均匀的，但使用语法过于麻烦…