Autovalor e autovetor de uma matriz

[Thiago Heidrich]

Pessoal, bom dia.

Estou com um baita problema p resolver. Tenho um trabalho p entregar sobre autovalor e autovetor de uma matriz, consegui um código que acha esses dados q desejo. Ele não estava depurando, consegui depurar até metade do código porém tá dando erro.

 Estou usando o microsoft visual studio. O erro está em anexo.

Esse é o link com o código original:https://www.vivaolinux.com.br/script/Metodo-de-Power-para-calcular-o-autovelor-dominante-de-uma-matriz

Esse aqui é o código.

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include<malloc.h>
#define MAX(a,b) c = (absval(a) >= absval(b)) ? a : b // usado  no calculo da norma infinita.

FILE * fp; // arquivo texto com o autovetor e seu autovalor associado

double* cria_vetor(int tam) { /* rotina para alocar vetores dinamicamente */
    double* vetor;

    vetor = (double*)calloc(tam, sizeof(double));
    if (!vetor) {
        printf("Falta memoria para alocar o vetor de ponteiros");
        exit(1);
    }
    return vetor;
}

double absval(double val) // valor absoluto com precisao dupla
{
    if (val >= 0) return val;
    else return -val;
}

void dados(int dim, int* N, double* TOL, double** A, double* x) {  /* recebe os dados da
                                                                     matriz de entrada */
    for (int i = 0; i < dim; i++)
        for (int j = 0; j < dim; j++)
            A[j] = 0.0;

    /* entrada da matriz */
    for (int i = 0; i < dim; i++) {
        printf("\n Linha #%d: \n", i + 1);
        for (int j = 0; j < dim; j++) {
            printf("\n    Coluna #%d: ", j + 1);
            scanf_s("%lf", &A[j]);
        }
    }
    printf("\n Entre com o vetor X inicial:"); /* vetor inicial X(0) */
    for (int i = 0; i < dim; i++) {
        printf_s("\n x#%d: ", i + 1);
        scanf_s("%lf", &x);
    }
    printf("\n Entre com a Tolerancia: ");
    scanf_s("%lf", TOL);
    printf("\n Entre com o numero maximo de iteracoes: ");
    scanf_s("%d", N);
    system("clear");
    fprintf(fp, "Matriz de entrada\n");
    for (int i = 0; i < dim; i++) {
        for (int j = 0; j < dim; j++)
            fprintf(fp, "%lf ", A[j]);
        fprintf(fp, "\n");
    }
    fprintf(fp, "\n");
} // fim dados

double norma_inf(double* x, int dim) {
    int i;
    double c;                        // norma infinita do vetor
    i = 0;
    c = x[0];
    for (i = 1; i < dim; i++)
        c = MAX(c, x);
    if (c == 0) {
        printf("\n A matriz \"A\" possui um autovalor 0, selecione um novo vetor \"x\" e reinicie!");
        fprintf(fp, "\nA matriz \"A\" possui um autovalor 0, selecione um novo vetor \"x\" e reinicie!");
        exit(1);
    }
    return c;
} // fim norma_inf

void exibe(double* x, double lambda, int dim, int k) { // exibe os dados computados
    int i;

    if (k == 0) {
        printf("n\t lambda(n)\t\t\t X(n)' \n");
        printf("-----------------------------------------------------------------\n");
        fprintf(fp, "n\t lambda(n)\t\t\t X(n)' \n");
        fprintf(fp, "-----------------------------------------------------------------\n");
    }
    printf("%d\t ", k + 1);
    fprintf(fp, "%d\t ", k + 1);
    printf("%lf\t [ ", lambda);
    fprintf(fp, "%lf\t [ ", lambda);
    for (i = 0; i < dim; i++) {
        printf("%lf ", x);
        fprintf(fp, "%lf ", x);
    }
    printf("]'\n");
    fprintf(fp, "]'\n");
} // fim exibe

void power(int dim, int N, double TOL, double** A, double* x) { // calcula o autovalor e o autovetor
    double* y, lambda, c1, * aux, erro;
    int i, j, k, FLAG;
    // alocacao dinamica de memoria
    y = cria_vetor(dim);
    aux = cria_vetor(dim);
    i = 0;
    j = 0;
    k = 0;
    lambda = 0;
    FLAG = 0;
    c1 = norma_inf(x, dim);
    for (i = 0; i < dim; i++)
        x = x / c1;
    while ((k < N) && (FLAG == 0)) {
        for (i = 0; i < dim; i++) {
            y = 0;
            for (j = 0; j < dim; j++)
                y += A[j] * x[j];
        }
        c1 = norma_inf(y, dim);

        for (int l = 0; l < dim; l++)
            aux[l] = x[l] - y[l] / c1;
        erro = norma_inf(aux, dim);
        for (i = 0; i < dim; i++)
            x = y / c1;
        FLAG = 0;
        if (fabs(erro) < TOL) FLAG = 1;
        lambda = c1;
        exibe(x, lambda, dim, k);
        k++;
    }
    free(y);
} // fim power

int main() {
    double TOL;
    double** A, * x;
    int dim, N;

    system("clear");
    printf(" ----------------------------------------------\n");
    printf("  Calcula o autovalor dominante de uma matriz  \n");
    printf("  quadrada e o seu autovetor associado pelo    \n");
    printf("  metodo da potencia.                          \n");
    printf(" ----------------------------------------------\n");
    printf("\n Dimensao da matriz ? "); scanf_s("%d", &dim);
    // alocacao dinamica de memoria
    A = (double**)calloc(dim, sizeof(double)); // matriz de entrada
    for (int i = 0; i < dim; i++) A = (double*)calloc(dim, sizeof(double));
    x = cria_vetor(dim);  // vetor inicial
    dados(dim, &N, &TOL, A, x);

    power(dim, N, TOL, A, x);
    free(A);
    free(x);
    fclose(fp);
    
    
    return 0;
}