Ошибка LNK2019 - ссылка на неразрешенный внешний символ

typedef double(*funcX)(double);

#include <iostream>
#define _USE_MATH_DEFINES
#include <math.h>
 
#define eps 0.00001
#define step 0.001
 
using std::cout;
using std::cin;
using std::endl;
 
typedef double(*funcX)(double);
 
double integr_left_rect(double, double, funcX, double);
double func1_1(double);
double func1_2(double);
double func1_3(double);
double func2_1(double);
double func2_2(double);
double func3_1(double);
double func3_2(double);
double func4_1_1(double);
double func4_2_1(double);
double func4_3_1(double);
double func4_1_2(double);
double func4_2_2(double);
double func5_1(double);
double func5_2(double);
double func5_3(double);
 
int main()
{
    cout << "Добро пожаловать\n";
    cout << "ВВедите координаты точки x, y:\n";
    int x, y;
    cin >> x >> y;
 
    int i;
    if (x > -2 && x < 0 && y > 0 && y > (x + 1)* (x + 1) + y * y > 1 && x * x + (y - 1)*(y - 1) > 1)
    {
        cout << "Точка принадлежит M1\n";
        i = 1;
    }
    if (x * x + y * y < 1 && (x + 1) * (x + 1) + y * y < 1 && x * x + (y - 1)*(y - 1) < 1)
    {
        cout << "Точка принадлежит M2\n";
        i = 2;
    }
    if (x * x + y * y > 1 && x < 1 && (x - 1) * (x - 1) + y * y < 1)
    {
        cout << "Точка принадлежит M3\n";
        i = 3;
    }
    if (x * x + y * y > 1 && x * x + (y + 1)*(y + 1) < 1 && y < 0)
    {
        cout << "Точка принадлежит M4\n";
        i = 4;
    }
    if (x < 2 && x > 0 && x * x + (y + 1)*(y + 1) > 1 && y < -1 && y > -2)
    {
        cout << "Точка принадлежит M5\n";
        i = 5;
    }
    else
    {
        cout << "Не принадлежить областям.\n ";
    }
 
    cout.precision(5);
    switch (i)
    {
        double M1, M2, M3, M4, M5, sum;
    case 1:
        cout << "Геометрическая площадь" << (M1 = 1 + 2 * (1 - M_PI / 4));
        cout << "Интегральная площадь" << (M1 = integr_left_rect(-2, -1, func1_3, sum) - integr_left_rect(-1, 0, func1_2, sum) + integr_left_rect(-1, 0, func1_1, sum)) << endl;
        break;
    case 2:
        cout << "Геометрическая площадь" << (M2 = sqrt(3) / 4);
        cout << "Интегральная площадь" << (M2 = integr_left_rect(-sqrt(3)/2, 0, func2_1, sum) - integr_left_rect(-sqrt(3) / 2, 0, func2_2, sum))<< endl;
        break;
    case 3:
        cout << "Геометрическая площадь" << (M3 = M_PI / 4 - sqrt(3) / 4);
        cout << "Интегральная площадь" << (M3 = integr_left_rect(0.5 , 1, func3_1, sum) - integr_left_rect(0.5, 1, func3_2, sum)) << endl;
        break;
    case 4:
        cout << "Геометрическая площадь" << (M4 = M_PI / 4 + M_PI / 4 - sqrt(3) / 4);
        cout << "Интегральная площадь" << (M5 = (integr_left_rect(-1, -sqrt(3)/2, func4_1_1, sum) + integr_left_rect( -sqrt(3) / 2, 0, func4_2_1, sum) - integr_left_rect(-1, 0, func4_3_1, sum)) + (integr_left_rect(0, 2, func4_2_2, sum) + integr_left_rect( 0, 1, func4_2_1, sum))) << endl;
        break;
    case 5:
        cout << "Геометрическая площадь" << (M5 = 1 + (1 - M_PI / 4));
        cout << "Интегральная площадь" << integr_left_rect(0, 1, func5_1, sum) + integr_left_rect(1, 2, func5_2, sum) - integr_left_rect(0, 2, func5_3, sum) << endl;
        break;
 
    }
    system("pause");
}
 
double integral_left_rectangle(double a, double b, funcX func, double sum)
    {
        sum = 0;
        for (int i = a; i < b; i = i + step)
        {
            sum += func(i);
        }
        return sum * step;
    }
 
    double func1_1(double x)
    { 
        return sqrt(1 - (x - 1) * (x - 1)); 
    } 
 
    double func1_2(double x)
    {
        return sqrt(1 - x * x);
    }
 
    double func1_3(double x)
    {
        return 2;
    }
 
    double func2_1(double x)
    {
        return sqrt(1 - x * x);
    }
 
    double func2_2(double x)
    {
        return  -sqrt(1 - x * x) + 1;
    }
 
    double func3_1(double x)
    {
        return -sqrt(1 - x * x);
    }
 
    double func3_2(double x)
    {
        return -sqrt(1 - (x - 1) * (x - 1));
    }
 
    double func4_1_1(double x)
    {
        return sqrt(1 - x * x) - 1;
    }
 
    double func4_2_1(double x)
    {
        return -sqrt(1 - x * x);
    }
 
    double func4_3_1(double x)
    {
        return -1;
    }
 
    double func4_1_2(double x)
    {
        return - sqrt(1 - x * x) - 1;
    }
 
    double func4_2_2(double x)
    {
        return -1;
    }
 
    double func5_1(double x)
    {
        return -sqrt(1 - x * x) - 1;
    }
 
    double func5_2(double x)
    {
        return -1;
    }
 
    double func5_3(double x)
    {
        return -2;
    }