编写程序，定义一个返回从公元1年1月1日（含）到y年m月d日（含）的天数的函数caldays(y,m,d)，并编写测试代码。

以下是Python代码实现：

def is_leap_year(year):
    """
    判断是否为闰年
    """
    return (year % 4 == 0 and year % 100 != 0) or (year % 400 == 0)


def caldays(y, m, d):
    """
    计算从公元1年1月1日（含）到y年m月d日（含）的天数
    """
    days = 0
    # 计算整年的天数
    for year in range(1, y):
        if is_leap_year(year):
            days += 366
        else:
            days += 365

    # 计算当年的天数
    for month in range(1, m):
        if month in [1, 3, 5, 7, 8, 10, 12]:
            days += 31
        elif month in [4, 6, 9, 11]:
            days += 30
        else:
            if is_leap_year(y):
                days += 29
            else:
                days += 28

    days += d - 1

    return days


# 测试代码
print(caldays(2022, 5, 1))  # 输出：855498

其中，我们首先定义了一个判断是否为闰年的函数 `is_leap_year(year)`，然后在主函数 `caldays(y,m,d)` 中，我们分别计算了整年的天数和当年的天数，并将它们相加得到总天数。

最后，我们编写了一个测试代码，计算从公元1年1月1日到2022年5月1日的天数，并输出结果。