Thomas Richard

Mr. Thomas Richard

3255 Reputation

13 Badges

15 years, 96 days
Maplesoft Europe GmbH
Technical professional in industry or government
Aachen, North Rhine-Westphalia, Germany

MaplePrimes Activity


These are answers submitted by Thomas Richard

Maple 2016 returns 0 for both integrals.

If you don't want to upload your files to MaplePrimes, please send them to support@maplesoft.com. Otherwise the problem cannot be reproduced, I'm afraid.

Please include version numbers of Excel, Maple, Windows. And wordsize - since Office 2010, there are both 32-bit and 64-bit Excel variants.

It's probably best to output to EPS and convert that to PDF (via GhostScript / GSview e.g.).

That export can also be done programmatically, i.e. non-interactively. Some code to get you started:

currentdir(kernelopts(homedir));
plotsetup(cps,plotoutput=`myplot.eps`,plotoptions=`landscape`):
plot(cos(z^2),z=0..5);
plotsetup(default); # to switch back to screen output
plot(cos(z^2),z=0..5);

Please see ?plotsetup, ?plot,device and ?plot,ps for more details and options.

There is a colon missing in your assignment of [eqd1, eqd2, eqd3, l(0) = 50, g(0) = 50, s(0) = 0.02] to sysdif. Use the := operator for that.

A single = character forms an equation, so Maple won't issue any warning message here.

After fixing that typo, you can obtain plots such as

plots:-odeplot(soleqd,[[t,g(t)],[t,l(t)],[t,s(t)]],t=0..0.002,legend=["g","l","s"]);

For me, Maple 2016 was added automatically in the correct position (right under Maple 2015). I suggest that you send your install log file to support@maplesoft.com.

Do you have any desktop icon / taskbar icon for Maple 2016?

This can be achieved by simplify with the 'symbolic' option (see ?simplify,details):

simplify(f,'symbolic');

Please see the docs under ?odetest: for verifying series solutions, you will need to add the 'series' option. It will then return FAIL and a warning message if you supply an ODE only. This is because it needs initial conditions. Simplest case:

ic1 := T(0) = 0:
ic2 := D(T)(0) = 0:
ivp := [ode, ic1, ic2]:
Sol__ivp := dsolve(ivp, T(r), 'series');
odetest(Sol__ivp, ivp, 'series');

What did you enter in Maple? Here's my input (skipped trailing zeros, and simplified exponents a bit):

tf_original := (0.1299e-4*(17.174*Omega^2-1.5702e6))*Omega^2/(-196.12708*Omega^4+3.95412129e7*Omega^2-1.8771741e12);
tf := subs(Omega=s, tf_original);
with(DynamicSystems):
sys := TransferFunction(tf);
BodePlot(sys);


 

That's not possible, for several technical reasons. Unlike the Maplets system, Embedded Components are under active development, so that's your best alternative. Please see ?EmbeddedComponents for a start. In the Plotter component, all of Maple's plot toolbar is available.

It should be possible, albeit with some effort. Check out e.g. this cam model in our Model Gallery. But as you can see, it requires programming a custom component in Modelica, and the STL geometry was created in Maple. So you wouldn't win much, I'd say. In short, MapleSim is designed for physical modeling, whereas Maple is the more natural environment for such geometrically motivated animations (which I wouldn't really call simulations, but that's debatable).

I suppose the equations come from MapleSim.  There is a command for the Jacobian w.r.t. the generalized coordinates. Please see ?Kinematic_Exports. You can easily add this call to the Multibody Analysis template:

MB:-xJacPos();

To see the position constraints, modify the first line in section "Multibody Analysis", replacing the trailing colon by a semicolon:

vPosCons := MB:-GetPosCons();

Not exactly what you had asked for, but perhaps it's useful as well:

length(eq_liaison);
sel := simplify(eq_liaison,'size');
length(sel);

AFAIK, using VS 2015 requires Maple 2015 with the latest update 2015.2.

It's pretty straightforward, but the result is monstrous, so Maple won't display it by default:

eq1 := b1 = a*x + b*y + c*z:
eq2 := b2 = d*x^2 + e*y^2 + f*z^2:
eq3 := b3 = g*x^3 + h*y^3 + k*z^3:
sys := [eq1,eq2,eq3]:
vars := [x,y,z]:
sol := solve(sys,vars);

If you really want to see the result, uncheck "Tools > Options > Precision > Limit expression length to" and do the solve call again. You will notice that the result contains RootOfs. Explicit solutions will be even longer, so I didn't try.

Edit: this becomes more manageable:

sol := solve(sys,vars): length(sol);
ssol := simplify(sol,size):  length(ssol);

To be precise, that warning is generated when parsing a procedure definition, not when executing (calling) it.
You can set

interface(warnlevel=1):

but that will suppress other warnings as well (see ?interface), so I would keep its default value 3.

First 19 20 21 22 23 24 25 Last Page 21 of 43