Maple 18 Questions and Posts

These are Posts and Questions associated with the product, Maple 18

When I use listcontplot, the tickmarks on the axis show the ordinal of the point plotted, so if I have 20x50 points, it shows (1 to 21)x(1 to 51).

Is there any way to reescalate the axis to show the actual units of the points?, so for example (0.3 to 0.7)x(2.2 to 3.7).

Thanks.

hi every one

 im new in maple, i have this problem , the output result (nod) dont disply after the for loop how can i disply it ?( sorry for this question, im new in this prog.)

the 2nd question is the (Typesetting:-delayDotProduct) how can i resolve this???

plz help 

 

for i from 2 to 3 do for j from 2 to 3 do Typesetting[delayDotProduct](rod, [i, j], true) := s*imp*l*i*f*y*subs*({cx^2 = 1-cy^2-cz^2}, Typesetting[delayDotProduct](rod, [i, j], true)) end do end do;
Typesetting[delayDotProduct](rod, [1, 1], true) := collect(rod[1, 1], cos(theta));

Typesetting:-delayDotProduct

Implementation of Maple apps for the creation of mathematical exercises in
engineering

In this research work has allowed to show the implementation of applications developed in the Maple software for the creation of mathematical exercises given the different levels of education whether basic or higher.
For the majority of teachers in this area, it seems very difficult to implement apps in Maple; that is why we show the creation of exercises easily and permanently. The purpose is to get teachers from our institutions to use applications ready to be evaluated in the classroom. The results of these apps (applications with components made in Maple) are supported on mobile devices such as tablets and / or laptops and taken to the cloud to be executed online from any computer. The generation of patterns is a very important alternative leaving aside random numbers, which would allow us to lose results
onscreen. With this; Our teachers in schools or universities would evaluate their students in parallel on the blackboard without losing the results of any student and thus achieve the competencies proposed in the learning sessions.
In these apps would be the algorithms for future research updates and integrated with systems in content management. Therefore what we show here is extremely important for the evaluation on the blackboard in bulk to students without losing any scientific criteria.

FAST_UNT_2018.mw

FAST_UNT_2018.pdf

Lenin Araujo Castillo

Ambasador of Maple

 

There seems to be a bug: in some of the cases, the intersection command works, in others it doesn't.

I'd appreciate any hints as to why this is the case.

 

intersection_of_two_lines.mw

Hello everyone, 

please need some brain power here to help me with this problem. I'm stacked here and I cannot find a solution. 

I would like to find the intersection between a plane and a quadric surfaces. Please notice that the two equations are written like f(x,y) and g(x,z) respectively. My first problem is that I cannot visualize the two surfaces correctly as implicitplot3d of the quadric surfaces displays nothing in return. 

I want then to find the centre of the ellipse and plot the point and the two surfaces all together in the same 3d plot. Somehow I canont make it. What is wrong with me? 

I am pasting here my MWE. 

Any kind of help would be really appreciated. Thank you so much. 

Andrea

MWE_maple_help.mw
 

restart:

with(plottools):
with(plots):

with(Student:-MultivariateCalculus):
with(Student:-Precalculus):
with(geometry):

with( Optimization ):

my_plane := 2024.30587449691-.341275505799078*x-3.89936179341114*y;
my_quadric := -10595.4104931095+6.73749956241827*x+42.1022818380012*z-0.654818649508000e-1*x^2-0.421174257681115e-1*z^2;

2024.30587449691-.341275505799078*x-3.89936179341114*y

 

-10595.4104931095+6.73749956241827*x+42.1022818380012*z-0.654818649508000e-1*x^2-0.421174257681115e-1*z^2

(1)

my_quadric_plot_3d := implicitplot3d(my_quadric, x=-10000..10000,
                                       y = -100..100,
                                      z=-10000..10000,color="DodgerBlue"):
