*sigh* Mathematica,

how do I plot the bifurcation diagram with Maple of the function:

dx/dt = rx - sin(x) ?

BIFURCATION DIAGRAM PROCEDURE
restart:
with(plots):
*xexpr is the logistic function to be iterated (we always start off at x=1/2, which will eventually attract).
*[ra,rb] is the range of the parameter.
*acc is the number of points sampled in [ra,rb]
Bifurcation := proc(initialpoint,xexpr,ra,rb,acc)
   local p1,hr,A,L1,i,j,phi:
   global r,L2:
   hr := unapply(xexpr,x);
   A := Vector(600):
   L1 := Vector(acc*500):

The LinearAlgebra`MatrixManipulation` package is now obsolete with the new version of Mathematica.

Say you had an existing script with commands:

TakeRows, TakeColumns, AppendColumns, ZeroMatrix, AppendRows (which are also obsolete)

What do you change these commands to?

Hello,

if my Fourier coefficient b_n := 4/(n*pi) which is only true for n odd, how do I program this in Maple?

For[n=0,n<1000,n++,

   If[Mod[2, n],

      c = b_n Sin[n w_0 x]

]

That was my attempt, it's wrong of course. What do I need to add/change?

Thanks!