Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I am trying to expand the ridgid beam model to flexible beam model to 1.) use fixed walled end to calulate cantilever deflection from a constant load, and visualize deflection; and then 2.) incorporate a compliant walled joint (rotary joint) with stiff rotational spring. To get angles of rotation; and finally simulate a periodic force and calc/plot all reactions/theta(t).

 

There is a video on line of a round beam with fixed-fixed ends and a translating weight with deflection visulized, but only a you tube vid can be found of this set-up.   It is close to what I need.

 

Lastly, if anyone knows how to send the diagram to maple as in A:=MapleSim:-___()

This method to transfer from Sim to maple with diagram and current session(#?) was done in a tutorial vid from MapleSim session. I am using Maple/MapleSim 2017 which is difficult to follow changes in workspace and commands from vids/tutorials.

PS. I have connected through opening a new worksheet and explicitly opening by filename. 

ENTIRE_FUNCTION_CHECK_GOOD_COPY.mwSorry for the silly question, but i have actually tried to find it in the help pages i swear.

 

how do i change the font size of the output?

 

I am making some educational worksheets and the problem ive faced is that all of people i know
 


 

 

online are mathematica users or others, so i am designing the worksheets in such a way that it requires 0 understanding of maple for the student, but can  be readily used as an educational tool. 

 

They are all going to be aimed  for undergraduate level, and designed to investigate questions that are not answered by plugging the function into wolfram alpha. 

A draft of one attached, any general advice for this project appreciated.

 


I am trying to find for the equilibrium but why my solutions lost?

restart

interface(imaginaryunit = j);

I

(1)

unprotect(Pi)

lambda := k*tau*(C*Upsilon+I)/N;

k*tau*(C*Upsilon+I)/N

(2)

eqn1 := (1-p)*Pi+phi*V+delta*R-(mu+lambda+vartheta)*S;

(1-p)*Pi+phi*V+delta*R-(mu+k*tau*(C*Upsilon+I)/N+vartheta)*S

(3)

eqn2 := p*Pi+vartheta*S-(epsilon*lambda+mu+phi)*V;

p*Pi+vartheta*S-(epsilon*k*tau*(C*Upsilon+I)/N+mu+phi)*V

(4)

eqn3 := rho*lambda*S+rho*epsilon*lambda*V+(1-q)*eta*I-(mu+beta+chi)*C;

rho*k*tau*(C*Upsilon+I)*S/N+rho*epsilon*k*tau*(C*Upsilon+I)*V/N+(1-q)*eta*I-(mu+beta+chi)*C

(5)

eqn4 := (1-rho)*lambda*S+(1-rho)*epsilon*lambda*V+chi*C-(mu+alpha+eta)*I;

(1-rho)*k*tau*(C*Upsilon+I)*S/N+(1-rho)*epsilon*k*tau*(C*Upsilon+I)*V/N+chi*C-(mu+alpha+eta)*I

(6)

eqn5 := beta*C+q*eta*I-(mu+delta)*R;

beta*C+q*eta*I-(mu+delta)*R

(7)

mu := 0.1e-1;

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-2

 

.2

 

11610

(8)

Equilibria := solve({eqn1 = 0, eqn2 = 0, eqn3 = 0, eqn4 = 0, eqn5 = 0}, {C, I, R, S, V});

Warning, solutions may have been lost

 

{C = 0., I = 0., R = 0., S = 5946.585366, V = 5663.414634}

(9)

``


 

Download Equilibria.mw
 

restart

interface(imaginaryunit = j);

I

(1)

unprotect(Pi)

lambda := k*tau*(C*Upsilon+I)/N;

k*tau*(C*Upsilon+I)/N

(2)

eqn1 := (1-p)*Pi+phi*V+delta*R-(mu+lambda+vartheta)*S;

(1-p)*Pi+phi*V+delta*R-(mu+k*tau*(C*Upsilon+I)/N+vartheta)*S

(3)

eqn2 := p*Pi+vartheta*S-(epsilon*lambda+mu+phi)*V;

p*Pi+vartheta*S-(epsilon*k*tau*(C*Upsilon+I)/N+mu+phi)*V

(4)

eqn3 := rho*lambda*S+rho*epsilon*lambda*V+(1-q)*eta*I-(mu+beta+chi)*C;

rho*k*tau*(C*Upsilon+I)*S/N+rho*epsilon*k*tau*(C*Upsilon+I)*V/N+(1-q)*eta*I-(mu+beta+chi)*C

(5)

eqn4 := (1-rho)*lambda*S+(1-rho)*epsilon*lambda*V+chi*C-(mu+alpha+eta)*I;

(1-rho)*k*tau*(C*Upsilon+I)*S/N+(1-rho)*epsilon*k*tau*(C*Upsilon+I)*V/N+chi*C-(mu+alpha+eta)*I

(6)

eqn5 := beta*C+q*eta*I-(mu+delta)*R;

beta*C+q*eta*I-(mu+delta)*R

(7)

mu := 0.1e-1;

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-2

 

.2

 

11610

(8)

Equilibria := solve({eqn1 = 0, eqn2 = 0, eqn3 = 0, eqn4 = 0, eqn5 = 0}, {C, I, R, S, V});

Warning, solutions may have been lost

 

{C = 0., I = 0., R = 0., S = 5946.585366, V = 5663.414634}

(9)

``


 

Download Equilibria.mw

 

It is a truth universally acknowledged, that a single man in possession of an algorithm for calculating a Groebner basis, must be in want of an algorthim for calculating a reduced Groebner basis.

It seems odd that i can't find something in the Groebner package - if there isn't something there, I assume that there is a well known piece of code for doing this!

I want to determine components of a vector in explicit form.  Let I have vector nn[mu] with components (1,0,0,0). Is it possible to explain to Maple it?

I tried to do the next:

with(Physics);
Setup(dimension = 4);
Setup(metric = {(1, 1) = -1, (2, 2) = -1, (3, 3) = -1, (4, 4) = 1});
Setup(mathematicalnotation = true);
Coordinates(X);
Setup(tensors = {nn[mu](X)});
nn[mu] = Matrix(1, 4, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 0, (1, 4) = 1});
nn[mu]*k[mu];
SumOverRepeatedIndices(%)

 