my_plane_plot_3d := implicitplot3d(my_plane, x=-10000..10000,
                                       y = -100..100,
                                      z=-10000..10000,color="DarkOrange"):

my_ellipse := CompleteSquare( subs(y = solve(z = my_plane, y), my_quadric));

-0.4211742577e-1*(z-499.8202178)^2-0.6548186495e-1*(x-51.44553815)^2+99.68250

(2)

ellipse(e1, my_ellipse, [x,z]):
XY_centre_e1     := coordinates(center(o,e1)):
my_center_point := [XY_centre_e1[1],XY_centre_e1[2],
                       subs(x=XY_centre_e1[1], y=XY_centre_e1[2], my_plane)]:
my_center_point_plot := pointplot3d(my_center_point, color = "red"):

display(my_center_point_plot, my_plane_plot_3d, my_quadric_plot_3d);

 


 

Download MWE_maple_help.mw

 

 

Does any body have example of maple code for solving optimal control problem in deterministic model using Pontryagin's maximum (or minimum) principle?

Hi

I d like to limit my solution to a real (non complex) solution
There should be a simple solution to my calculation but maple can t process the solutions in one of my "solve" commands.

How can do I tell maple to limit itself to one
 

NULL

restart

with(Student[Calculus1]):


#geometry [mm]

b := 250:

h := 720:

ds := 70:

d := h-ds:

As := 3000:


#concrete [MPa]

fck := 30:

fcm := fck+8;

38

(1)

Ecm := 33000;

33000

(2)

`ϵc1` := 2.2*(1/1000);

0.2200000000e-2

(3)

eta := `ϵc`/`ϵc1`;

454.5454545*`ϵc`

(4)

Ec1 := fck/`ϵc1`;

13636.36364

(5)

k := 1.05*Ecm/Ec1;

2.540999999

(6)

sigma := fcm*(-eta^2+eta*k)/(1+(k-2)*eta);

38*(-206611.5702*`ϵc`^2+1154.999999*`ϵc`)/(1+245.9090904*`ϵc`)

(7)


#steel [MPa]

Es := 200000:

fsy := 400:

fsu := 600:

`ϵy` := fsy/Es;

1/500

(8)

`ϵsh` := 0.9e-2:

`ϵsu` := 0.75e-1:

P := 4:

`ϵs` := `ϵcm`*(d-c)/c;

`ϵcm`*(650-c)/c

(9)

i := 1;

1

(10)

for `ϵcm` from .1*(1/1000) by .1*(1/1000) to 10*(1/1000) do `ϵs` := `ϵcm`*(d-c)/c; T[1] := `ϵs`*Es*As; T[2] := fsy*As; T[3] := (fsu+(fsy-fsu)*((`ϵsu`-`ϵs`)/(`ϵsu`-`ϵsh`))^P)*As; C := b*c*(int(sigma, `ϵc` = 0 .. `ϵcm`))/`ϵcm`; eq[1] := T[1] = C; `ϵl`[1] := `ϵy`; eq[2] := T[2] = C; `ϵl`[2] := `ϵsh`; eq[3] := T[3] = C; cc := max(solve(eq[1], c)); `ϵss` := subs(c = cc, `ϵs`); Ta := subs(c = cc, T[1]); if `ϵss` >= `ϵl`[1] then cc := max(solve(eq[2], c)); Ta := subs(c = cc, T[2]); `ϵss` := subs(c = cc, `ϵs`) end if; if `ϵss` >= `ϵl`[2] then cc := max(`assuming`([solve(eq[3], c, useassumptions)], [c::real])); Ta := subs(c = cc, T[3]); `ϵss` := subs(c = cc, `ϵs`) end if; M[i] := b*cc^2*fcm*(int(sigma*`ϵc`, `ϵc` = 0 .. `ϵcm`))*10^(-6)/`ϵcm`^2+T*(d-cc)*10^(-6); phi[i] := `ϵcm`/cc; cd[i] := cc/d; print(`ϵcm`, `ϵss`, Ta/As); i := i+1 end do

