## 14 Badges

15 years, 323 days

## Notation...

If you expect 1+y'=0 as result, you have to explicitly describe y as dependent on x:

x+y(x)=1;

By the way 1=0! is correct :o)

## don't know...

Well I'm not a developer, I can only speculate. Maybe it is because the falling bar in MapleSim  has its masses concentrated at its ends and the simulink model assumes otherwise.

Do you have a barbell in reach, so that you could perform some falling experiments?

:-)

## Damping...

I think one has to find the proper values for damping and the elastic constant. Also, I added a horizontal reactive force to the ground of 2.msim, as you can see in the custom component. Show you roll of wrapping paper model.

2.msim

## Falling rod...

Hi,
here an example how to model the falling rod. It can be easily extended to more complex examples in 3D, like cube, pyramid or something, just add some more mass points with their own attached subsystem, which checks if the position of the rigid body component drops below ground and applies an reactive force to it, if required.

1.msim

Here the model of a rolling dice:

http://www.maplesoft.com/products/maplesim/ModelGallery/detail.aspx?id=123

## Prismatic and Force Component...

How about this? There is a signal force step in this model that works on a prismatic. After 2 seconds, the force goes to zero, and the "rope" is released.

Pendulum_with_Releas.msim

## Your worksheet opens fine here (Maple 16...

Your worksheet opens fine here (Maple 16.02). There is a recursive assignment error in there, but all is readable.
Maybe a cold start would just do the trick.

## example...

restart:
f:=piecewise(x^2+y^2<3,y-x^2,undefined);
plots[implicitplot](f,x=-2..2,y=-2..2,scaling=constrained,view=[-2..2,-2..2]);

## Formal Power Series Option for Convert...

convert((1+x)^n,FormalPowerSeries,x);

## plot3D...

restart:c:=50;beta:=10;A:=0.01;L:=0.5;
y:= A*exp(-beta*(x-L-c*t)^2/(L^2));
plot3d(y,x=0..1.5,t=0..0.003,axes=framed,grid=[80,80]);

## keyboard shortcuts...

Use Strg+J to insert a new execution group after cursor (Strg+K before cursor). Just type ?shortcuts to get an overview and select the help page you want from the list.

overview of packages: ?index,package

example to check content: with(LinearAlgebra);

you can ask for specific help pages to a command by entering it with a leading question mark, like ?inifcns

## DEPlot...

In case you don't want or can't solve the DE explicitly, use DEplot:

restart;
k:=2:M:=10*exp(-t/5)+3*(1-exp(-t/5)):
DETools[DEplot]( diff(y(t),t)=k*(1-y(t)/M)*y(t), y(t), t=0..20, [y(0)=1] );

## Try this...

There are surely different ways to do this. start with creating the sequence. My starting point: 7/9 =0.777...
Then multiply that with 10^n and cut off the decimal places.

A:=proc(n) local m; convert([seq(floor(7/9*10^m),m=1..n)],`+`);  end proc;

## Another way...

collect(expand(A*cos(w*t+p1)+B*cos(w*t+p2)),[sin(w*t),cos(w*t)]);
convert(%,phaseamp,t);

## solve numerically...

Note that not all equations can be solved explicitly. Maple returns RootOfs in this case. See my sheet for a numerical solution.

Max_difference_of_vo.mw

## difference of volumes -> max...

You already have two explicit formulas V(t) for both cases. Now you have to find the maximum of V1(t)-V2(t).

So calculate the first derivative of it, then...

 1 2 3 4 5 6 7 Last Page 1 of 16
﻿