Расхождение результатов в GUI

function btn_Callback(hObject, eventdata, handles)
% hObject    handle to btn (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
m=str2num(get(handles.edit1,'String'));
m0=str2num(get(handles.edit2,'String'));
I=str2num(get(handles.edit3,'String'));
R=str2num(get(handles.edit4,'String'));
z=str2num(get(handles.edit5,'String'));
N=str2num(get(handles.edit6,'String'));
B = ((m*m0*I*R^2)./(2*((R^2+z.^2).^1.5)))*N; % вычисление магнитной индукции
deltaz = 0.1;
ind_z = find (z == deltaz) - 1 % находим номер элемента массива
i =1;
while i <= length(B) - ind_z % объявляем предел массива В
    x(i)=-( B(i + ind_z) - B(i) )./deltaz; % вычисляем гардиент магнитного поля
    y(i) = z(i);
    i = i + 1;
end
plot(y,x) % выводим график
xlabel('z, м'); % именуем ось абсцисс 
ylabel('dB/dz, Тл/м'); % именуем ось ординат 
grid on % выводим сетку
axis % масштабируем оси графика

m = 5*10^3; % магнитная проницаемость материала сердечника 
m0 =4*pi*10^(-7); % коэффициент магнитной проницаемости
I = 1.5; % сила тока
R = 0.15;% радиус соленоида
z = 0:0.01:0.7 % диапазон значений по оси z
N=200; %количество витков
B = ((m*m0*I*R^2)./(2*((R^2+z.^2).^1.5)))*N; % вычисление магнитной индукции
deltaz = 0.1;
ind_z = find (z == deltaz) - 1 % находим номер элемента массива
i =1;
while i <= length(B) - ind_z % объявляем предел массива В
    x(i)=-( B(i + ind_z) - B(i) )./deltaz; % вычисляем гардиент магнитного поля
    y(i) = z(i);
    i = i + 1;
end
plot(y,x) % выводим график
xlabel('z, м'); % именуем ось абсцисс 
ylabel('dB/dz, Тл/м'); % именуем ось ординат 
grid on % выводим сетку
axis % масштабируем оси графика    
B
x