Tutorial

ylims=solve(z2-z1,y) solve. has found four solutions instead of two. This should not surprise us too much since the equation in question is quartic. If we look carefully, we can see that the last two solutions involve square roots of negative quantities, so that the limits we want are the first two. ................
................