MaplePrimes Questions

I try to fill an Array with a procedure to have pointpairs for a plot. Here is my problem. It doesn't work I wont it to.

 


 


 

Can someone help my please? The first procedure works, I tested it with some numbers, but the second one is the Problem.

Hi dear maple team. i have a question on integration and i need a "real" and "finite" solution with any assumption or options. thanks for the help.


 

restart

f := ((1 - a)^2 + a^2*((1 - exp(-y))*(1 - exp(-x)) - 2 + exp(-x) + exp(-y)) + a*(2 - exp(-x) - exp(-y) + (1 - exp(-y))*(1 - exp(-x))))/(1 - a*exp(-x)*exp(-y))^3;

((1-a)^2+a^2*((1-exp(-y))*(1-exp(-x))-2+exp(-x)+exp(-y))+a*(2-exp(-x)-exp(-y)+(1-exp(-y))*(1-exp(-x))))/(1-a*exp(-x)*exp(-y))^3

(1)

a := 0.3;f

.3

 

(.91+.39*(1-exp(-y))*(1-exp(-x))-.21*exp(-x)-.21*exp(-y))/(1-.3*exp(-x)*exp(-y))^3

(2)

s := 2*evalf(int((int(f*exp(-x)*exp(-y), x = 0 .. y + t,AllSolutions)), y = 0 .. infinity,AllSolutions)) assuming real ;

 

 


 

Download stat1.mw

Hello, Probably there is a way to do this easily but I do not quickly find it within the help.

I want

rand(0..1)

to give a true! random number and not always the same number; otherwise it should be called

predefinedlist()

Seed is deprecated, not sure it would help though. So how do I go abouts?

In this code created by Acer a long while back. I can't figure out how to increase the length of the text rotating so the letters don't collide with one another.

Here is the code.

restart:

#
# One (of many) STL format collections of the English alphabet:
#     http://www.stlfinder.com/model/all-alphabet-letters-a-z
#
basedir:=cat(kernelopts(homedir), "/My Documents/stl/All_Alphabet_Letters_A-Z"):

getletter:=L->Import(cat(basedir, "/Letter_", L, ".stl"),
                     orientation=[-90,0,0],title=""):

letters := [m,a,p,l,e]:

for letter in letters do
  if not assigned(AZ[letter]) then
    AZ[letter]:=getletter(letter);
  end if;
end do:

with(plottools): with(plots):

Rplot:=display(AZ[letters[1]],
               seq(translate(AZ[letters[i]], 40*(i-1)+15, 0, 0),
                   i=2..nops(letters))):

Rplot:=transform((x,y,z)->[z,x,y])(scale(Rplot, 5/6*2*Pi/200, 1/5, 1/20, [0,0,10])):

display(changecoords(Rplot,cylindrical), scaling=constrained, axes=none,
        orientation=[0,90,0], viewpoint=[circleright], frames=50);

Here is the file with the letters

All_Alphabet_Letters_A-Z.zip

Good day everyone and happy new year in advance. I'm trying to write a code for a series but its not substituting. Below is the attached file. Thanking you very much in anticipation for your help.

 

Testing.mw

I went to save a file in Maple and it had large amounts of calculated data. 

Maple didn't ask me if I also wanted to save that data like it usually does - it just saved it.  Saved a 13Mb worksheet that took 5 minutes to relinquish control on my laptop. 

Hello,

is it possible to send compiled Maple functions (C function) to the nodes with the Grid toolbox? I can run the non-compiled function on the nodes but when I try to send a compiled version of the same function I get segmentation error. I have run the compiled function just fine in other nonparallel codes, and also with the Threads package. So I was wondering if there is something else you need to do to get it to work? Maybe compile on each node?

 

Dear Maple users

I just created a sunflower in Maple using the Golden angle. See attached file. I am using a scatterplot with solidcircles as symbols. But what I would like now is to plot integers instead of the solidcircles in order to investigate the order of these different solidcircles. To make it clear: I have 500 solidcircles. They have been plottet using two vectors containing the x and y coordinates of the individual solidcircles. The n'th solidcircle I want to replace with the integer n. How can I create that different plot?

NB! The plot might need to be larger or the integers smaller, but that is fine. 

 

Regards

Erik V

sunflower.mw

A flexible rope of constant linear density hangs from one fixed end.

The lower end is pulled aside and then released.

The rope swings back and forth with a whip-like action.

What equation of motion can be used to animate this action in Maple?

I cannot find any references through Google to the appropriate math. 

Is there any alternative to Mathematica's "RegionFunction" in Maple, with applications of controlling the plotting domain for surface plots ? This question has been asked previosly, and given a workaround by "Markiyan Hirnyk", suggesting to leave the function that is to be plotted undefined outside of the plot region. Something like:

