-
Notifications
You must be signed in to change notification settings - Fork 4
/
platesD.sci
39 lines (33 loc) · 906 Bytes
/
platesD.sci
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
// Una función que crea con filas arbitrarias
// la matriz A de coeficientes y el vector b
// filas: El número de platos
// Vm: El multiplicador a la corriente de gas
// Lm: EL multiplicador a la corriente de aceite
function [A, b] = platesD(filas, Vm, Lm)
// Constantes
V=100*Vm
L=500*Lm
y0=0.09
x0=0
m=0.12
// Vector b
b = zeros(filas,1)
b(1) = V*y0
// Matriz de coeficientes
_A = zeros(filas, filas+2)
// La fila que forma el resto de filas
fila = [-V*m L+m*V -L]
fila3 = [-V*m L*1.25+m*V -L*1.25] // La fila 3 es especial
// Repitiendo la misma fila pero con un padding de 0s
for i = 1:filas
if(i == 3)
filaExt = [zeros(1, i-1) fila3 zeros(1, filas+2-3-(i-1))]
else
filaExt = [zeros(1, i-1) fila zeros(1, filas+2-3-(i-1))]
end
_A(i,:) = filaExt
end
// En la primera fila el -V*m no existe y en la última fila
// el -L tampoco.
A = _A(:,2:filas+1)
endfunction