Skip to main content

Gauss C++ (codigo)

Image
By 

% A es la matriz de coeficientes
% b es el vector de terminos independientes
% x es el vector solucion
% n es el orden del sistema

[n,m] = size(A);

% eliminacion
for i=1:n
  for k=i+1:n
    factor = (A(k,i)/A(i,i));
    A(k,i:n) = A(k,i:n) - factor*A(i,i:n);
    b(k) = b(k) - factor*b(i);
  end
end

%sustitucion inversa
x = zeros(1,n);
x(n) = b(n)/A(n,n);
for k=n-1:-1:1
  x(k) = (b(k) - sum(x(k+1:n).*A(k,k+1:n)))/A(k,k);
end
x

Comments

Post a Comment

Popular posts from this blog

Maquina de Mealy Ejemplo: Ascensor

By
entity behavioral is port ( Planta_Baja : in std_logic; Piso_Uno    : in std_logic; Piso_Dos    : in std_logic; salidas     : out std_logic_vector(1 downto 0) -- para polarizar el motor, tierra y voltaje por eso se require los dos bits ); end behavioral architecture AscMealy of behavioral is type ascEstados is (PB, 1, 2); signal Estados, prox_estados AscEstados; begin Process(Planta_Baja, Piso_Uno, Piso_Dos, Estado) begin Case PB => if (Planta_Baja = '1')then Prox_estado <= 1; else Prox_estado <= PB; endif; case 1 => if (Piso_Uno = '1')then Prox_estado <= 2; else prox_estado <= 1; endif; case 2 => if (Piso_Dos = '1')then Prox_estado <= 1; else prox_estado <= 2; endif; end Process; Process (reloj, Reset) begin if (Reset = '1')then estado <= PB; elsif (clk'event and clk = '1') then ...
Post a Comment
Read more

Llamar a 1 800 en USA desde Mexico

By
Cómo llamar números gratis en USA [toll free] desde México Ok, resulta que alguna vez los números gratis, se llamaban números 1- 800 (01-800 desde México). Pero lo malo de los números es que se acaban, Sobre todo en USA que les encanta vender todo. Ahora los Números gratis pueden ser 1 - 800, 1 - 877, 1 - 866 Para los americanos es cuestion de acostumbrarse, pero para nosotros los usuarios de telmex ¿cómo le hacemos si queremos marcar? Pues en primera pagando, los números gratis en USA usualmente no son gratis desde otros paises. Y en segunda, marcando correctamente. Muchos ya se saben el de los 800, es nada mas cuestion de cambiar el 01 - 800 por 001 - 881, pero y los 1 - 866 ? y los 1 - 877 ? Ahhh pues en mi otro blog ya tenía una entrada con una tablita que lo explica: Si el numero comienza con 1 888 se marca: 001 881 Si el numero comienza con 1 877 se marca: 001 882 Si el numero comienza con 1 866 se marca: 001 883 ¿Por qué? No tengo la más mínima idea, supongo que...
Post a Comment
Read more

Maquina Moore Ejemplo: Ascensor

By
entity behavioral is port ( Planta_Baja : in std_logic; Piso_Uno    : in std_logic; Piso_Dos    : in std_logic; salidas     : out std_logic_vector(1 downto 0) -- para polarizar el motor, tierra y voltaje por eso se require los dos bits ); end behavioral architecture AscMoore of behavioral is type ascEstados is (PB, 1, 2); signal Estados, prox_estados AscEstados; begin Process(init, PB_UP, P1_UP, P1_DOwN, P2, P2_DOWN) signal estado, prox_edo : AscMoore; begin Process (estado,Planta_Baja, Planta_Un, Planta_Dos) begin case estado is when init => if Planta_Baja='1' then Prox_estado <= PB_UP else Prox_estado <=init; endif; when PB_UP => prox_estado <= P1_UP; when P1_UP => if Piso_Uno='1' then Prox_estado <= P1_UP; else Prox_estado <= P1_DOWN; endif; when P1_DOWN => Prox_estado <= init; Process (estado) ...
Post a Comment
Read more