写一个C++程序,在一个半径为 9.5 米的 90 度扇形中放置 16 个互不重叠的半径为 1 的圆，输出他的坐标,满足1<=x,y<8.5

#include <stdio.h>
#include <math.h>

#define PI 3.14159265358979323846 // 定义圆周率

int main()
{
double r1 = 1.0, r2 = 9.5; // 圆的半径和扇形的半径
double angle = 90.0; // 扇形的角度
int n = 16; // 圆的个数
double x, y; // 圆的坐标
int i;

for (i = 0; i < n; i++) {
double a = i * (angle / (n - 1)); // 计算圆心的角度
x = r2 * cos(PI * a / 180.0) + r1 * cos(PI * (a * (n - 1) / angle) / 180.0); // 计算圆心的横坐标
y = r2 * sin(PI * a / 180.0) + r1 * sin(PI * (a * (n - 1) / angle) / 180.0); // 计算圆心的纵坐标
printf("圆 %d 的坐标：%lf, %lf\n", i + 1, x, y); // 输出圆的坐标
}

return 0;
}