embed

for i=1:length(Bmec)
    if abs(Bmec(i))>=2 
        for j=1:(length(Bfou))
            if abs(Bfou(j))>=2
                while abs((Bfou(j)-Bmec(i))./Bmec(i)) > 0.05 %hata payi yuksek
                    g=g+g*0.1;%hava araligini arttir %10
                    [reg]=Fourier_open_circuit(g,delta_z,Rout_stat,Rin_stat,La,h_mag,h_rot_tooth,ts_y,Nr,Ns,rpm,mu_mr,Br,N,I);
                    Bfou=reg(6).Bp;
                end
            elseif (abs(Bfou(j)-Bmec(i))./Bmec(i)<=0.05) %hata payi uygun if dongusunden cikip bir sonraki Bmec degerine bak