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
3D Studio Max 2013
Espai
3D Studio MAX es una de las plataformas 3D más extendidas para la infografía profesional. Actualmente es utilizada...
Solicita InformaciÓn


publicidad

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