Thank you! Now it works great! Always the little things...

Just two last small questions:

- How do I find the maximum value in a function? Currently I just plot it and measure the peak, but there must be a better way. Just using something like max or extrema doens't really work.

- (More important): how do I plot a derative? I need to plot the moment, to obtain the maximum moment. The moment is defined as the derative of theta (I believe). This does not work:

> sol3:-plot(E*J*(diff(theta(x, t), x)), t = 0.188e-1, labels = ['x*[m]', 'M*[Nm]'], title = 'total*moment');

It gives me an empty graph. Changing diff(theta(x, t), x) to D[1](theta)(x, t) does not make a difference. How do I plot a derative?

shear_beam_pin_pin_t.mw

Thank you! Now it works great! Always the little things...

Just two last small questions:

- How do I find the maximum value in a function? Currently I just plot it and measure the peak, but there must be a better way. Just using something like max or extrema doens't really work.

- (More important): how do I plot a derative? I need to plot the moment, to obtain the maximum moment. The moment is defined as the derative of theta (I believe). This does not work:

> sol3:-plot(E*J*(diff(theta(x, t), x)), t = 0.188e-1, labels = ['x*[m]', 'M*[Nm]'], title = 'total*moment');

It gives me an empty graph. Changing diff(theta(x, t), x) to D[1](theta)(x, t) does not make a difference. How do I plot a derative?

shear_beam_pin_pin_t.mw