역행렬 1000*1000 을구하려고하는데요 고쳐주세요~
새누
2023.11.23
#include stdio.h
#define N 3 /* 원의 수 */
int main(int argc, char* argv[])
{ /* 계수행렬 */
static double a[N][N+1];
double p,d;
int i,j,k;
for (i=0;iN;i++)
{
for(j=0;jN+1;j++)
{
printf(%lf,&a[i][j]);
}
}
for (k=0;kN;k++) {
p=a[k][k]; /* 피벗 계수 */
for (j=k;jN+1;j++) /* 피벗행을 p로 나눈다. */
a[k][j]=a[k][j]/p;
for (i=0;iN;i++) { /* 피벗열 소거 */
if (i!=k){
d=a[i][k];
for (j=k;jN+1;j++) /* ① */
a[i][j]=a[i][j]-d*a[k][j];
}
}
}
for (k=0;kN;k++)
printf(x%d=%f\n,k+1,a[k][N]);
return 0;
}
이소스로 1000*1000까지 구해야하는데 어디를손봐야하나욤