CalculoNumerico/LU.c at master · comth/CalculoNumerico · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
//LU

#include<stdio.h>

int main(){
	printf("Fatoracao LU")
    int k,i,n,j,aux,aux1;
    printf("Quantidade de variaveis: ");
    scanf("%d",&n);

    float a[n][n],m[n][n];

    for(aux=0;aux<n;aux++){
        for(aux1=0;aux1<n;aux1++){
            printf("Valor a[%d][%d]: ",aux+1,aux1+1);
            scanf("%f",&a[aux][aux1]);
        }
    }

    for(aux=0;aux<n;aux++){
        for(aux1=0;aux1<n;aux1++){
            if(aux==aux1)
                m[aux][aux1] = 1.0;
            else
                m[aux][aux1] = 0.0;
        }
    }

    for(k=0;k<(n-1);k++){
        for(i=k+1;i<n;i++){
            m[i][k] = a[i][k]/a[k][k];
            for(j=k+1;j<n;j++)
                a[i][j] = a[i][j]-(m[i][k]*a[k][j]);
            a[i][k] = 0;
        }
    }

    printf("\n\n");

    printf("Matriz U:\n");
    for(aux=0;aux<n;aux++){
        for(aux1=0;aux1<n;aux1++){
            printf("\t%.3f",a[aux][aux1]);
        }
        printf("\n");
    }

    printf("Matriz L:\n");
    for(aux=0;aux<n;aux++){
        for(aux1=0;aux1<n;aux1++){
            printf("\t%.3f",m[aux][aux1]);
        }
        printf("\n");
    }

    printf("\n\n");

    return 0;

}