Найти точки пересечения двух линий заданных функциями

#include <iostream>
#include <cmath>
 
using namespace std;
/*y = ax^4 + bx^3 + cx^2 + dx + f и y = bx^3 + mx^2 + dx + p
a x^4  + (c-m)x^2  + f-p =  0
t=x2
b=c-m
c=f-p
at^2+bt+c=0
d=-b-4*a*c
t1=(b^2+sqr(d)) /2a
t2=(b^2-sqr(d)) /2a
 
x1=sqrt(t1)
x2=-sqrt(t1)
x3=sqrt(t2)
x4=-sqrt(t2)y = ax^4 + bx^3 + cx^2 + dx + f и y = bx^3 + mx^2 + dx + p
a x^4  + (c-m)x^2  + f-p =  0
t=x2
b=c-m
c=f-p
at^2+bt+c=0
d=-b-4*a*c
t1=(b^2+sqrt(d)) /2a
t2=(b^2-sqrt(d)) /2a
 
x1=sqrt(t1)
x2=-sqrt(t1)
x3=sqrt(t2)
x4=-sqrt(t2) */
void intersection(double a, double b, double c, double m, double d, double f, double p);
 
int main()
{   intersection(-1,-2,1,1,1,-1,-1.1);
    return 0;
}
void intersection(double a, double b, double c, double m, double d, double f, double p){
    double t1, t2, x1, x2, x3, x4;
    int i=0;
    b=c-m;
    c=f-p;
    d=-b-4*a*c;
    if(d<0){
        cout << "не пересекаются" << endl;
        return;
    }
    t1=(b*b+sqrt(d))/2*a;
    t2=(b*b-sqrt(d))/2*a;
    if(t1>=0){
        x1=sqrt(t1);
        x2=-sqrt(t1);
        i++;
        cout << "пересечение " << i << " x= "<< x1 << endl;
        i++;
        cout << "пересечение " << i << " x= "<< x2 << endl;
    }
    if(t2>=0){
        x3=sqrt(t2);
        x4=-sqrt(t2);
        i++;
        cout << "пересечение " << i << " x= "<< x3 << endl;
        i++;
        cout << "пересечение " << i << " x= "<< x4 << endl;
    }
    return;
}

d=b*b-4.0*a*c;

t1=(b*b+sqrt(d))/(2.*a);