Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Greetings to all,

I want to make a "circular" list but I cannot find a way.

To be more precise, let K=[x1,x2,x3,x4], is any way to set K[0]:=K[4] and K[5]:=K[1]?

 

Thanx a lot,

Giorgos

There is million items in list, 

what is the difference between seq and for loop ?

is seq faster than for loop?

it is very slow when running code below, how to speed up this part of code?

 

https://gist.github.com/LovelyYanki/c0b61fbb9d5954b34e03#file-gistfile1-txt

HSKeyIn := Table();
for ij from 1 to 1685159 do
s := solve([mm[ij][1,1]=1,mm[ij][1,2]=1,mm[ij][1,3]=1],[a,b,c]):
if nops(s) > 0 then
if (rhs(s[1][1]) = 1 or rhs(s[1][1]) = 0) and (rhs(s[1][2]) = 1 or rhs(s[1][2]) = 0) and (rhs(s[1][3]) = 1 or rhs(s[1][3]) = 0) then
#print(lhs(indices(T3, pairs)[ij])):
#print("***"):
h := HilbertSeries([mm[ij][1,1],mm[ij][1,2],mm[ij][1,3]], {a,b,c}, z):
#print(h):
if not assigned(HSkeyIn[h]) then
if mod(ij, 100) = 0 then
print(ij):
end if:
HSkeyIn[h] := [[mm[ij][1,1],mm[ij][1,2],mm[ij][1,3]]]:
else
if mod(ij, 100) = 0 then
print(ij):
end if:
HSkeyIn[h] := [op(HSkeyIn[h]), [mm[ij][1,1],mm[ij][1,2],mm[ij][1,3]]]:
end if:

end if:
end if:
od:

Hello,

 

In a course I am looking at, the system of units is mostly in FPS.  Wich I never use.  But they ask sometime to convert it to SI.  But with Maple, I have a problem:

For Force, it use lb that I cannot convert into N

So what unit should I use to convert FPS to SI for the force.

 

I know that slug*ft*s^(-2) convert to N with no problem.

 

Thank you for your help.

 

 

--------------------------------------
Mario Lemelin
Maple 2015 Ubuntu 14.04 - 64 bits
Maple 2015 Win 7 - 64 bits messagerie : mario.lemelin@cgocable.ca téléphone :  (819) 376-0987

Cheers!

I'm having a problem with my student work, about to have a solution of 6 equations... Can help me in this file? i dont know how to solve this... this had-me a null solve...

 

 


Thanks for the help =)

restart

M1 := 0.15e5;

0.15e5

 

0.60e5

 

0

 

0.12e5

 

21000.00000

 

3

 

1

 

2.5

 

1

 

3

(1)

`σadm` := 175*10^6;

175000000

 

(1/300000)*L

 

210000000000

(2)

Atria := (3.5*12)/(LBC+LCD)

12.00000000

(3)

Ctria := LAB+LBC+(1/3)*(2*(LCD+LDE))

6.333333334

(4)

AiXil := Atria*Ctria

76.00000001

(5)

C := AiXil/Atria

6.333333334

(6)

``

``

``

SumFX := FAx;

FAx

(7)

SumFY := FAy+FCy+FEy-F5-QTria;

FAy+FCy+FEy-81000.00000

(8)

SumMA := FCy*(LAB+LBC)-F5*(LAB+LBC)+FEy*(LAB+LBC+LCD+LDE)+M1-MA-QTria*Ctria;

4*FCy-358000.0000+7.5*FEy-MA

(9)

NULL

``

``

EIYac := EIYo+`EIθo`*x+M1*(x+0)^3/factorial(3);

EIYo+`EIθo`*x+2500.000000*x^3

(10)

EIYce := EIYac+FCy*(x-4)^3/factorial(3)-F5*(x-4)^3/factorial(3)-q5*(x-4)^5/((3.5)*factorial(5));

EIYo+`EIθo`*x+2500.000000*x^3+(1/6)*FCy*(x-4)^3-10000.00000*(x-4)^3-28.57142857*(x-4)^5

