想做国债收益率曲线，用三次样条函数构造，需要同时进行三个函数方程的回归，用到广义最小二乘法，请问如何进行三个方程的回归呢？

请好心人帮忙，多谢！