metoda banachiewicza

Question

Jak zakodowac ten(E suma od s=1 do j-1, lis*usj) fragment algorytmu w c++? Napisalem taki kod ale nie dziala

#include <iostream>
using namespace std;
/*
4x1 +2x2 +2x3 -2x4 +2x5=-2
4x1 -5x2 +0.5x3 +2x4 +0.5x5=0.5
-5x1 +0.5x2 +2x3 -4x4 +0.5x5=-2
-2x1 +2x2 +0.5x3 -5x4 -4x5=0.5
0.5x1 +2x2 -4x3 -3x4 +2x5=0.5
*/
int main()
{
	const int n = 5;
    int i,j,s;
    float sum,a=0;
    float A[n][n]={{4,2,2,-2,2},{4,-5,0.5,2,0.5},{-5,0.5,2,-4,0.5},{-2,2,0.5,-5,-4},{0.5,2,-4,-3,2}};
    float L[n][n]={{0,0,0,0,0},{0,0,0,0,0},{0,0,0,0,0},{0,0,0,0,0},{0,0,0,0,0}};
	float U[n][n]={{1,0,0,0,0},{0,1,0,0,0},{0,0,1,0,0},{0,0,0,1,0},{0,0,0,0,1}};
	float B[n]={-2,0.5,-2,0.5,0.5};
	
	for(i=0;i<n;i++)
	{
		L[i][1]=A[i][1];
	}
	
	for(j=1;i<n;i++)
	{
		U[1][j]=A[1][j]/L[1][1];
	}
	
	for(j=1;i<n;i++)
	{
		L[i][2]=A[i][2]-L[i][1]*U[1][2];
	}
	
	for(j=2;i<n;j++)
	{
				for(i=1;i<j-1;i++)
				{	
				
					for(s=0;s<i-1;s++)
					{
						a=a+L[i][s]*U[s][j];
					}
					U[i][j]=(A[i][j]-a)/L[i][i];
				}a=0;
				
				for(i=j;i<n;i++)
				{	
				
					for(s=0;s<j-1;s++)
					{
						a=a+L[i][s]*U[s][j];
					}
					L[i][j]=A[i][j]-a;
				}a=0;
		
	}
   
	return 0;
}

 

Comments

Powiem szczerze, że nie chcę mi się czytać tego kodu. Ale spróbuję Ci doradzić jak to ruszyć. Przeprowadź kawałek tych obliczeń na kartce, a następnie sprawdź debuggerem, czy Twój program daje takie same wyniki. Jeśli nie wiesz co to debugger i nie chcesz się dowiedzieć, to możesz wszystko wypisywać na ekran, np.

a=a+L[i][s]*U[s][j];
std::cout << "a=a+L[i][s]*U[s][j];" << std::endl;
std::cout << "a = " << a << std::endl;
std::cout << "i = " << i << std::endl;
// ... i tak dalej ...
std::cout << "U[s][j] = " << U[s][j] << std::endl;

Trochę się napracujesz, ale błąd powinieneś znaleźć :)