clear all; clc; format long;

dolni_vnitrni = 1; %buzení vnitřního vinutí dolního deflektoru (Az)
dolni_vnejsi = 1; %buzení vnějšího vinutí dolního deflektoru (Az)

vychylka_x = 2; %požadovaná výchylka svazku ve směru x
vychylka_y = 2; %požadovaná výchylka svazku ve směru y

horni_vnitrni_1 = -0.897819076*dolni_vnitrni; %buzení vnitřního vinutí horního deflektoru v prvním kroku(Az)
horni_vnejsi_1 = 0.314655821*dolni_vnitrni; %buzení vnějšího vinutí horního deflektoru v prvním kroku (Az)
gamma_vnitrni = 1.0699591E-01 -3.0951453E-02j; %vychýlení v prvním kroku
gamma_vnitrni_real = real(gamma_vnitrni); %vychýlení v prvním kroku ve směru x (mm)
gamma_vnitrni_imag = imag(gamma_vnitrni); %vychýlení v prvním kroku ve směru y (mm)

horni_vnitrni_2 = -0.248364845*dolni_vnejsi; %buzení vnitřního vinutí horního deflektoru v prvním kroku(Az)
horni_vnejsi_2 = -0.898919383*dolni_vnejsi; %buzení vnějšího vinutí horního deflektoru v prvním kroku (Az)
gamma_vnejsi =  2.9076359E-02 +9.5537210E-02j; %vychýlení v druhém kroku
gamma_vnejsi_real = real(gamma_vnejsi); %vychýlení v druhém kroku ve směru x (mm)
gamma_vnejsi_imag = imag(gamma_vnejsi); %vychýlení v druhém kroku ve směru y (mm)

A = [gamma_vnitrni_real gamma_vnejsi_real; gamma_vnitrni_imag gamma_vnejsi_imag];
B = [vychylka_x; vychylka_y];

C=A\B;
magnituda_dolni_vnitrni = C(1,1) %potřebné buzení vnitřního vinutí dolního deflektoru pro minimalizaci komy (Az)
magnituda_dolni_vnejsi = C(2,1) %potřebné buzení vnějšího vinutí dolního deflektoru pro minimalizaci komy (Az)

magnituda_horniho_vnitrniho = horni_vnitrni_1*magnituda_dolni_vnitrni + horni_vnitrni_2*magnituda_dolni_vnejsi %potřebné buzení vnitřního vinutí horního deflektoru pro minimalizaci komy (Az)
magnituda_horniho_vnejsiho = horni_vnejsi_1*magnituda_dolni_vnitrni + horni_vnejsi_2*magnituda_dolni_vnejsi %potřebné buzení vnějšího vinutí horního deflektoru pro minimalizaci komy (Az)