0.1000000000e-3, 0.1955232439e-3, 39.10464877

 

0.2000000000e-3, 0.3845102290e-3, 76.90204580

 

0.3000000000e-3, 0.5671741821e-3, 113.4348364

 

0.4000000000e-3, 0.7437144096e-3, 148.7428819

 

0.5000000000e-3, 0.9143174400e-3, 182.8634880

 

0.6000000000e-3, 0.1079158043e-2, 215.8316087

 

0.7000000000e-3, 0.1238400148e-2, 247.6800297

 

0.8000000000e-3, 0.1392197667e-2, 278.4395334

 

0.9000000000e-3, 0.1540695238e-2, 308.1390476

 

0.1000000000e-2, 0.1684028897e-2, 336.8057793

 

0.1100000000e-2, 0.1822326682e-2, 364.4653363

 

0.1200000000e-2, 0.1955709188e-2, 391.1418377

 

0.1300000000e-2, 0.2226921078e-2, 400

 

0.1400000000e-2, 0.2583745185e-2, 400

 

0.1500000000e-2, 0.2954159196e-2, 400

 

0.1600000000e-2, 0.3336139462e-2, 400

 

0.1700000000e-2, 0.3727768551e-2, 400

 

0.1800000000e-2, 0.4127227927e-2, 400

 

0.1900000000e-2, 0.4532791280e-2, 400

 

0.2000000000e-2, 0.4942818378e-2, 400

 

0.2100000000e-2, 0.5355749454e-2, 400

 

0.2200000000e-2, 0.5770100024e-2, 400

 

0.2300000000e-2, 0.6184456101e-2, 400

 

0.2400000000e-2, 0.6597469817e-2, 400

 

0.2500000000e-2, 0.7007855358e-2, 400

 

0.2600000000e-2, 0.7414385220e-2, 400

 

0.2700000000e-2, 0.7815886744e-2, 400

 

0.2800000000e-2, 0.8211238846e-2, 400

 

0.2900000000e-2, 0.8599369121e-2, 400

 

0.3000000000e-2, 0.8979250971e-2, 400

 

Error, complex argument to max/min

 

`ϵcm` := 0.35e-2; `ϵs` := `ϵcm`*(d-c)/c; T[1] := `ϵs`*Es*As; T[2] := fsy*As; T[3] := (fsu+(fsy-fsu)*((`ϵsu`-`ϵs`)/(`ϵsu`-`ϵsh`))^P)*As; C := b*c*(int(sigma, `ϵc` = 0 .. `ϵcm`))/`ϵcm`; eq[1] := T[1] = C; `ϵl`[1] := `ϵy`; eq[2] := T[2] = C; `ϵl`[2] := `ϵsh`; eq[3] := T[3] = C; cc1 := max(solve(eq[1], c)); `ϵss1` := subs(c = cc, `ϵs`); Ta1 := subs(c = cc, T[1]); cc2 := max(solve(eq[2], c)); `ϵss1` := subs(c = cc, `ϵs`); Ta2 := subs(c = cc, T[2]); cc3 := max(solve(eq[3], c)); `ϵss3` := subs(c = cc, `ϵs`); Ta3 := subs(c = cc, T[3]); M[i] := b*cc^2*fcm*(int(sigma*`ϵc`, `ϵc` = 0 .. `ϵcm`))*10^(-6)/`ϵcm`^2+T*(d-cc)*10^(-6); phi[i] := `ϵcm`/cc; cd[i] := cc/d; print(`ϵcm`, `ϵss`, Ta/As); i := i+1

7501.663386*c

 

308.9753080

 

0.1055633997e-1

 

6333803.979

 

159.9645223

 

0.1055633997e-1

 

1200000

 

Error, complex argument to max/min

 

0.1055633997e-1

 

1254623.652

 

0.35e-2, 0.9349901112e-2, 400

(11)

