How to plot inseparable variables.

How to plot the following equation for values of d ranging from 0 to 0.04

m=0.025, w=10, k=8*10^-9. I have obtained the following graph

The required graph is as follows

The code is as follows

function answer

%solve implicit function and plot

npoints = 40;

d = linspace(0,0.04,npoints);

x = zeros(1,npoints);

guess = 0.04;

for i = 1:npoints

x(i) = fzero(@equation,guess,[],d(i));

guess = x(i);

end

plot(d,x)

xlabel('d')

ylabel('x')

function dy = equation(x,d)

dy = 2.5*x-(2.4*10^-8*x)/((x^2+d^2)^2.5);

Answer 1

[VBBV]

npoints = 50;
d = linspace(0,0.04,npoints);
x = zeros(1,npoints);
equation = @(x,d) 2.5*x-(2.4*10^-8*x)/((x^2+d^2)^2.5);
guess = -0.04;
for i = 1:npoints
x(i) = fzero(equation,guess,[],d(i));
guess = x(i);
end;
guess = 0.04;
for i = 1:npoints
xx(i) = fzero(equation,guess,[],d(i));
guess = xx(i);
end
plot(d,x,d,xx,'-b')
xlabel('d')
ylabel('x')

Use positive and negative values