I think that the option should be added to the plotsetup command as a plotoptions for the default plot device. The reason is that all the commands that produce plots actually produce an algebraic structure that should be device independent; device specific parameters should not be stored in that structure.

@emma hassan Are you saying that you want the numbers themselves actually printed on the plot?

@emma hassan Are you saying that you want the numbers themselves actually printed on the plot?

@maxou Once you take the derivative with respect to a variable appearing in exponents (rho in this case), the expression becomes vastly complicated. I have no idea on how to approach it.

@maxou Once you take the derivative with respect to a variable appearing in exponents (rho in this case), the expression becomes vastly complicated. I have no idea on how to approach it.

@Michael_Ormstrup To expand on it, I need to know what a and b represent. I would've thought that they were the plot range, but you don't use them in the plot command in your procedure.

@Michael_Ormstrup To expand on it, I need to know what a and b represent. I would've thought that they were the plot range, but you don't use them in the plot command in your procedure.

@spradlig Can you show an example of a simple one-line function on which D doesn't work? The unapply solution will only work if the application of the function to a generic argument results in an algebraic expression. It's a somewhat specialized situation.

@spradlig Can you show an example of a simple one-line function on which D doesn't work? The unapply solution will only work if the application of the function to a generic argument results in an algebraic expression. It's a somewhat specialized situation.

How do you define the eigenvalues of a system of nonlinear equations?

@jschulzb You may be able to save memory by not using parameters (assign their values in the original dsolve), and using the output= Array form of output. With this form you specify the values of the independent value at which the system is to be evaluated and dsolve returns an Array with those evaluations. This should have an effect similar to interpolate= false. See ?dsolve,numeric .

@goli A better way to correct the second plot is to get better spacing of the points by using logarithmic scaling on the parameter:

plot([eval(ns, t= 10^t), eval(r, t= 10^t), t= 11..17], 0.8..1.2, 0..1);

@goli A better way to correct the second plot is to get better spacing of the points by using logarithmic scaling on the parameter:

plot([eval(ns, t= 10^t), eval(r, t= 10^t), t= 11..17], 0.8..1.2, 0..1);

@goli By redisplaying the plots with style= point, you can see that in the second plot, the points are all bunched up at the bottom. If you redo the second plot with numpoints= 21000, then it will look like the first plot.

You are almost right that all the points are inside the view. In each plot, there is one point outside the view.

Download_13.3.mw

@goli By redisplaying the plots with style= point, you can see that in the second plot, the points are all bunched up at the bottom. If you redo the second plot with numpoints= 21000, then it will look like the first plot.

You are almost right that all the points are inside the view. In each plot, there is one point outside the view.

Download 13.3.mw