``

``

``

``

``

NULL

NULL

``

``

``

``

``

``

``

``

``

``

``

``

``

NULL

``

``

``

``

``

``

``

``

``


 

Download mathias.mw

simple non complex solutions ?

Hi

I m trying to do a nested loop but for some reason maple tells me my for loop is unterminated 

I d love some feedback on my code and on how to make it work
 

NULL

restart

with(Student[Calculus1]):


#geometry [mm]

b := 250:

h := 720:

ds := 70:

d := h-ds:

As := 3000:


#concrete [MPa]

fck := 30:

fcm := fck+8;

38

(1)

Ecm := 33000;

33000

(2)

`ϵc1` := 2.2*(1/1000);

0.2200000000e-2

(3)

eta := `ϵc`/`ϵc1`;

454.5454545*`ϵc`

(4)

Ec1 := fck/`ϵc1`;

13636.36364

(5)

k := 1.05*Ecm/Ec1;

2.540999999

(6)

sigma := fcm*(-eta^2+eta*k)/(1+(k-2)*eta);

38*(-206611.5702*`ϵc`^2+1154.999999*`ϵc`)/(1+245.9090904*`ϵc`)

(7)


#steel [MPa]

Es := 200000:

fsy := 400:

fsu := 600:

`ϵy` := fsy/Es;

1/500

(8)

`ϵsh` := 0.9e-2:

`ϵsu` := 0.75e-1:

P := 4:

`ϵs` := `ϵcm`*(d-c)/c;

`ϵcm`*(650-c)/c

(9)

i := 1;

1

(10)

"for epsiloncm from 0.1/(1000) by 0.1/(1000)to (3.5)/(1000) do  print(epsiloncm);     epsilons:=(epsiloncm)/(c)*(d-c):  T[1]:=epsilons*Es*As:  T[2]:=fsy*As:  T[3]:=fsu+(fsy-fsu)*((epsilonsu-epsilons)/((epsilonsu-epsilonsh)))^(P):  T[4]:=0:  C:=(b*c)/(epsiloncm)*int(sigma,epsilonc=0..epsiloncm):  eq[1]:=T[1]=C:  epsilonl[1]:=epsilony:  eq[2]:=T[2]=C:  epsilonl[2]:=epsilonsh:  eq[3]:=T[3]=C:  a:=1          for j from 1 by 1 to 3 do          if a=1 then          cc:=max(solve(eq[j],c)):          epsilons:=subs(c=cc,epsilons):          if epsilons<=epsilonl[j] then T:=subs(c=cc,T[j]): a:=0:          end if:          end if:          end do:  M[i]:=(b*cc^(2)*fcm)/(epsiloncm^(2))*int(sigma*epsilonc,epsilonc=0..epsiloncm)*10^(-6)+T*(d-cc)*10^(-6):  phi[i]:=(epsiloncm)/(cc):  T[i]:=T:  cd[i]:=(cc)/(d):    print(M[i],epsiloncm,phi[i]);  i:=i+1:  end do:                                        "

Error, unterminated loop

"for epsiloncm from 0.1/1000 by 0.1/1000to 3.5/1000 do  print(epsiloncm);     epsilons:=epsiloncm/c*(d-c):  T[1]:=epsilons*Es*As:  T[2]:=fsy*As:  T[3]:=fsu+(fsy-fsu)*((epsilonsu-epsilons)/(epsilonsu-epsilonsh))^P:  T[4]:=0:  C:=(b*c)/epsiloncm*int(sigma,epsilonc=0..epsiloncm):  eq[1]:=T[1]=C:  epsilonl[1]:=epsilony:  eq[2]:=T[2]=C:  epsilonl[2]:=epsilonsh:  eq[3]:=T[3]=C:  a:=1         for j from 1 by 1 to 3 do          if a=1 then  cc:=max(solve(eq[j],c)):  epsilons:=subs(c=cc,epsilons):          if epsilons<=epsilonl[j] then T:=subs(c=cc,T[j]): a:=0:          end if:          end if:          end do:  M[i]:=(b*cc^2*fcm)/(epsiloncm^2)*int(sigma*epsilonc,epsilonc=0..epsiloncm)*10^(-6)+T*(d-cc)*10^(-6):  phi[i]:=epsiloncm/cc:  T[i]:=T:  cd[i]:=cc/d:  print(M[i],epsiloncm,phi[i]);  i:=i+1:  end do:                                        "

 

