Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 31 vistas

    Tpnº 13 Método de Eliminacion de Gauss

    Cargado por

    Facundo Figueroa Paz Sosa
    El documento presenta un trabajo práctico sobre el método de eliminación de Gauss para resolver sistemas de ecuaciones lineales. Incluye el código fuente en FORTRAN que implementa el algoritmo de eliminación de Gauss para resolver un sistema de 3 ecuaciones con 3 incógnitas. Finalmente, muestra la salida del programa con la solución del sistema.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd

    Tpnº 13 Método de Eliminacion de Gauss

    Cargado por

    Facundo Figueroa Paz Sosa
    31 vistas4 páginas
    El documento presenta un trabajo práctico sobre el método de eliminación de Gauss para resolver sistemas de ecuaciones lineales. Incluye el código fuente en FORTRAN que implementa el algoritmo de eliminación de Gauss para resolver un sistema de 3 ecuaciones con 3 incógnitas. Finalmente, muestra la salida del programa con la solución del sistema.

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento presenta un trabajo práctico sobre el método de eliminación de Gauss para resolver sistemas de ecuaciones lineales. Incluye el código fuente en FORTRAN que implementa el algoritmo de eliminación de Gauss para resolver un sistema de 3 ecuaciones con 3 incógnitas. Finalmente, muestra la salida del programa con la solución del sistema.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    31 vistas4 páginas

    Tpnº 13 Método de Eliminacion de Gauss

    Cargado por

    Facundo Figueroa Paz Sosa
    El documento presenta un trabajo práctico sobre el método de eliminación de Gauss para resolver sistemas de ecuaciones lineales. Incluye el código fuente en FORTRAN que implementa el algoritmo de eliminación de Gauss para resolver un sistema de 3 ecuaciones con 3 incógnitas. Finalmente, muestra la salida del programa con la solución del sistema.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    Está en la página 1de 4

    Universidad Católica de Salta Cátedra: Análisis Numérico

    Facultad de Ingeniería e Informática


    Alumno: Facundo, Figueroa Paz Sosa DNI: 39537093
    Ingeniería Civil

    Trabajo Práctico Nº 13
    Tema: Sistemas de ecuaciones lineales “Método de
    ELIMINACION DE GAUSS”

    Actividades a realizar
    a) Realizar el programa en FORTRAN , dado el algoritmo de “ELIMINACION DE
    GAUSS” que resuelva:

    0 -1 2 x1 0
    -2 2 -1 x2 = 0
    -2 4 3 x3 1

    Valores iniciales

    x1=x2=x3=0

    Solución

    x1=1.8750 x10-1 ; x2=2.5000 x10-1 ; x3=1.2500 x10-1


    Universidad Católica de Salta Cátedra: Análisis Numérico
    Facultad de Ingeniería e Informática
    Alumno: Facundo, Figueroa Paz Sosa DNI: 39537093
    Ingeniería Civil

    CODIGO FUENTE

    PROGRAM ELIMINACIONDEGAUSS
    /* CSL/F6 - 1FOR ELIMINACION DE GAUSS
    DIMENSION A(10,11)
    PRINT*, ""
    PRINT*, 'CSL/F6-1 ELIMINACION DE GAUSS'
    DATA N/3/ /* N ES EL ORDEN DE LA MATRIZ
    DATA (A(1,J),J=1,4)/0,-1,2,0/ /*INICIALIZA LOS ELEMENTOS DE LA
    MATRIZ
    DATA (A(2,J),J=1,4)/-2,2,-1,0/ /*INICIALIZA LOS ELEMENTOS DE LA
    MATRIZ
    DATA (A(3,J),J=1,4)/-2,4,3,1/ /*INICIALIZA LOS ELEMENTOS DE LA
    MATRIZ
    PRINT*,""
    PRINT*, "MATRIZ AUMENTADA"
    PRINT*, ""
    DO I=1, N
    PRINT 61,(A(I,J), J=1,4)
    61 FORMAT (1X,1P6E12.4)
    END DO
    PRINT*,""
    CALL GAUSS(N,A)
    65 PRINT*, ""
    68 PRINT*, "SOLUCION"
    69 PRINT*, "---------------------------------"
    PRINT*, " I X(I)"
    70 PRINT*, "---------------------------------"
    DO I=1,N
    72 FORMAT (5X, I5, 1PE16.6)
    PRINT 72, I, A(I, N+1)
    END DO
    75 PRINT*, "-------------------------------"
    80 PRINT*, ""
    STOP
    END
    /******************
    SUBROUTINE GAUSS(N,A) /*ELIMINACION DE GAUSS
    INTEGER PV /* INDICE DE PIVOTEO
    DIMENSION A(10,11)
    EPS=1.0 /* SE CALCULA EL EPSILON DE LA MAQUINA
    10 IF (1.0+EPS.GT.1.0) THEN
    EPS=EPS/2.0
    GOTO 10
    END IF
    Universidad Católica de Salta Cátedra: Análisis Numérico
    Facultad de Ingeniería e Informática
    Alumno: Facundo, Figueroa Paz Sosa DNI: 39537093
    Ingeniería Civil

    EPS=EPS*2
    PRINT*, " EPSILON DE LA MAQUINA=", EPS
    EPS2=EPS*2
    1005 DET=1 /* INICIALIZACION DEL DETERMINANTE
    DO 1010 I=1, N-1
    PV=I
    DO J=I+1,N
    IF(ABS(A(PV,I)) .LT. ABS (A(J,I))) PV=J
    END DO
    IF (PV.EQ.I) GOTO 1050
    DO JC=1,N+1
    TM=A(I,JC)
    A(I,JC)=A(PV,JC)
    A(PV,JC)=TM
    END DO
    1045 DET=-1*DET /*CADA VEZ QUE SE REALICE UN
    PIVOTEO, CAMBIA EL SIGNO DE DET
    1050 IF(A(I,I).EQ.0) GOTO 1200 /* UNA MATRIZ SINGULAR SE
    A(I,I)=0
    DO JR=I+1, N /* ELIMINACION POR DEBAJO DE LA
    DIAGONAL
    IF (A(JR,I).NE.0) THEN
    R=A(JR,I)/A(I,I)
    DO KC=I+1,N+1
    TEMP=A(JR,KC)
    A(JR,KC)=A(JR,KC)-R*A(I,KC)
    IF (ABS(A(JR,KC)).LT.EPS2*TEMP) A(JR,KC)=0.0
    /* SI EL RESULTADO DE LA RESTA ES MENOR QUE
    /* EL DOBLE DEL EPSILON DE LA MAQUINA POR EL VALOR
    /* ORIGINAL, SE CAMBIA SU VALOR A CERO
    END DO
    END IF
    1060 END DO
    1010 CONTINUE
    DO I=1,N
    DET=DET*A(I,I)
    END DO
    PRINT*, ""
    PRINT*,"DETERMINANTE=", DET
    PRINT*, ""
    IF(A(N,N).EQ.0) GOTO 1200
    A(N,N+1)=A(N,N+1)/A(N,N)
    DO NV=N-1,1,-1
    VA=A(NV,N+1)
    DO K=NV+1,N
    VA=VA-A(NV,K)*A(K,N+1)
    END DO
    A(NV,N+1)=VA/A(NV,NV)
    Universidad Católica de Salta Cátedra: Análisis Numérico
    Facultad de Ingeniería e Informática
    Alumno: Facundo, Figueroa Paz Sosa DNI: 39537093
    Ingeniería Civil

    END DO
    RETURN
    1200 PRINT*, "LA MATRIZ ES SINGULAR"
    STOP
    END

    SALIDA

    GRAFICO

    NO SE GRAFICA PORQUE ES MATRIZ

    También podría gustarte