Graphical Programming with MapleSim in Vector Mechanics to Structures 2D

At the present time before constructing or starting to develop a mechanical structures project it is necessary to model it using graphic programming; In this opportunity and used MapleSim as a computational tool belonging to the company Maplesoft. The modern approach to modeling and simulation makes the fabrication of complex designs easy to solve. We will cover some examples taken from the engineering being implemented in Maplesim with insertion of physical objects; To be seen in real time through video output; Then integrates with Maple to analyze the equations and data through the static and dynamic behavior of the fabricated. Solved methods of physical block components include functionality for many domains: rotational and translational mechanics, multi-body dynamics, logic, and structural blocks; With techniques like: Drag-and-Drop Physical Modeling Environment and Create Custom Components Directly From Their Equations, thus the systems that would take hours or days to build from equations; In principle they can be created in a fraction of time using MapleSim, so it can incorporate significantly more complex graphical algorithms. In MapleSim, I use the revolutionary multibody technology that perfectly combines advanced multi-domain modeling tools to provide all the functionality you need in one environment.

FAST_UNT_2017.pdf

Lenin Araujo Castillo

Ambassador Maple - Perú

Let:

f:=x->1/sqrt(2*Pi)*exp(-x^2/2);

I.e. f is a standard Gaussian PDF.

Then (in Maple 2016.1):

Int(convert(f(x)*f(y)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity):
evalf(%);

Returns:

1.692568751

However (again in Maple 2016.1):

int(convert(f(x)*f(y)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity):
evalf(%);

Returns:

-0.5641895835

This is clearly incorrect, as the integral of a positive function must be positive.

This also seems to be a problem in which ever version of Maple is used behind the scenes on this forum.

gives:

It is suggested  

hypergeom([1/3, 2/3], [3/2], (27/4)*z^2*(1-z)) = 1/z

if z > 1. Here is my try to prove that with Maple:

Download simplification.mw

I want to get solutions of this system ,can anyone help me ?solutions.mw

Can someone help me with this:

reverse_eng.mw

happy new year!

Friends in Maple

I have a Vehicle Routing Problem I wish to cast as a integer model. VRP is a kind of TSP and knapsack problem hybrid  I would be grateful if someone can finish it. answers included. This isn't homework BTW.

VRP_IP.mw

I use maple 12 on my dual core laptop and i plan to buy maple 2016 and an 8 core system to get advandage of multiple processors. 

Does maple 2016 functions get advandage of multiprocessor systems or it will be the same as having one processor?

Here is one that the students generated which caused confusion. 

a := 0.76;
eq1 := 2*cot(a*sqrt(2*E)) = (2*E-5.4)/(sqrt(E*(5.4-E));
solve(eq1, E)

And the results are: 0., 4.411954070, 2.423743792

The problem is with the second answer because it does not exist. If we plot the LHS and the RHS of eq1 vs E where E=0..5.4

And it gets more interesting, if we calculate:

solve(evalf(eq1), E)

The answers are: 0., 2.423743793, 14.33807304+27.39159712*I

where the 3rd answer is again incorrect.

Finally, if a = 0.8 or larger, the incorrect answers disappear. 

Note - fsolve does handle this problem correctly. And despite my attempts to remind them to use fsolve, they see the solve command as the universal truth. Apparently this will be another teaching moment for next year.

So any thoughts about why this happens and why there is a difference in the outcomes between 0.76 and 0.8 for the value of a?

hi .how i can fsolve 8 equations, in which obtained from

for k from 2 to m do eq1[k_] end do

AND

for k from 2 to m do eq2[k_] end do  and so on........

And gain unknown functions as f3[2], f4[4],....

thanks?

fdm-maple.mw
 

Download fdm-maple.mw

hi.

how i can dsolve this differential equations?

thanks

ich.mw
 

Download ich.mw

In the test package below I want the the Value of Clr to be "B" when the package is loaded. How do I achieve that?
 

Download Module_Test_global_Var.mw

When I was using Maple 2015, I can start Maple in Linux Mint Sarah, if Maple 2015 was installed in another Linux Distribution (that is, like Debian Jessie), as follows.   Suppose I have installed Maple 2015 in /home/programfiles/maple/maple2015, which is the location in Debian Jessie.  Start Linux Mint Sarah, then mount  the /home partition of Debian on the mountpoint /debian-home.  Now copy the file /debian-home/programfiles/maple/maple2015/bin/xmaple to be xmaple-a, and change the path in xmaple-a. Then copy the file  /debian-home/programfiles/maple/maple2015/bin/maple to be maple-a, and make some changes in maple-a. Then open a terminal, run the command  ./debian-home/programfiles/maple/maple2015/bin/xmaple-a. Thus I can start Maple2015 in Linux Mint Sarah, althought I have reinstalled Maple 2015 in Linux Mint.

But the problem occurs, when I upgrade to Maple 2016, together with Linux Mint from Sarah to Serena, because the method above does not work.  It gives information saying that the Host ID is wrong. But when I start Debian Jessie, I can start Maple 2016, and it works very well. So, what should I do to start Maple2016 in Linux Mint Serena, without installing Maple2016 in Mint, because I have installed in Debian?  

Let L = {a1, a2, ..., an} (or L = [a1, a2, ..., an]) with is homogeneous polynomial. How to create function createDegree() return elements of degree <= n from L.

Expamle. L = {a-2b, b^2, (a+c)^2}.

createDegree(L,1) return {a - 2b}.

createDegree(L,2) return {b^2, (a - 2b)^2, (a+c)^2, a - 2b}.

createDegree(L,3) return {b^2, (a - 2b)^2, (a - 2b)^3, (a+c)^2, (a - 2b)b^2, (a - 2b)(a+c)^2, a - 2b}.

Thanks you very much.

As you see above, when I want to evaluate the definite integral, Maple can not provide the desired result. Why does this happen? How can I do so that Maple can calculate this integral?

I want to plot the following equation

((1/2)*x*tan(log((1/5)*(2*(x^2+y^2)))) = y

I tried:

with(plots, implicitplot)

implicitplot((1/2)*x*tan(log((1/5)*(2*(x^2+y^2)))) = y, x = -5 .. 5, y =-5..5, numpoints = 1000, gridrefine = 6, color = red)

but result is strange...