MaplePrimes Questions

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:

mm := [[a, 1], [a, b]];
mmseq := {seq(mm[ij],ij=1..2)};

how to filter above in seq when any items in the [] has 1

criteria is below

if mm[ij][1] = 1 or mm[ij][2] = 1 or mm[ij][3] = 1 then
Dummy := 0:
else
mm2 := [op(mm2), mm[ij]]:
end if

 

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

 

What is the best source of learning maple for an abecedarian to become a professional?

Hi

I can't understand difference between plots!

Please expailn it.

Thanks.

plot({sqrt(x+2*sqrt(x-1))+sqrt(x-2*sqrt(x-1)), sqrt(x-2*sqrt(x-1)), sqrt(x+2*sqrt(x-1))}, x = -3 .. 3)

 

 

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.

Hi,

I have a problem with the adaptive question designer: when I use the multiple choice question type then occasionally parts of the question environment appear multiple times in the text, duplicating each time I reopen the question to edit. This happens in particular if the answers are a bit longer (4-5 lines each). So far I couldn't figure out how to fix this, does anyone have a similar problem? 

Many thanks for your help!

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)

Maple 2015:

simplify(1-2*sin(x)^2);  gives 2*cos(x)^2-1

I looked at help trying to understand why Maple thinks 2*cos(x)^2-1 is simpler than 1-2*sin(x)^2 but did not see it. I was expecting to see cos(2*x) as a result.

Is there a place to understand more Maple's simplification rules other than the help page? http://www.maplesoft.com/support/help/maple/view.aspx?path=simplify%2fdetails

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.

First 1253 1254 1255 1256 1257 1258 1259 Last Page 1255 of 2434