I do a calculation o a simple integral

      int(cosh(cos(2*t)), t = 0 .. 2*Pi)

and get the answer 0 (zero), which I know is incorrect. The cosh function of a real argument is always positive.

If I represent the cosh function as a sum of two exponentials

(1/2)*(int(cosh(cos(2*t)), t = 0 .. 2*Pi) + int(cosh(-cos(2*t)), t = 0 .. 2*Pi))

I get    

2*Pi*BesselI(0, 1)

which looks much better.

What's the matter? Why Maple yields two different results? Victor.

I'm trying to create an iterator function

Iter := proc(ff, n)    
    local i,f;
    if n = 0 then return (x->x); end;
    f := x->ff:
    f := apply(unapply(ff))(x);
    for i from 1 to n-1 do
        f := ff(f);
    end;    
    return unapply(f,x);
end:

This seems to work except that it converts the independent variable to x, which can cause problems.

That is, if I pass something like n->sin(n) then it will output something like x->sin(sin(x))...

Of course, I do explicitly use x but it is because I don't know how to get the input's independent variable.

Ideally I'd like to be able to work on multiple arguments and iterate over the first, by default.

e.g., f#2(x,y,z) = f(f(x,y,z),y,z)

Ultimately the proc above doesn't work well as I want to be able to use it in all contexts(I could pass the value to Iter but I'd rather use function notation.

This must be so complicated process to describe the exact mechanism in human or animals ... and there are many reserchs on web for this subject ! Not expected Large answer !

By this way i am just curious to know by support of signal processing methods available in maple we can model a very basic vibrator as larynx and variable duct (what kind of filter we can call ?) or cavity that play the role of tongue or lips for changing of main stream to made various phonemes or frequencies

Here is a movie that show the section of human vocal system :

Anatomy

is it possible to get all the iteration values for a minimization problem

Hey,

I want to animate action of linear transformation on multiple polygons, in the end, it should be a grid of multiple parallel lines that will get morphed. The transformation is defined by a matrix A.

So far I got this code to work:

with(plots):
with(plottools):
setoptions(style=patch,view=[-5..5,-5..5], scaling=constrained):
basis_i:=polygon([[0,0],[1,0]],color=red):

N:=20:
F:=(x,y) -> (1-k/N)*<x,y> + k/N*(A.<x,y>):
L := transform((x,y)->convert(F(x,y),list)):
A:=<<2,1>|<-1,1>>;
                          

frames:=seq(L(basis_i),k=0..N):

display(frames, insequence=true);
 

It animates the basis vector nicely, but I can't seem to figure out, how to animate multiple lines, multiple polygons.

Anyone knows how to do this?

Any help will be much appreciated.

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.

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 ?