(11)

EIYef := EIYce+FEy*(x-7.5)^3/factorial(3)+(1/3)*q5*(x-7.5)^5/factorial(5);

EIYo+`EIθo`*x+2500.000000*x^3+(1/6)*FCy*(x-4)^3-10000.00000*(x-4)^3-28.57142857*(x-4)^5+(1/6)*FEy*(x-7.5)^3+33.33333333*(x-7.5)^5

(12)

`EIθac` := diff(EIYac, x);

`EIθo`+7500.000000*x^2

(13)

`EIθce` := diff(EIYce, x);

`EIθo`+7500.000000*x^2+(1/2)*FCy*(x-4)^2-30000.00000*(x-4)^2-142.8571428*(x-4)^4

(14)

`EIθef` := diff(EIYef, x);

`EIθo`+7500.000000*x^2+(1/2)*FCy*(x-4)^2-30000.00000*(x-4)^2-142.8571428*(x-4)^4+(1/2)*FEy*(x-7.5)^2+166.6666666*(x-7.5)^4

(15)

``

Mac := diff(`EIθac`, x);

15000.00000*x

(16)

Mce := diff(`EIθce`, x);

-45000.00000*x+FCy*(x-4)+240000.0000-571.4285712*(x-4)^3

(17)

Mef := diff(`EIθef`, x);

-45000.00000*x+FCy*(x-4)+240000.0000-571.4285712*(x-4)^3+FEy*(x-7.5)+666.6666664*(x-7.5)^3

(18)

``

Vac := diff(Mac, x);

15000.00000

(19)

Vce := diff(Mce, x);

-45000.00000+FCy-1714.285714*(x-4)^2

(20)

Vef := diff(Mef, x);

-45000.00000+FCy-1714.285714*(x-4)^2+FEy+1999.999999*(x-7.5)^2

(21)

``

x := 0:
``

`EIθo` = 0

 

EIYo = 0

(22)

x := 4:

EIYo+4*`EIθo`+160000.0000

(23)

x := 7.5:

EIYo+7.5*`EIθo`+610931.2500+7.145833333*FCy

(24)

SOL := solve({CF1, CF2, CF3, CF4, SumFY, SumMA}, {EIyo, FAy, FCy, FEy, MA, `EIyθo`});

"SOL:="

(25)

``

NULL

``

 

Download Equacoes_universais_T12_-_4.mwEquacoes_universais_T12_-_4.mw

Dear Maple users

I am delighted that Maple has builtin commands to plot so many polyhedrons in 3D. Here I am talking about the polyhedraplot command in the plots package. I was however disappointed that the socalled Truncated Icosahedron is not supported (not present in the supported list ...). My first question is:

1. Why isn't it supported?

