2740 Reputation

12 Badges

17 years, 47 days

Dr. Robert J. Lopez, Emeritus Professor of Mathematics at the Rose-Hulman Institute of Technology in Terre Haute, Indiana, USA, is an award winning educator in mathematics and is the author of several books including Advanced Engineering Mathematics (Addison-Wesley 2001). For over two decades, Dr. Lopez has also been a visionary figure in the introduction of Maplesoft technology into undergraduate education. Dr. Lopez earned his Ph.D. in mathematics from Purdue University, his MS from the University of Missouri - Rolla, and his BA from Marist College. He has held academic appointments at Rose-Hulman (1985-2003), Memorial University of Newfoundland (1973-1985), and the University of Nebraska - Lincoln (1970-1973). His publication and research history includes manuscripts and papers in a variety of pure and applied mathematics topics. He has received numerous awards for outstanding scholarship and teaching.

MaplePrimes Activity

These are replies submitted by rlopez

It appears that the Enter key was pressed after each of the three lines of the do-loop. The loop has to be in one single execution group, and this is done by use of Shift-Enter, not just Enter.

There was one more trig function without parentheses. Once I fixed that last guy, r-vec display properly.

There are syntax errors: all trig functions require round parentheses. It is an error to write sin x. It must always be sin(x).

Also, the red error message tells you about another error: The variable r has been assigned to. Therefore, it is no longer a free name, so you can't differentiate with respect to r.

I don't understand why the vector r_vec is printing as a column vector. I tried <x,y,z> and that prints in Basis Format with basis vectors displayed. When I execute About(r_vec) I see that it is being defined as a Matrix, not a vector. I don't understand why.

I tyically work in the Student version of VC because it is more forgiving with respect to defining and redefining coordinate systems. In that package, r_vec displayed in Basis Format, not as a column vector. That's as much as I've discovered so far.

@gmzsvsclk Please define what you mean by "the points (t,y)"? Which points? I gave you a method for finding those points on a graph that Maple could be  made to draw. Obviously, you didn't want that collection of points, so which points do you want? You have also been given recipes for finding arbitrary points by solving explicitly for y. So, which points do you want?

For a single variable, Newton's method takes x[n+1] as the x-intercept on the line tangent to the graph of f(x) at x[n]. The Secant method replaces this tangent line with a secant line drawn through x[n-1] and x[n].

In the multivariate case, the Jacobian matrix in Newton's method can be replaced with alternatives, such as in Broyden's method. I'm not sure that such methods follow from a simple geometry that generalizes a secant line. In fact, Newton's method in the multivariate case is derived via Taylor expansions of each function, and not through any obvious geometry. If it can be so based, I've never seen it, and would appreciate a pointer to it.


The Roots command will succeed if Digits is set to at least 17. As several responses indicate, the calculation is delicate and subject to roundoff. I guess I missed that the first time I looked. Thanks to all who provided a sharper analysis.


Delete the number=1 option, and don't specify a derivative condition as an initial condition for a first-order ODE.

With the DEtools package installed, the following syntax works:

phaseportrait(X(t) = diff(X(t), t), X(t), t = -1 .. 2.5,[[X(0)=1],[X(0)=-1]])

but not phaseportrait(X(t) = diff(X(t), t), X(t), t = -1 .. 2.5,[[X(0)=1],[D(X)(0)=-1]])

The dsolve command can obtain the solution with the derivative IC, which can be graphed separately, then joined to the plaseportrait via the display command in the plots package.

p1:=DEtools:-phaseportrait(X(t) = diff(X(t), t), X(t), t = -1 .. 2.5,[[X(0)=1]]):

@brian bovril 

Nice use of Maple to find a general pattern. I have to admit that I scratched around a bit before I came to the same conclusion. But I did end up using rsolve along the way.

I don't understand your statement that the loop should finish when the computed T is the same as the initial T. If it's the initial T that you want, why do any calculations? You already have it.

It would make more sense to terminate the loop when two successive computations of T are essentially the same.

Please clarify.

@Joe Riel 

The SpaceCurve command in either of the VectorCalculus packages will graph a curve defined parametrically by either a 2- or 3-dimensional vector. No need to switch graphing commands for different dimensions.

Have you followed the thread most recently appearing here:




The Physics:-Vectors package has some neat functionality for vectors and differentiation. Unfortunately, it has none for visualization or integration. The Student VectorCalculus package has wonderful functionality for integration and visualization.

The Physics:-Vectors package is the only place in Maple where one can manipulate "unprojected" vectors, that is, vectors that do not have to be given components. So, vector identities can be manipulated without having to specify components for the vectors.

The Student VectorCalculus package is far more lenient with respect to setting coordinates. The VectorCalculus package demands at every turn a declaration of the coordinate system and the coordinate-variable names. The Student package, recognizing only Cartesian, polar, cylindrical, and spherical coordinates, all with default coordinate-variable names, is much more user-friendly for someone not working in esoteric systems.

@Markiyan Hirnyk 

Interesting. I have the source code for Maple V Release 4 on a set of 10 floppy disks, the original media used for the distribution of this version. My Windows 10 laptop does not have a floppy-disk reader, nor does my Windows7 desktop. How did you manage to install MVR4 on a machine capable of running Windows 10?


Minimize in the Optimization package computes numerically. It is not able to work with the c's as parameters, and that accounts for the zero solution. Hence, the LagrangeMultipliers command, which can work symbolically.

@mahmood1800 Originally, the constraint was an equation. Now, you show that not to be the case. I'm therefore confused as to what your original problem actually was.

3 4 5 6 7 8 9 Last Page 5 of 18