Построить график

//-------------------------------------------------------------------------------------------------
#include <windows.h>
#include <iostream>
#include <cmath>
 
//-------------------------------------------------------------------------------------------------
double func(double x)
{
    // корень не определён для отрицательных чисел
    // подкоренное выражение становится отрицательным при X от минус бесконечности до единицы
    if (x < 1)
        return 3.0;
 
    return sqrt(pow(x, 3) - 1) + sin(pow(x, 2));
}
 
//-------------------------------------------------------------------------------------------------
int main()
{
    HWND hwnd = GetConsoleWindow();
    HDC hdc = GetDC(hwnd);
 
    // задаём диапазон изменений X
    double src_x_min = -2.5;
    double src_x_max = 2.5;
    double src_cx = src_x_max - src_x_min;
    // задаём диапазон изменений Y
    double src_y_min = -1;
    double src_y_max = 4;
    double src_cy = src_y_max - src_y_min;
 
    // динамически считаем диапазон Y (что бы включить рассчёт нужно заменить false на true)
    if (false)
    {
        // ищем минимум и максимум Y на заданном диапазоне
        double min_y = func(src_x_min);
        double max_y = func(src_x_min);
        for (double x = src_x_min; x < src_x_max; x += 0.005)
        {
            double y = func(x);
 
            if (min_y > y)
                min_y = y;
 
            if (max_y < y)
                max_y = y;
        }
        std::cout << "min: " << min_y << std::endl;
        std::cout << "max: " << max_y << std::endl;
 
        // посчитанный диапазон Y
        src_y_min = min_y;
        src_y_max = max_y;
        src_cy = max_y - min_y;
    }
 
    // размеры окна
    RECT rc;
    GetClientRect(hwnd, &rc);
    int width = rc.right - rc.left;
    int height = rc.bottom - rc.top;
 
    // посчитаем параметры области вывода (сделаем 10% отступы по краям)
    double dst_x0 = width / 10.;
    double dst_y0 = height / 10.;
    double dst_cx = width * 8 / 10.;
    double dst_cy = height * 8 / 10.;
 
    // нарисуем оси
    {
        HPEN pen = CreatePen(PS_SOLID, 1, RGB(128, 128, 128));
        HPEN prevPen = (HPEN)SelectObject(hdc, pen);
 
        double dst_x = (0 - src_x_min) * dst_cx / src_cx + dst_x0;
        MoveToEx(hdc, dst_x, 0, 0);
        LineTo(hdc, dst_x, height);
 
        double dst_y = (0 - src_y_min) * dst_cy / src_cy + dst_y0;
        dst_y = height - dst_y;
        MoveToEx(hdc, 0, dst_y, 0);
        LineTo(hdc, width, dst_y);
 
        SelectObject(hdc, prevPen);
        DeleteObject(pen);
    }
 
    // нарисуем график
    for (double src_x = src_x_min; src_x < src_x_max; src_x += 0.005)
    {
        double src_y = func(src_x);
 
        double dst_x = (src_x - src_x_min) * dst_cx / src_cx + dst_x0;
        double dst_y = (src_y - src_y_min) * dst_cy / src_cy + dst_y0;
        dst_y = height - dst_y;
 
        SetPixel(hdc, dst_x, dst_y, RGB(255, 255, 255));
    }
 
    ReleaseDC(hwnd, hdc);
    std::cin.ignore();
    return 0;
}
 
//-------------------------------------------------------------------------------------------------