reflect callback

figure(1);
            epsilon0 = 1/36/pi*10^(-9);
            mu0 = 4*pi*10^(-7);
            Epsilon = app.EpsilonEdit.Value*epsilon0;
            Epsilon1 = app.Epsilon1Edit.Value*epsilon0;
            Epsilon2 = app.Epsilon2Edit.Value*epsilon0;
            Mu = app.MuEdit.Value*mu0;
            Mu1 = app.Mu1Edit.Value*mu0;
            Mu2 = app.Mu2Edit.Value*mu0;
            Sigma = app.SigmaEdit.Value;
            Sigma1 = app.Sigma1Edit.Value;
            Sigma2 = app.Sigma2Edit.Value;
            Omega = 2*pi*app.FreqEdit.Value;
            phi =0:pi/500:pi/2;
            EpsilonComplex = Epsilon*(1-1i*Sigma/(Omega*Epsilon));
            Epsilon1Complex = Epsilon1*(1-1i*Sigma1/(Omega*Epsilon1));
            Epsilon2Complex = Epsilon2*(1-1i*Sigma2/(Omega*Epsilon2));
            k = Omega*sqrt(EpsilonComplex*Mu);
            k1 = Omega*sqrt(Epsilon1Complex*Mu1);
            k2 = Omega*sqrt(Epsilon2Complex*Mu2);
            phipr1 = asin((k/k1)*sin(phi));
            phipr2 = asin((k/k2)*sin(phi));
            W = sqrt(Mu/EpsilonComplex);
            W1 = sqrt(Mu1/Epsilon1Complex);
            W2 = sqrt(Mu2/Epsilon2Complex);


            if app.HorizPol.Value %normal
                R1 = (W1*cos(phi) - W*cos(phipr1)) ./ (W1*cos(phi) + W*cos(phipr1));
                R2 = (W2*cos(phi) - W*cos(phipr2)) ./ (W2*cos(phi) + W*cos(phipr2))
            end;
            if app.VertPol.Value %parall
                R1 = (W*cos(phi) - W1*cos(phipr1)) ./ (W*cos(phi) + W1*cos(phipr1));
                R2 = (W*cos(phi) - W2*cos(phipr2)) ./ (W*cos(phi) + W2*cos(phipr2));
            end;
            figure(1);
            plot(phi, abs(R1), phi, abs(R2), 'LineWidth',2);
            legend('-->1','-->2');
            figure(2);
            plot(phi, abs(angle(R1)), phi, abs(angle(R2)));