Program stopped working, when tried to run

I wrote a program for Lagrange Interpolation which do compiles. Unfortunately when I run the program in Code Blocks it only shows warning "Program stopped working".

In Dev-C++ it returns value over 32 000 000.

Can you check the code and tell me why? There you have code:

#include <iostream>
#include<stdlib.h>
using namespace std;

struct wielo
{
    double a, b, c, d, e, f;
};

double *mnozenie(double A[], double B[], int m, int n)
{
    double *prod = new double[m + n - 1];

    for (int i = 0; i < m + n - 1; i++)
        prod[i] = 0;

    for (int i = 0; i < m; i++)
    {
        for (int j = 0; j < n; j++)
            prod[i + j] += A[i] * B[j];
    }

    return prod;
    delete prod;
}

int main()
{
    int wierz = 3; // max 6
    int sizeW = wierz - 1;
    double wynik = 0;
    double szukana = 5;
    double tabX[wierz];
    double tabY[wierz];
    double *tabTT;
    double wynikowy[wierz];

    tabX[0] = 1;
    tabX[1] = 2;
    tabX[2] = 4;

    tabY[0] = 3;
    tabY[1] = 4;
    tabY[2] = 6;

    wielo w[] =
    {
    { 0, 0, 0, 0, 0, 0 },
    { 0, 0, 0, 0, 0, 0 },
    { 0, 0, 0, 0, 0, 0 } };

    wynikowy[0] = 0.0;
    wynikowy[1] = 0.0;
    wynikowy[2] = 0.0;
    int k = 0;
    for (int i = 0; i < wierz; i++)
    {
        double dzialanie = 0;
        dzialanie = tabY[i];

        for (int j = 0; j < wierz; j++)
        {
            if (tabX[i] != tabX[j])
            {
                dzialanie *= (szukana - tabX[j]);
            }
        }

        for (int j = 0; j < wierz; j++)
        {
            if (tabX[i] != tabX[j])
            {
                dzialanie /= (tabX[i] - tabX[j]);
            }
            if (tabX[i] != tabX[j] && i < j)
            {
                double *tabT1;
                double *tabT2;
                tabT1[0] = 1;
                tabT1[1] = tabX[j] * (-1);
                tabT2[0] = 1;
                tabT2[1] = tabX[i] * (-1);
                tabTT = mnozenie(tabT1, tabT2, sizeW, sizeW);
                w[k].a = tabTT[0];
                w[k].b = tabTT[1];
                w[k].c = tabTT[2];
                k++;
            }
        }
        wynik += dzialanie;
        cout << dzialanie << endl;
    }
    double x, y, z;
    x = tabY[0] / ((tabX[0] - tabX[1]) * (tabX[0] - tabX[2]));
    y = tabY[1] / ((tabX[1] - tabX[0]) * (tabX[1] - tabX[2]));
    z = tabY[2] / ((tabX[2] - tabX[0]) * (tabX[2] - tabX[1]));
    for (int i = 0; i < wierz; i++)
    {
        if (i == 0)
        {
            wynikowy[0] = wynikowy[0] + (w[i].a * x);
            wynikowy[1] = wynikowy[1] + (w[i].b * x);
            wynikowy[2] = wynikowy[2] + (w[i].c * x);
        }
        if (i == 1)
        {
            wynikowy[0] = wynikowy[0] + (w[i].a * y);
            wynikowy[1] = wynikowy[1] + (w[i].b * y);
            wynikowy[2] = wynikowy[2] + (w[i].c * y);
        }
        if (i == 2)
        {
            wynikowy[0] = wynikowy[0] + (w[i].a * z);
            wynikowy[1] = wynikowy[1] + (w[i].b * z);
            wynikowy[2] = wynikowy[2] + (w[i].c * z);
        }

    }
    cout << "Wzor to: " << wynikowy[0] << "X^2 " << wynikowy[1] << "X " << wynikowy[2] << endl;
    cout << "Wynik dla x=" << szukana << " " << "y=" << wynik << endl;

    return 0;
}