Ошибка "нарушение доступа для чтения"

#include <iostream>
using namespace std;
class polyline {
    //Вершина ломаной
    struct polypoint {
        //Порядковый номер вершины (по порядку добавления)
        int n;
        //Координаты вершины ломаной
        double x;
        double y;
        //Тангенс касательной справа от данной точки
        double k;
        //Указатели на соседние вершины ломаной
        polypoint* prev;
        polypoint* next;
        polypoint(int inn, double inx, double iny, double ink, polypoint* inprev, polypoint* innext) {
            n = inn;
            x = inx;
            y = iny;
            k = ink;
            prev = inprev;
            next = innext;
        }
    };
    polyline() = delete;
private:
    int N; //Число точек ломаной
    polypoint* begin;
 
    //Функция поиска вершины с наименьшим значением координаты y
    polypoint* FindMinPoint() {
        polypoint* min = begin;
        for (polypoint* current = begin; current != nullptr; current = current->next)
            if (min->y > current->y)
                min = current;
        return min;
    }
    //Функция поиска вершины с порядковым номером n1
    polypoint* FindNPoint(int n1) {
        polypoint* current = begin;
        while ((current != nullptr) && (current->n != n1))
            current = current->next;
        return current;
    }
public:
    //Ломаная создается по двум точкам с углом сегментов ломаных справа от этих точек
    polyline(double x1, double y1, double k1, double x2, double y2, double k2) {
        N = 2;
        begin = new polypoint(1, x1, y1, k1, nullptr, nullptr);
        polypoint* end = new polypoint(2, x2, y2, k2, begin, nullptr);
 
        begin->next = end;
    }
    ~polyline() {
        while (begin != nullptr) {
            polypoint* temp = begin;
            begin = begin->next;
            delete temp;
        }
    }
    int GetN() { return N; } //Возвращает номер последней добавленной вершины (количество вершин ломаной)
    //Расстояние по x между точками с номерами n1 и n2
    double XDistance(int n1, int n2) {
        return abs(FindNPointX(n1) - FindNPointX(n2));
    }
    //Порядковый номер вершины с наименьшим значением координаты y
    int FindMinN() {
        return FindMinPoint()->n;
    }
    //Функции возвращают координаты и угол касательной вершины с номером n1
    double FindNPointX(int n1) {
        return FindNPoint(n1)->x;
    }
    double FindNPointY(int n1) {
        return FindNPoint(n1)->y;
    }
    double FindNPointK(int n1) {
        return FindNPoint(n1)->k;
    }
    //Возвращает порядковый номер вершины справа от данной
    int FindPrevN(int n1) {
        return FindNPoint(n1)->prev->n;
    }
    //Замена вершины с номером n1 на новую
    void ReplaceNPoint(int n1, double inx, double iny, double ink) {
        polypoint* npoint = FindNPoint(n1);
        npoint->x = inx;
        npoint->y = iny;
        npoint->k = ink;
    }
    //Добавление новой вершины справа от данной
    void AddPointAfterN(int n1, double inx, double iny, double ink) {
        polypoint* npoint = FindNPoint(n1);
        polypoint* newpoint = new polypoint(++N, inx, iny, ink, npoint, npoint->next);
        if (npoint->next != nullptr)
            npoint->next->prev = newpoint;
        npoint->next = newpoint;
    }
};
 
inline double F(double x) {
    return x * x * x * x + 2 * x * x + 4 * x + 1;
}
inline double Fd1(double x) {
    return 4 * x * x * x + 4 * x + 4 + 4;
}
//Координата x точки пересечения двух касательных 
inline double FindIntersection(double x1, double y1, double k1, double x2, double y2, double k2) {
    return (x1 * k1 - y1 - x2 * k2 + y2) / (k1 - k2);
}
//Значение касательной с тангенсом наклона k, проходящей через точку (x0,y0), в точке x 
inline double TangentY(double x0, double y0, double k, double x) {
    return y0 + k * (x - x0);
}
 
int main()
{
    double a = -1, b = 0; // Заданный отрезок
    double TargetPrecision = 1e-6; // Заданная точность
 
    //Точка пересечения касательных, выходящих из точек функции на границах заданного отрезка
    double x0 = FindIntersection(a, F(a), Fd1(a), b, F(b), Fd1(b));
    double y0 = TangentY(b, F(b), Fd1(b), x0);
    //Ломаная строится по точкам (a,f(a)), (x0,y0)
    polyline P = polyline(a, F(a), Fd1(a), x0, y0, Fd1(b));
    while (abs(F(P.FindNPointX(P.GetN())) - P.FindNPointY(P.GetN())) > TargetPrecision) {// |F(xk)-P(xk)|
        //Находится минимум ломаной и его координаты
        int minN = P.FindMinN();
        double x = P.FindNPointX(minN);
        double y = P.FindNPointY(minN);
        //Вычисляется значение функции над данной точкой
        double Fx = F(x);
        //Вычисляется значение касательной в этой точке функции
        double Fd1x = Fd1(x);
        //Углы ломаной справа и слева от данной точки
        double Kright = P.FindNPointK(minN);
        double Kleft = P.FindNPointK(P.FindPrevN(minN));
        //Точки пересечения касательной с сегментами ломаной справа и слева от точки
        double Xk1 = FindIntersection(x, y, Kleft, x, Fx, Fd1x);
        double Xk2 = FindIntersection(x, y, Kright, x, Fx, Fd1x);
 
        //Левая точка пересечения добавляется вместо текущей минимальной точки (x,y)
        P.ReplaceNPoint(minN, Xk1, TangentY(x, Fx, Fd1x, Xk1), Fd1x);
        //Правая точка пересечения добавляется к ломаной справа от текущей
        P.AddPointAfterN(minN, Xk2, TangentY(x, Fx, Fd1x, Xk2), Kright);
    }
    //while (P.XDistance(P.GetN(), P.GetN() - 1) > TargetPrecision) {// |P(xk) - P(xk-1)|   }
 
    int k = P.GetN();
    double xmin = P.FindNPointX(k);
    double fmin = F(xmin);
 
    cout.precision(7);
    cout << fixed << "k= " << k << ", x= " << xmin << ", f(x)= " << fmin;
 
    return 0;
}