Métodos numéricos

Gauss-Seidel. Diagrama de Flujo. Gauss-Jordan. Newton-Raphson. Programas de métodos

  • Enviado por: Rockdrigo
  • Idioma: castellano
  • País: México México
  • 5 páginas

publicidad
cursos destacados
Curso de Keynote: Diseña y contruye tus presentaciones fácil y rápido
Curso de Keynote: Diseña y contruye tus presentaciones fácil y rápido
Keynote es el software para realizar presentaciones por computadora más poderoso del mercado. Con este curso...
Ver más información

Excel 2010 para Dummies
Excel 2010 para Dummies
¡Este curso es para TI si eres un DUMMY en Excel! Este curso básico, como su nombre lo indica, está...
Ver más información


Fórmula general del método.

Codificación (Programa).

DIMENSION X(5), T(5)

READ (5,30)N

FORMAT (12)

READ (5,31) (X(I), I=1,N)

FORMAT (5F3.0)

WRITE (6,32) (X(I),I=1,N)

FORMAT (5(2X,F3.0))

READ (5,31)(X(I),I=1,N)

WRITE(6,32)(X(I), I=1,N)

C=0

DO 6 I=1,N

C=C+X(I)*T(I)

CONTINUE

WRITE (6,33)C

FORMAT(E15.8)

CALL EXIT

END

Aproximaciones sucesivas.

Fórmula General

Codificación(Programa)

READ (2,1)N

1 FORMAT(I2)

K=0

READ (2,4) A,B E1, E2

4 FORMAT (2F4.0,2F9.6)

7 X=(A+B)/2

WRITE(3,5)X

5 FORMAT (5X, 2HX=, E15.8)

K=K+1

IF (ABS (X*2-25.))-E1)6,6,8

8 IF ((ABS(B-A))-E2)6,6,10

10 IF (N-K)6,6,12

12 IF ((X*2-25.)*(A*2-25.)) 13,6, 15

15 A=X

GO TO 7

13 B=X

GO TO 7

6 CALL EXIT

END

Newton-Raphson

Diagrama de bloques (pseudo código).

Diagrama de Flujo

Codificación(Programa)

WRITE (3,101)

FORMAT (1H1,//,10X,36H DETERMINACIÓN DE UNA RAÍZ APROXIMADA,/,10X,31 1H POR EL METODO DE NEWTON RHAPSON,//)

READ (2, 102) X, E, N

FORMAT (F5, 0,F10, 6, 13)

K=1

FX=X**3-4.3*X*2+4.*X-17.2

DFX=3.*X**2-8. 6*X+4

XNUEVA=X X-(FX/DFX)

WRITE (3,103)XNUEVA

FORMAT (10X,7HX NUEVA=E15,8)

IF (ABS (XNUEVA-X)-E 8,89

IF(N-K)8,8,7

K=K+1

X=XNUEVA

GO TO 3

STOP

END

Operaciones de Matrices

SUMA

Codificación (Programa)

DIMENSION A(3,3), B(3,3), C(3,3)

READ (5,1)M,N

FORMAT(2I2)

READ (5,2) ((A(I,J),J=1, N), I=1, M

READ (5,2) ((B(I,J),J=1, N), I=1, M

FORMAT (3F.0)

WRITE (6,3)((A(I,J), J=1, N) , I=1, M)

WRITE (6,3)((B(I,J), J=1, N) , I=1, M)

FORMAT (3(2X,F5.0))

DO 6 I=1, M

DO 6 J=1, N

C(I,J)=A(I,J)+B(I,J)

CONTINUE

WRITE (6,3)((C(I,J),J=1,N),I=1,M

CALL EXIT

END

Multiplicación

Codificación (Programa)

DIMENSION A(3,3),B(3,3),C(3,3)

READ(5,10)M,N

FORMAT (2I2)

READ(5,1)((A(I.,J), J=1, N), I.=1, M)

FORMAT (3F3.0)

READ (5,1) ((B(I,J), J=1, N), I=1, M)

WRITE (6,2) ((A(I,J),J=1, N) ,I=1, M

WRITE (6,2) ((B(I,J),J=1, N) ,I=1, M

FORMAT (3(2X,F5.0))

Gauss-Jordan

DIAGRAMA DE FLUJO DE GAUSS JORDAN

INICIO

A,B,C,D [ 1..10, 1...20]

J,M,N,[INT]

MATRIZ CUADRADA?

1 AL 10

N=M+1 FOR I= 1 TO M DO

M>=1 FOR J = 1 TO N DO

N<=20

J<=0

ELEMENTOS DE MATRIZ

INDEPENDIENTE A [ I,J] FOR J= 1 TO N DO

A,,B,C,..[1,J],[2,J]....

GOTO XY (X+4,Y)

VALOR DE X1,X2,X3...

OTRO CALCULO(S/N)

Codificación(Programa)

DIMENSION A(20,21)

READ (5,10)N,M,E

FORMAT(2I3,E8.1)

NA=N+M

READ (5,11)((A (I,J),J=1NA),I=1, N)

FORMAT (7F5.0)

WRITE (6,12)((A(I,J),J=1,NA),I=1, N)

FORMAT (2X,7F5.0)

DET=1

DO 2 K=1, N

DET = DET* A (K,K)

IF ((ABS (A(K,K))).GT.E)GO TO 50

WRITE (6,13)

FORMAT (2X,*PIVOTE PE QUEÑO*)

GO TO 5

IP1=K+1

DO 3 J=IP1, NA

A(K,J)=A(K,J)/A(K,K)

CONTINUE

A(K,K)=1

DO 2 I=1, N

IF (I.EQ.K.OR.A(I,K).EQ.0.)GO TO 2

DO 4 J=IP1, NA

A(I,J)=A(I,J)-A(I,K)*A(K,J)

A(I,K)=0.

CONTINUE

WRITE(6,14)DTE

FORMAT (2X, 4HDET=,E15.8)

WRITE (6,15)((A(Di,J), J=1, NA), I=1, N)

FORMAT (7(2X,F8.5))

CALL EXIT

END

Gauss-Seidel

DIAGRAMA DE FLUJO DEL MÉTODO GAUSS - SEIDEL

INICIO

VAR: A:ARRAY [1...20],1...20]

I,J,K,N,X,Y,I2[ INT]

T,M,S,N [R]

NUMERO DE ECUACIONES

ELEM =9

FOR I = 1 TO N DO

GOTOXY (1,Y,X:3)

FOR J: = 1 TO N +1

GOTO XY READ [A[ I , J ] ]

IF J =N+1 IF A [I,O] IF A[ 1,1]=0

I = 0 I2=J

T=A[ I2,J]

T=[I,J]

IF I <>K

A[I,J]=A[1,J]- T*{K,J]

OTRO CALCULO (S/N)

FIN