dlambda1dmuX2mpl2
# FILE dlambda1dmuX2mpl2.txt
phi2 := z-> arctan(1/sqrt(z))/sqrt(z);
plot([phi2(z), 1/(z+1)], z = .125 .. 8, color = [red, green])
#
# Next bit for N=2 Theorem 3
factor(diff(phi2(z), z)+phi2(z)/(2*z));
dmuYdmuX := simplify(solve(diff(phi2(muX)*phi2(muY(muX)), muX) = 0, diff(muY(muX), muX)));
dlambda1dmuX:= simplify(solve(diff(phi2(muX)*phi2(lambda1(muX)-muX), muX) = 0,
diff(lambda1(muX), muX)));
with(plots, implicitplot, display, odeplot);
ans := dsolve([diff(lambda1(muX), muX) = dlambda1dmuX, lambda1(1) = 2],
lambda1(muX), range = 1 .. 10, type = numeric);
odeplot(ans);