``

``

``

``

``

``

NULL

NULL

``

``

``

``

``

``

``

``

``

``

``

``

``

NULL

``

``

``

``

``

``

``

``

``


 

Download mathias.mw

 

Thanks in advance

Hello,  after almost a decade not using Maple, I am trying to get going to solve a set of differential equations and I have some problems that are probably pretty basic.

 

1.  I have the following simplification equations

which rightfully gives

and

which I hoped would give

but gives

"Error, invalid left hand side in assignment".  What am I doing wrong?

 

2.  The above simplification equations will be used to simplify this vector:

Once I get passed the "Invalid left hand side assignm,ent" issue, I would like to reorganize each line of the above vector and simplify using the small equations above.  What command do I need to use to do this?  The Phi & Beta variables are not necessarily next to one another so the equations need to be reorganized and the substitution for Omega & Omega_dot done.

Any help is appreciated! Thanks

Hi!

I'm doing an assignment where we are supposed to create a random walk with 3 options (go straight forward, turn left, turn right), and the first step is to be taken to [1,0]. I am supposed to plot the result afterwards, but when I plot it, it's empty...

(The first # is forward, the second # is to the right, and the third # is to the left)

for i to M0 do

X[i, 0] := 1; Y[i, 0] := 0;

for j from 2 to M1 do

r := R3();

if r = 1 then X[i, j] := 2*X[i-1, j]-X[i-2, j]; Y[i, j] := 2*Y[i-1, j]-Y[i-2, j];

elif r = 2 then X[i, j] := X[i-1, j]+Y[i-1, j]-Y[i-2, j]; Y[i, j] := Y[i-1, j]-X[i-1, j]+X[1-2, j];

else X[i, j] := X[i-1, j]-Y[i-1, j]+Y[i-2, j]; Y[i, j] := Y[i-1, j]+X[i-1, j]-X[i-2, j];

end if;

end do;

end do:

Then, for the plot:

plot([[[([X[i, j], Y[i, j]], $j = 0 .. M1)], $i = 1 .. M0)], scaling = CONSTRAINED);

This is what I got for the plot:

Where have I gone wrong? 

Thanks in advance!

Hello everyone,

I am experiencing a strange result. The Theta_double_dot got equated to zero in line 1.7. Please see the red arrow.

(I am using the Physics package)

How can I solve this issue?

Dear Maple useres,

I have a polynomial in the form f(s)/g(s) where f and g are polynomial of order 4 with coefficients in symbolic form. I want to convert it in a form of (s-a)*(s-b)*(s-c)*(s-d)/(s-a1)*(s-b1)*(s-c1)*(s-d1)

When i provide the numic value of  symbolic coefficients, Maple find it easy to factor it in the desired form. Otherwise i get a messy solution in the form of root of... expression. Even with assume command, I see no difference in the result. Is there a way to get the above simple form rather than root of expression?

thanks.

 

Hello,

I used the interface(typesetting = extended) command to change the notation used for time derivative of o variable.

But it doesn't seem to be working well when the variable is in the supercript form (check the last line).

How can I solve this issue?

Mapleprimesquestion.mw

 

What is the meaning of the output of the command Weierstrassform?

E.g. what does mean x0^3 - 7*x0 + 88 +y0^2 ?

Maybe y0^2 = x0^3 + 7*x0 -88?

First 31 32 33 34 35 36 37 Last Page 33 of 86