MaplePrimes Questions

I'd want to create numerical vectors or matrices such that each entry follow the same normal law $N(a,b)$ (the entries are independent). How to do ?

Thanks in advance.

Dear Maple users

Sometimes I would like to place certain Maple content on my website, for example making it possible for the user to rotate 3d objects (created in Maple) directly on the webpage, in the browser. I did that on one of my own webpages, written in Danish though:

http://www.matematiksider.dk/soapbubbles.html

If you scroll to the Catenoid shape, almost at the bottom of the page, one should be able to click this object to have a Javaapplett started. I made this years ago with the help of the great tools created by Konrad Polthier:

http://www.javaview.de/maple/index.html

It wasn't easy, but I made it work. Things have changed since that time, however. Java is more and more "in bad standing". I was blocked when I tried opening it because of security. Banks is looking for alternatives, also because of security. And the mostly used browser now, Chrome, has stopped supporting it.

Now my question: Is there an alternative? Will HTML5 make it possible to make some kind of "Maple plugin" for Interactive Maple objects to be placed on webpages and interacted with in a Webbrowsers?

Erik

 

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

http://www.maplesoft.com/applications/view.aspx?SID=5083&view=html&L=F

Error, (in hanoi) cannot determine if this expression is true or false: 0 < n

N:=n:moves:=[];
hanoi:=proc(N,source,dest,aux) global moves;
if N>0 then
hanoi(N-1,source,aux,dest):
moves:=[moves[],[N,source,dest]]:
hanoi(N-1,aux,dest,source);
end if:
end proc:
sol:=hanoi(n,1,3,2);

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)

`&sigma;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&theta;o`*x+M1*(x+0)^3/factorial(3);

EIYo+`EI&theta;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&theta;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&theta;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&theta;ac` := diff(EIYac, x);

`EI&theta;o`+7500.000000*x^2

(13)

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

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

(14)

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

`EI&theta;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&theta;ac`, x);

15000.00000*x

(16)

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

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

(17)

Mef := diff(`EI&theta;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&theta;o` = 0

 

EIYo = 0

(22)

x := 4:

EIYo+4*`EI&theta;o`+160000.0000

(23)

x := 7.5:

EIYo+7.5*`EI&theta;o`+610931.2500+7.145833333*FCy

(24)

SOL := solve({CF1, CF2, CF3, CF4, SumFY, SumMA}, {EIyo, FAy, FCy, FEy, MA, `EIy&theta;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.

First 1248 1249 1250 1251 1252 1253 1254 Last Page 1250 of 2429