Создание графика по формуле

plot(x,y)

clear;
clc;
 
B = [1.2 0.66 0.46 0.28 0.24 0.2] ; 
z = [0 0.1 0.2 0.3 0.4 0.5 0.6]; 
z1 = 0.1;
 
ind_z = find (z == z1) - 1
i = 1;
while i <= length(B) - ind_z
    x(i) = ( B(i + ind_z) - B(i) )./z1;
    y(i) = z(i);
    i = i + 1;
    
end
 
plot(y,x)

m = 10^5;
m0 = 0.000001253;
I = 10;
R = 0.1;
z = 0:0.1:10 % диапазон значений z
B = ( m*m0*I*R^2 ). / (2*((R^2+z.^2).^1.5) );
B % чтобы вывести значения

m = 10^5;
m0 = 0.000001253;
I = 10;
R = 0.1;
z = 0:0.1:0.6 % диапазон значений z
B = ( m*m0*I*R^2 ). / (2*((R^2+z.^2).^1.5) );
z1 = 0.1;
ind_z = find (z == z1) - 1
i = 1;
while i <= length(B) - ind_z
    x(i) = ( B(i + ind_z) - B(i) )./z1;
    y(i) = z(i);
    i = i + 1;
end
plot(y,x)
xlabel('z, м');
ylabel('dB/dz, Тл/м');
grid on
x

A^2

A.^2