Написать программу на си для аппроксимации спалайн-функциями. У меня есть код, нужно дописать функцию spline

#include <stdio.h>
#include <math.h>
#include <conio.h>
 
void spline(int n, float x[n], float y[n], float r[n][3], float s[n-1][4]);
 
int main(int argc, char **argv)
{
    int n;      //  количество узлов аппроксимации
    int i, j;       // рабочие переменные для ввода исходных данных
    char c0;
    printf("\nUzlov approksimacii  n= ");               
    scanf("%d", &n);        
    float x[n], y[n];           // исходные узлы аппроксимации
    float r[n][3], s[n-1][4];   // массивы коэффициентов
// аппроксимирующего кубического сплайна
                    // по возрастанию степеней 
    printf("\nVvedite vector argumenta (%d znacheniy):\n", n);
    for (i=0; i<n; i++)
        scanf("%f", &x[i]);
    printf("\n");
    printf("\nVvedite vector funkcii (%d znacheniy):\n", n);
    for (i=0; i<n; i++)
        scanf("%f", &y[i]);
    printf("\n");
    printf("\nApproksimacija kubicheskimi splajnami\n");
        
    spline(n, x, y, r, s);
    
    printf("\nApproksimir. kubich. splajny\n");
    printf("\nUchastok\t\tkoefficienty approks. splajna\n");
    for (i=0; i<n-1; i++)
    {
        printf("\n   %d\t\t", i+1);
        for (j=0; j<4; j++)
            printf("%10.5f", s[i][j]);
    }
    printf("\n");
    c0=getch();
    return 0;
}