The answer that maple gives is:

nn[1]*k[`~1`]+nn[2]*k[`~2`]+nn[3]*k[`~3`]+nn[4]*k[`~4`]

 

So, doesn't Maple understand that nn[1]=nn[2]=nn[3]=0?

I am still working on the worksheet linked below.  But I have run into trouble solving the equation symbolicly within for the variable t.  I think this is due to the fact there are multiple solutions and the commands I employed through the GUI interface is not capable of handling this issue?  For example, the solutions for sin(pi*t/T) would be N*T.  MAPLE is simply stating t=0.  So I think this is why my solutions are failing to produce results.

What other commands should I be employing?

solving_transcendental.mw

Imagine three isosceles triangles with coordinates of each stored in a matrix such as:

> coord1:=Matrix(3,2,[0,0,5,0,2.5,4]);

> coord2:=Matrix(3,2,[2,0,7,0,4.5,4]);

> coord3:=Matrix(3,2,[4,0,9,0,6.5,4]);

and plotted together as follows (only the first is shown)

> PLOT(CURVES([[0,0],[5,0],[2.5,4],[0,0]]),COLOR(HUE,1));

I want to create an animation with stationary (original) isosceles triangles in the background along with new isosceles triangles generated by incrementally (say 4 increments) moving only the apex of each triangle until it touches the base of the triangle; obtained by multiplying the apex’s vertical coordinate by (say ¼) in each increment using a do loop. 

Your help is greatly appreciated.

Thank you in advance.

  

 

Dear sir,

 

I request the Maple experts to help me on the above cited subject.

 

 

With best regards.

 

Mr M ANAND,

Associate Professor in Mathematics

Hi everyone,

 

I have below integral. The result must be real. I would like to simplify and get to real part but I couldnt.

> fs := (1-1/sqrt((Uo+U*sin(x))^2+a^2))*(Uo+U*sin(x)); fss := `assuming`([2*(int(fs*sin(x), x = 0 .. Pi))/Pi], [U::real, U::positive, a::real, a::positive, Uo::real, Uo::positive]);

 

I used evalc but it didn't give good results as well.

 

thanks

Maple 2015

Using with(combinat) the permutation of {a,b,c} is determined.

>restart:
>with(combinat):
>permute({a, b, c})
                  [[a, b, c], [a, c, b], [b, a, c], [b, c, a], [c, a, b], [c, b, a]]

The tree diagram of this permutation is

    

In Maple, using with(combinat) and with(GraphTheory), when I attempt to draw the permutation I get the following error:

>L := permute({a, b, c});
       L := [[a, b, c], [a, c, b], [b, a, c], [b, c, a], [c, a, b], [c, b, a]]
>DrawGraph(L);
  Error, invalid input: GraphTheory:-DrawGraph expects its 1st argument, H, to be of type       {GRAPHLN, list(GRAPHLN), set(GRAPHLN)}, but received [[a, b, c], [a, c, b], [b, a, c], [b, c,      a], [c,   a, b], [c, b, a]]

On Maple, again using with(combinat) and with(GraphTheory) the command permute(3) is used.  The results are manually configured as node-connection lines.  A fair representation of the tree diagram is configured by Maple, although the diagram has numeric instead of alpha configurations, and the a,b,c structure shown above is not easily recognized.

Any suggestions on developing a procedure that will graph (draw) an alpha-labeled permutation welcomed.  Thanks!  WC44_Permutation_Graph.mw

