Пример #1
0
int main()
{
    int n;
    int niz[10][10];

    printf("Unesite n: ");
    scanf("%d", &n);

    matrica (niz, n);
    return 0;
}
Пример #2
0
int main(){
	int i, j, k, a, b,step=1;
	float matr[n][n+n], matrob[n][n], x[n],xp[n], F[n], per[n], delit, koof, tmp,E, tmpE[2];
    printf("Vvedite Epsilon: ");
        scanf("%f", &E);
	x[0] = 2;	x[1] = 1;
	printf("\nX0 = %7.1f\nX1 = %7.1f\n", x[0], x[1]);	
//***********************************

do{
	matrica (&matr[0][0], x[0], x[1]);
// обратная матрица 
// приписываем единичную матрицу	
	for(i=0;i<n;i++){
		for(j=n; j<n+n; j++)
			if(i == j-n) matr[i][j] = 1.0;
				else matr[i][j] = 0.0;
	}
/*	for(i=0;i<n;i++){
		for(j=0;j<n+n;j++)
			printf("%7.2f ",matr[i][j]);
		printf("\n");
	}
*/	
// считаем обратную матрицу
	for(i=0; i<n; i++){
		delit = matr[i][i];
		for(j=0; j<n; j++){
			if(i==j) continue;
			koof = matr[j][i]/delit;
			for(k=i; k<n+n; k++){
				matr[j][k] -= matr[i][k]*koof;
				matr[i][k] /= delit;
			}
		}
	}
//***************************
	printf("\nObratnaya:\n");
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			matrob[i][j] = matr[i][j+n];
//			printf("%7.2f ", matrob[i][j]);
		}
//		printf("\n");
	}

	funct(&F[0], &F[1], x[0], x[1]);

// перемножаем матрицы
	xp[0] = x[0]; xp[1] = x[1];
    for(i=0; i<n; i++){
		tmp = 0.0;
		for(j=0; j<n; j++){
 			tmp += matrob[i][j]*F[j];
		}
			x[i] -= tmp;	
    }
    tmpE[0] = fabs(x[0]-xp[0]);
    tmpE[1] = fabs(x[1]-xp[1]);
	printf("\nX0 = %7.12f\nX1 = %7.12f", x[0], x[1]);
    printf("\nstep=%d",step++);	
    tmpE[0] = (tmpE[0] > tmpE[1] ? tmpE[0] : tmpE[1]);
    printf("\nPogrewhost: %f", tmpE[0]);
	printf("\n*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/\n");
}while(tmpE[0] > E);
	getch();
	return 0;
}