double wektorP[4];
	double wek1p1[4];
	double wek1p2[4];
	double sumawek1[4];
	double wek2p1[4];
	double wek2p2[4];
	double sumawek2[4];
	double wek3p1[4];
	double wek3p2[4];
	double sumawek3[4];
	double wek4p1[4];
	double wek4p2[4];
	double sumawek4[4];

	for (i = 0; i < 4; i++)
	{
		wek1p1[i] = alfa * N1pow1[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		wek1p2[i] = alfa * N2pow1[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		sumawek1[i] = wek1p1[i] * wek1p2[i];
	}

	for (i = 0; i < 4; i++)
	{
		wek2p1[i] = alfa * N1pow2[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		wek2p2[i] = alfa * N2pow2[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		sumawek2[i] = wek2p1[i] * wek2p2[i];
	}

	for (i = 0; i < 4; i++)
	{
		wek3p1[i] = alfa * N1pow3[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		wek3p2[i] = alfa * N2pow3[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		sumawek3[i] = wek3p1[i] * wek3p2[i];
	}

	for (i = 0; i < 4; i++)
	{
		wek4p1[i] = alfa * N1pow4[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		wek4p2[i] = alfa * N1pow4[i] * talfa * tabJ[i];
	}

	for (i = 0; i < 4; i++)
	{
		sumawek4[i] = wek4p1[i] * wek4p2[i];
	}




	for (i = 0; i < 4; i++)
		{
			wektorP[i] = tabPow[0] * sumawek1[i] + tabPow[1] * sumawek2[i] + tabPow[2] * sumawek3[i] + tabPow[3] * sumawek4[i];
		}