LU разложение матрицы

[Артегор]

хотел написать программу для LU разложения матрицы.
взял формулы с вики , не правильно находит элемент (3,3) в матрице U.
(пробовал только матрицу 3х3 ,иожет при больших размерах ещё больше неточностей будет )
что не так ?
#include <iostream>


int main(int argc, char* argv[])
{
int i,j,k, n=3;
float sum;

float ** A=new float *[n];
float ** L=new float *[n];
float ** U=new float *[n];

for (int i=1; i<=n; i++) {
A[i]=new float [n];
L[i]=new float [n];
U[i]=new float [n]; }

std::cout<<"A matr \n";

for (i=1;i<=n;i++){
for (j=1;j<=n;j++){
A[i][j]=rand()%9;
L[i][j]=0;
U[i][j]=0;
L[i][i]=1;
std::cout<<A[i][j]<<" ";
}std::cout<<"\n";}

/////////// L и U марицы


for (j=1;j<=n;j++) {
U[1][j]=A[1][j];}
for (j=2;j<=n;j++) {
L[j][1]=A[j][1]/A[1][1]; }


for (i=2;i<=n;i++){
for (j=i;j<=n;j++){
sum=0;
for (k=1;k<i;k++){sum=sum+L[i][k]*U[k][j];}
U[i][j]=A[i][j]-sum;

}}



for (i=2;i<=n;i++){
for (j=i+1;j<=n;j++){
sum=0;
for (k=1;k<i;k++){sum=sum+L[j][k]*U[k][i];}
L[j][i]=(A[j][i]-sum)/U[i][i];
L[i][i]=1;
}}



std::cout<<"L matr \n";
for (i=1;i<=n;i++){
for (j=1;j<=n;j++){
std::cout<<L[i][j]<<" ";
}std::cout<<"\n";}

std::cout<<"U matr \n";
for (i=1;i<=n;i++){
for (j=1;j<=n;j++){
std::cout<<U[i][j]<<" ";
}std::cout<<"\n";}


system("pause");
return 0;

}

[Ezhik Kihze]

Начните с того, что прочитайте как нумеруются массивы в с++...

[Артегор]

Сообщение от Ezhik Kihze

Начните с того, что прочитайте как нумеруются массивы в с++...

о боже,ну с нуля ,что с того то ?это всего лишь смещение
разницы в подсчете не будет

[whatever]

Сообщение от Артегор

о боже,ну с нуля ,что с того то ?это всего лишь смещение
разницы в подсчете не будет

Я очень поверхностно знаю C++, но тем не менее:
Вот тут ты задаешь размерность массива

A[i]=new float [n];

То есть если n=3 то, А[i] состоит из А[i][0], A[i][1] и A[i][2]
Далее

for (i=1;i<=n;i++){
 for (j=1;j<=n;j++){

пробегаются элементы A[i][1], A[i][2] и A[i][3], который не существует. Так что это не "всего лишь смещение".