1. Solución de línea de comandos
syms abcyx
[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[x,y])
2. Resuelve en el editor
borrar
syms A1 A2 B1 B2 r1 r2 t41 t42
k = 0.01
n = 1;
para d = 0: 300
para q =
1:10 eq1 = A1 + A2 + B1 + B2-sqrt (2) * (1 + r1);
eq2 = A1-A2 + B1-B2-sqrt (2) * (1 + r2);
eq3 = 1 + r1-t41;
eq4 = 1 + r2-t42;
eq5 = A1 * exp (i * (qk) * d) + A2 * exp (i * (q + k) * d) + B1 * exp (-i * (q + k) * d) + B2 * exp ( -i * (qk) * d);
eq6 = A1 * exp (i * (qk) * d) -A2 * exp (i * (q + k) * d) + B1 * exp (-i * (q + k) * d) -B2 * exp ( -i * (qk) * d);
eq7 = t41 * (kq) + r1 * (q + k) + (A1 * (kq) -A2 * (q + k) + B1 * (q + k) + B2 * (qk)) / sqrt (2) + (kq) * t41;
eq8 = t42 * (- kq) + r2 * (qk) + (A1 * (kq) + A2 * (q + k) + B1 * (q + k) -B2 * (qk)) / sqrt (2) - (k + q) * t42;
[A01 (n), A02 (n), B01 (n), B02 (n), r01 (n), r02 (n), t041 (n), t042 (n)] = resolver (eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, A1, A2, B1, B2, r1, r2, t41, t42);
n = n + 1;
fin
fin
disp ('A01 A02 B01 B02 r01 r02 t041 t042')
X = [doble (A01); doble (A02); doble (B01); doble (B02); doble (r01); doble (r02); doble (t041) ; doble (t042)] '