Написать программу на си

Насколько я понял тему сначала нужно выяснить координаты точек P0, P1, P2 и P3 для двух кубических кривых и прикинуть принадлежат ли каждая из этих точек четырехугольнику, в который вписывается другая кривая.

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

// Функция для вычисления значения кубического полинома
double cubic(double a, double b, double c, double d, double x) {
return a * pow(x, 3) + b * pow(x, 2) + c * x + d;
}

// Производная кубического полинома
double cubic_derivative(double a, double b, double c, double x) {
return 3 * a * pow(x, 2) + 2 * b * x + c;
}

// Метод Ньютона-Рафсона для нахождения корня
double newton_raphson(double a, double b, double c, double d, double x0, double epsilon) {
double x = x0;
double f_x, f_prime_x;
do {
f_x = cubic(a, b, c, d, x);
f_prime_x = cubic_derivative(a, b, c, x);
x = x - f_x / f_prime_x;
} while (fabs(f_x) > epsilon);
return x;
}

// Функция для нахождения всех корней кубического уравнения
void find_roots(double a, double b, double c, double d, double epsilon) {
double roots[3];
int num_roots = 0;

// Начальные приближения для корней
double x0 = -10.0;
double x1 = 0.0;
double x2 = 10.0;

// Найдем корни
roots[0] = newton_raphson(a, b, c, d, x0, epsilon);
roots[1] = newton_raphson(a, b, c, d, x1, epsilon);
roots[2] = newton_raphson(a, b, c, d, x2, epsilon);

// Проверка на уникальность корней
for (int i = 0; i < 3; i++) {
if (i == 0 || fabs(roots[i] - roots[i - 1]) > epsilon) {
printf("x = %.6f\n", roots[i]);
num_roots++;
}
}

if (num_roots == 0) {
printf("Кривые не пересекаются.\n");
} else {
printf("Кривые пересекаются в следующих точках:\n");
for (int i = 0; i < num_roots; i++) {
double y = cubic(a + l, b, c + w, d + v, roots[i]);
printf("Точка пересечения: (x = %.6f, y = %.6f)\n", roots[i], y);
}
}
}

int main() {
double p, o, n, m;
double l, w, v;
double epsilon = 1e-6;

// Ввод коэффициентов кривых
printf("Введите коэффициенты первой кривой (p, o, n, m): ");
scanf("%lf %lf %lf %lf", &p, &o, &n, &m);

printf("Введите коэффициенты второй кривой (l, w, v): ");
scanf("%lf %lf %lf", &l, &w, &v);

// Коэффициенты кубического уравнения
double a = p - l;
double b = o;
double c = n - w;
double d = m - v;

// Поиск корней
find_roots(a, b, c, d, epsilon);

return 0;
}