Modelos discretos ejercicios

[pic 1]

Determinar

1. Los kilogramos necesarios de cada material para preparar la formula
2. El costo total de la formula alimenticia.

X1= # Kgr de Maiz, X2= # Kgr de Desperdicio, X3= # Kgr de Alfafa, X4= # Kgr de Cebada

80*X1 + 15*X2 + 35*X3 + 60*X4 = 230 (Ec1)

28*X1 + 72*X2 + 57*X3 + 25*X4= 180 (Ec2)

20*X1 + 20*X2 + 12*X3 + 20*X4 = 80 (Ec3)

50*X1 + 10*X2 + 20*X3 + 60*X4= 160 (Ec4)

Proceso GaussJordan

        Dimension V(100,100);

        Dimension V1(100);

              Dimension C1(100);

        Definir T como real;

        Escribir  "Ingrese el número de componentes de la formula (n):";

        leer n;

              para i<-1 hasta n con paso 1 hacer

                escribir " Costo ", i, ":";

                Leer C1(i);

        FinPara

        escribir "Ingrese la incógnita del componente ";

        para i<-1 hasta n con paso 1 hacer

                escribir " de la variable ", i, ":";

                Leer V1(i);

        FinPara

        Escribir "Ingrese el termino independiente de la mezcla ";

        Para i<-1 Hasta n Con Paso 1 Hacer

                Para j<-1 Hasta n Con Paso 1 Hacer

                        Escribir " la variable ", V1(j), " de la ecuacion ", i, ":";

                        Leer V(i,j);

                FinPara

                Escribir "ingrese la constante de la ecuacion ", i, ":";

                Leer V(i,n+1);

        FinPara

        Para i<-1 Hasta n Con Paso 1 Hacer

                si V(i,i)=0.0 Entonces

                        Escribir "error: división entre cero.";

                FinSi

                Para k<-1 hasta n Con Paso 1 Hacer

                        si no (k=i) Entonces

                                x<-V(k,i)/V(i,i);

                                para j<-1 Hasta n+1 Con Paso 1 Hacer

                                        V(k,j)<-V(k,j)-x*V(i,j);

                                FinPara

                        FinSi

                FinPara

        FinPara

        Escribir "la solucion al sistema de ecuaciones es:";

              T<-0;

        Para i<-1 Hasta n Con Paso 1 Hacer

                x<-V(i,n+1)/V(i,i);

                Escribir V1(i), " = ", x;

                            T<-T+ x*C1(i);

        FinPara

            Escribir "Costo Total:", T;        

FinProceso

Gauss Seidel

        Algoritmo sin_titulo

        Escribir 'Ingrese el error relativo deseado';

        Leer ErD; // ingresar 1 que corresponde al 1%

        X1 <- 0;

        X2 <- 0;

        X3 <- 0;

        X4 <- 0;

        //X5 <- 0;

        Cont <-0;

        Repetir

                X1a <-X1;

                Escribir X1a;

                X2a <-X2;

                Escribir X2a;

                X3a <-X3;

                Escribir X3a;

                X4a <-X4;

                Escribir X4a;

                X1 <- (230 - 15*x2 - 35*X3 - 60*X4)/80;

                Escribir X1;

                X2 <- (180 - 28*X1 - 57*X3 - 25*X4)/72;

                Escribir X2;

                X3 <- (160 - 50*X1 - 10*X2 - 60*X4)/20;

                Escribir X3;

                X4 <- (80 - 20*X1 - 20*X2 - 12*X3)/20;

                Escribir X4;

                ErrRelX1 <- Abs(X1-X1a);

                Escribir ErrRelX1;

                ErrRelX2 <- Abs(X2-X2a);

                Escribir ErrRelX2;

                ErrRelX3 <- Abs(X3-X3a);

                Escribir ErrRelX3;

                ErrRelX4 <- Abs(X4-X4a);

                Escribir ErrRelX4;

                Cont <- Cont +1;

        Hasta Que (ErrRelX1 < ErD) & (ErrRelX2 < ErD) & (ErrRelX3 < ErD)& (ErrRelX4 < ErD) //(ErrRelX5 < ErD)

...