It seems more relevant than many of the other polyhedrons which are supported. It is a member of the Archimedian Solids (see https://en.wikipedia.org/wiki/Truncated_icosahedron). Besides it is the basic structure for soccer footballs. I found out that a TruncatedIcosahedron command is available in the geom3d package. This command is able to deliver data for the faces and more. With this command I succeded in writing a small program to actually plot this polyhedron in 3D:

-----------------------------------------------------

restart;
with(geom3d):
with(plots):
TruncatedIcosahedron(football,point(C,(0,0,0)),1):

PlotFootball:=proc(tr::float)
    
    local i::integer,
         FootballFace::Vector(datatype=float[8]),
         plotFace::Vector(datatype=float[8]);    
    
    for i from 1 to 32 do
        FootballFace[i]:=faces(football)[i];
        plotFace[i]:=polygonplot3d(Matrix(FootballFace[i]),axes=none,scaling=constrained,transparency=tr);
    end do;
    
    display(seq(plotFace[i],i=1..32)):
    
end proc:

PlotFootball(0.00);

-----------------------------------------------------

Since I am not really experienced in programming in Maple, here is my last question:

2. Can I simplify something in my code above?

 

Best wishes,

Erik

 

People im with error to show a plot, this is a wor for universty about mechanic materials, and i'm with this error... do not show the plot... before, i've a error because the variable x has values, and i've to unassign to have this.. without unassign, has error in domin of the plot... i dont know how to do this ! And ive to send the work today. 

Someone can help me about this?


Value of variables:
functions
plot algoritm...
plot algoritm

Thanks !

in page 137 of an introduction to groebner bases

how to eliminate the redundant solution (y^2-x*z, 0, -x^2+y*w)

from 3 of them?
eliminate({y,y^2-x*z,-z}, {x, y, z, w});
eliminate({-x,0,y}, {x, y, z});
eliminate({w,-x^2+y*w,-x}, {x, y, z, w});

ma := allstructs(Permutation([1, 1, 1, 2, 2, 2, 3, 3, 3]), size = 3);

above is fast

but below is very slow.
ma2 := allstructs(Permutation(ma), size = 3);

just for all combinations of matrix , replicateM in haskell is the fastest.

in maple, ma2 := allstructs(Permutation(ma), size = 3); is very slow

 

 

Let us define a piecewise-linear continuous function:

restart; VP := Vector[row](16, {(1) = 10, (2) = 177.9780267, (3) = 355.9560534, (4) = 533.9340801, (5) = 711.9121068, (6) = 889.8901335, (7) = 1067.868160, (8) = 1245.846187, (9) = 1423.824214, (10) = 1601.802240, (11) = 1779.780267, (12) = 1957.758294, (13) = 2135.736320, (14) = 2313.714347, (15) = 2491.692374, (16) = 2669.670400}); VE := Vector[row](16, {(1) = 5.444193931, (2) = .4793595141, (3) = .3166653569, (4) = .2522053489, (5) = .2123038784, (6) = .1822258228, (7) = .1544240625, (8) = .1277082078, (9) = .1055351619, (10) = 0.8639065510e-1, (11) = 0.6936612570e-1, (12) = 0.5388339810e-1, (13) = 0.3955702170e-1, (14) = 0.2612014630e-1, (15) = 0.1338216460e-1, (16) = 0.1203297900e-2}); for i to 15 do p[i] := VE[i+1] < x and x <= VE[i], (VP[i+1]-VP[i])*(x-VE[i])/(VE[i+1]-VE[i])+VP[i] end do; g := unapply(piecewise(seq(p[i], i = 1 .. 15)), x);

for i to 15 do print(fsolve(g(x) = VP[i])) end do;

Why doesn't the fsolve command work if i = 4, 7, 9, 11, 14? There are workarounds:

print(DirectSearch:-SolveEqutions(g(x) = VP[i]));

and/or

VP := convert(VP, rational); VE := convert(VE, rational); print(solve(g(x) = VP[i]));

 How to explain such behavior of the fsolve command? That was asked but not answered in http://forum.exponenta.ru/viewtopic.php?t=13524&sid=025a140e7e00b99803c86060a5c0c33c .

NULL

 

strange_behavior.mw

Edit. Replaced worksheet.

Please help me to solve this integration

restart; with(LinearAlgebra); int(exp(-(ln(z/(snr*B^2))+4*sigma^2)^2/(32*sigma^2))*eta^2*(y/z)^((1/2)*eta^2-1)/(z*sqrt(32*Pi*sigma^2)*(2*sqrt(y*z))*(2*A[o]^(eta^2))), z, z = y/A[o]^2 .. infinity);

 

 

 

restart; with(LinearAlgebra); int(exp(-(ln(z/(snr*B^2))+4*sigma^2)^2/(32*sigma^2))*eta^2*(y/z)^((1/2)*eta^2-1)/(z*sqrt(32*Pi*sigma^2)*(2*sqrt(y*z))*(2*A[o]^(eta^2))), z, z = y/A[o]^2 .. infinity)

I need to show that the following expression,
a^3b-a^3c+a^3z+a^3x+a^3y-a^2bx+a^2by+a^2cx-a^2cy-a^2zx+a^2zy-a^2x^2+a^2y^2-abcz-abcx-aczx-acx^2+b^2c^2+2bc^2x+c^2x^2-b^2c-2bcx-cx^2,

is positive

given that:

1. a,b,c,x,y,z are positive real numbers

2. a>b+x

3. c<b+y

I know a priori that the expression is indeed positive, but I do not know how to show it, or how to use Maple to do it?

Specifically, how can I use Maple to **partially factorize** the expression in terms of the expressions a-b-x and c-b-y?

Thanks for any help.

I have created and saved a MAPLE module in an .mla archive. The module contains three procedures A, B, C, where

A calls, B and C.  

Once the module library has been loaded, A acccepts inputs and generates outputs.

Is it possible to create a MAPLE player worksheet which calls the module and share it with a Maple Player (only) user, so that they can then supply the inputs and observe the outputs from A using the Maple Player programme components?

Can anyone help?

MRB

 

 

 

I have created and saved a MAPLE module in an .mla archive. The module contains three procedures A, B, C, where

A calls, B and C.  

Once the module library has been loaded, A acccepts inputs and generates outputs.

Is it possible to create a MAPLE player worksheet which calls the module and share it with a Maple Player (only) user, so that they can then supply the inputs and observe the outputs from A using the Maple Player programme components?

Can anyone help?

MRB

 

 

 

For the past thirty years, I have used several mathematical packages for problem solving and graphing. It all started with spreadsheet software that really helped speedup calculations compared to calculators. As many people do, once I had one tool I then started looking for another that would offer even more capabilities and features. I tested several of the very early math software but none really did all that I wanted until I came across Maple while I was working at SPAR Aerospace in Canada. For me, the rest is history. As long as I had a copy of Maple, it was all that I needed.

On occasions when I did not have a copy of this amazing software, I resorted to spreadsheets once more to complete fairly large and complex projects involving large databases and large numbers of calculations, especially when performing What-If scenarios. One distinct disadvantage of using a spreadsheet was the cryptic form of equation writing. I had to divide one long equation into several sections in different cells and then add them all up, which clearly is not good for documentation of the calculations. It is also very confusing for other engineers to know what that equation is or what it does. The development of the full engineering spreadsheet took months to complete, debug and verify. During this process, when I had errors, it was often very difficult to track exactly where the problem was, making the debugging process time consuming and sometimes very frustrating.

Having worked with Maple before, I remembered how easy it was to enter equations in a very familiar, readable math format. The real power of this software is that it allows you to write the equation(s) anyway you like and solve for any given parameter, unlike spreadsheets where you have to solve the problem first, by hand, for the parameter you want and then get the spreadsheet to calculate the value. I remember one time a few years ago when I wrote nine or ten simultaneous differential equations all in symbolic form and asked Maple to calculate certain parameters in a fully symbolic form. To my utmost disbelief, the answer came back within few minutes. With results in hand, I was able to quickly finish my research, and the results were published at PCIM Europe 2005 in “Distributed Gate ESR and its Effect on Shoot Through Performance at the Die Level”. I would never have gotten the results I needed if I was using a spreadsheet.

Even with much simpler systems of equations, finding solutions with a paper and pencil was never an easy task for me. It took a very long time, and even then there was no guarantee that I did not make copying errors, accidentally leave out a term, or make a calculation error. After I found the correct solution, I then had the problem of plotting the results, which I often needed in 3-D. Plotting allowed much deeper insights into the interdependency of all the parameters and made it easy for me to concentrate on the important ones without wasting any time. I was very happy when I could pass all these tasks onto Maple, which could do them much faster and more reliably then I ever could. Maple is a software that allows me to go beyond routine engineering calculations and gives me the tools to reach levels of insight and understanding that were completely out of reach of the average engineer until a few years ago.

For the record, I have no business affiliations with Maplesoft. I’m writing this article because Maple makes such a difference in my work that I feel it is important to share my experiences so other engineers can get the same benefits.

Hi everybody:

I'm going to learn programming with maple 18, are there any good and new pdf files for learn it?

with regards...

 

First 1221 1222 1223 1224 1225 1226 1227 Last Page 1223 of 2224