用c语言写一个已知一条直线和圆的方程，求它们的公切圆方程的代码

以下是用C语言编写已知一条直线和圆的方程，求它们的公切圆方程的代码：

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

int main()
{
    double a, b, c, h, k, r, d, e, f;
    printf("请输入直线的方程系数a,b,c：");
    scanf("%lf%lf%lf", &a, &b, &c);
    printf("请输入圆的方程的圆心坐标h,k和半径r：");
    scanf("%lf%lf%lf", &h, &k, &r);
    d = a * a + b * b;
    e = -2 * a * h - 2 * b * k + 2 * c;
    f = h * h + k * k - r * r + c * c - 2 * c * k;
    if (e * e - 4 * d * f < 0)
    {
        printf("无实根\n");
    }
    else
    {
        double x1, x2, y1, y2;
        x1 = (-e + sqrt(e * e - 4 * d * f)) / (2 * d);
        y1 = (-a * x1 - c) / b;
        x2 = (-e - sqrt(e * e - 4 * d * f)) / (2 * d);
        y2 = (-a * x2 - c) / b;
        printf("公切圆方程的参数为：\n");
        printf("圆心坐标为：(%lf, %lf)\n", (x1 + x2) / 2, (y1 + y2) / 2);
        printf("半径为： %lf\n", sqrt(((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2))) / 2);
    }
    return 0;
}

该程序首先要求用户输入直线的方程系数a,b,c以及圆的方程的圆心坐标h,k和半径r。然后，程序计算出公切圆方程的参数，并输出圆心坐标和半径。