f1 := proc(x,y) if x^2 + y^2 > 0 and x^2 + y^2 < 4 then 8 - x^2 - y^2 end if end proc:
f2 := proc(x,y) if x^2 + y^2 > 0 and x^2 + y^2 < 4 then x^2 + y^2 end if end proc:
plot3d({f1,f2}, -2..2, -2..2) 

The result is clearly not good looking. You can increase the number of points, but we need a better way to do this.

BTW, the similar question was asked in 2010. Are there any better alternatives, ~ 10 years later ?

Please explain the logic used to construct the plot below.

plots:-implicitplot(r >= cos(theta), r = 0 .. 1/3, theta = 0 .. 2*Pi, filledregions, coords = polar, numpoints = 5000, scaling = constrained)

I found that while using maple gui cases of corruption, stuck while typing etc. Hope this problem will get solved in coming version

Hi 

I have a differential equation. I can solve it using dsolve and plot it with odeplot. For example I obtain P(x).

I also have another equation. For example V in terms of x. So, I can plot it with plot(V).

Now I want to plot V in terms of P. How can I do that?

Thanks for your help

Dear Users!

Hoped everything going fine with you. I want to make animation of ten solutions as given bellow but fail to do that. Please see it fix the problem. I shall be very thankful to u.
SOLNSuy[1, 1] := 2.5872902469406659197*10^(-20)-.65694549571241255901*y+1.9708364871372376767*y^2-1.3138909914248251176*y^3-1.6010739356637904911*10^(-19)*y^4;
SOLNSuy[2, 1] := -4.002204462000*10^(-20)-1.7879176897079605225*y+5.3637530691192141414*y^2-3.5758353794044226250*y^3-6.8309939211286845440*10^(-12)*y^4;
SOLNSuy[3, 1] := -1.1953264450000*10^(-19)-3.2481690589079594122*y+9.7445071767154794599*y^2-6.4963381177952273213*y^3-1.2292726248071398400*10^(-11)*y^4;
SOLNSuy[4, 1] := -2.6720465500000*10^(-19)-4.9239979672954025921*y+14.771993901873204315*y^2-9.8479959345587718955*y^3-1.9029826928878336000*10^(-11)*y^4;
SOLNSuy[5, 1] := 3.416928541000*10^(-20)-6.7268498492441931137*y+20.180549547714413714*y^2-13.453699698443639810*y^3-2.6580790570532587008*10^(-11)*y^4;
SOLNSuy[6, 1] := -2.554122292000*10^(-20)-8.5884528335125514887*y+25.765358500514014457*y^2-17.176905666966875698*y^3-3.4587270427710613504*10^(-11)*y^4;
SOLNSuy[7, 1] := -9.206107680000*10^(-20)-10.456823708331499352*y+31.370471124965259849*y^2-20.913647416590986491*y^3-4.2774005353527132160*10^(-11)*y^4;
SOLNSuy[8, 1] := 1.9644186790000*10^(-19)-12.293003938471349390*y+36.879011815379230436*y^2-24.586007876856948223*y^3-5.0932823222176363520*10^(-11)*y^4;
SOLNSuy[9, 1] := -3.775112769000*10^(-19)-14.068404975282556550*y+42.205214925807397100*y^2-28.136809950465931724*y^3-5.8908824448577377280*10^(-11)*y^4;
SOLNSuy[10, 1] := 1.146281780000*10^(-19)-15.762658869974768890*y+47.287976609878780960*y^2-31.525317739837422477*y^3-6.6589592851037286400*10^(-11)*y^4;
plots[animate](plot, [SOLNSuy[A, 1], y = 0 .. 1], A = 1 .. 10);

Special request:
@acer @Carl Love @Kitonum @Preben Alsholm

Simple example to illustrate the desired functionality:
Say we have a 2D vector function which describes the position of a particle

r := t -> <5*cos(Pi*t), 5*sin(Pi*t)>

We want to define the velocity and acceleration as functions, so we could do something like

v := t -> <diff(r(t)[1], t), diff(r(t)[2], t)>

The problem now is that we cannot call our velocity function with numeric arguments.
A simple solution is to call the function via "subs", as in

subs(t = 2, v(t))

but IMHO, this is not very elegant and I guess inefficient. Is there a command that enables for pulling out the evaluated result from diff such that it can be used directly as a functional expression ? I.e., I want to be able to call

v(2)

directly, without having to do substitutions.

 

EDIT:

I found that you can do

v := <diff(r(t)[1], t, diff(r(t)[2], t)>
v := unapply(v, t)

but please provide your recommendations. Thanks

First 603 604 605 606 607 608 609 Last Page 605 of 2428