Am trying to teach myself "circle inversion" to solve circle geometry problems.  Towards this end I wrote a Maple 7 program to draw a Pappus chain of circles.  Pappus lived around 200AD, and I read that he considered this problem ancient!  

   The output diagram may not come out on this forum, but it comes out as a collection of tangential circles, and I have no problems with that.  In drawing the circles I was thinking of using the infinity symbol.  The following line uses the infinity, but with a minus sign.  Maple is "happy" with this and outputs square braces, suggesting an empty list?  I changed the program by omitting the -ve sign - and Maple came up with an error message.  I would have expected the output to include the value zero  ie [0]  I think it is just that the infinity symbol is not allowed in the seq command. 

[seq( i, i=0..-infinity )];  #???  why the -?,

 

In the next part I define X, Y and rad as the coordinates and radii of a sequence of circles: subscript i, ranging from start to nterm, (for the number of terms).  These are set at 1 and 286.   The program works OK with these values, but if the variable start is changed to zero, 0, it doesn't work.  I thought I might have a division by zero, but can't find one.  I also drew a circle with a radius of zero - no problem!  (...except it's difficult to see:-))   

   My main puzzle is why does it come up with the error message "invalid subscript selector" when the variable start is set to 0?

With i=start=0 I am asking it to draw a circle of radius zero at some specified coords - but it doesn't like it!  

  Any comments would be appreciated.  Merry Christmas.

David

 

 

restart:

# # # # # # # # # # # # # # # # # # # # # # # # # # # #

# Pappus Chain of circles

# # # # # # # # # # # # # # # # # # # # # # # # # # # #

with(plots):

with(plottools):

#Maple 7

#r is ratio of AC/AB = diam of large circle, to a smaller internally tangential one

r:=1.35:

#Why will start=0 not work?

start:=1:

#nterm=286 is smallest value where axes have "true" origin of (0.0)

#nterm<= 285 puts "x-axis" below y=-0.6

nterm:=286:

X := [seq( 0.5*r*(1+r)/(((i*(1-r))^2)+r), i=start..nterm )]:

 

Y := [seq( i*r*(1-r)/(((i*(1-r))^2)+r), i=start..nterm )]:

rad:=[seq(0.5*r*(1-r)/((i*(1-r))^2+r), i=start..nterm)]:

[seq( i, i=0..-infinity )];  #???  why the -?

#Gives [] output with -ve sign:  empty list?  Why is 0 not included?

# Not allowed to use with seq

#ccc:=[seq(circle([X[i],Y[i]],rad[i]), i=1..5)]:  #nops(X))]:

ccc:=[seq(circle([X[i],Y[i]],rad[i]), i=start..nterm)]:

#c := circle([1,1], 0, color=red):

#plots[display](c, scaling=constrained);

plots[display](ccc, scaling=constrained);

 

OUTPUT

 

I am working on a simple PDE and am trying to find the values of the constant coefficients in the PDE such that a function of a particular form satisfies the equation.  The function is a sum of general functions of x, y and z in the argument: for example, one of them is $f(x-2y+z)$.  I then have a PDE with constant coefficients which contains a sum of some mixed and some single second partial derivatives of that function with constant coefficients in front of some of the derivatives.  Is there any relatively easy package with Maple which would enable me to find the values of those coefficients such that the function of the form I have will satisfy the PDE?  I can provide more exact details of the PDE if necessary but didn't want to get bogged down in details as interested in the principle.

I am attempting to use Maple to solve a matrix equation of the form aB^2+bB+cI=B^(−2), where B is a 3×3 matrix whose entries are constants which have been specified, and I is the 3×3 identity matrix (i.e. to find the values a, b and c which satisfy the equation).

My idea was to use a piece of code like

A:=Matrix([[7,4,-2],[4,7,5],[2,-3,8]]);
d:=Vector([[8],[5],[2]]); 

for matrix inversion to solve a system of three linear equations but that would assume that the right-hand side of the equation is a column vector when it is actually another 3×33×3 matrix. One other approach I took was to equate the left and right hand sides so that you add the three matrices on the LHS and then equate the resulting each entry with the corresponding entry of the matrix on the right hand side, however this gives me a set of 9 linear equations for 3 unknowns.  When I try to solve these with fsolve it says that there are more equations than there are unknowns. 

Hey

Can someone figure out why the output does not produce hyperlinked:

A1

B1

A2

B2 

(in the same order as Games).

hyperlink_do-loop.mw

Hi,

I need help. My Maple 2017 do not launch in my MacBook Air (11-inch, Mid 2011, 1.8 GHz Intel Core i7, 4 GB 1333 MHz DDR3, Intel HD Graphics 3000 384 MB) . Whe I open it, I only see the icon of Maple 2017 on the top left on my screen.

 

Now I have macOs High Sierra but with the older operating system it did not launch either.

 

Can you help me, please?

 

Best Regards,

Clara

First 867 868 869 870 871 872 873 Last Page 869 of 2215