I use transition matrix guess the next note of Happy birthday, but it tends to the same note in the final few notes

what is wrong?

why?

restart;
with(LinearAlgebra):
with(AudioTools):

MK := proc(f::numeric)
AudioTools:-Create(duration=0.5,rate=44100,x->evalhf(sin(x*2*Pi*f/44100))):
end proc:

C4 := 261.62:
D4 := 293.67:
E4 := 329.63:
F4 := 349.23:
G4 := 392.00:
A4 := 440.00:
B4 := 493.88:
C5 := 523.25:

initstate := Matrix([[0,0,0,0,0,0,0],[0,0,0,0,0,0,0],[0,0,3,2,0,1,1],[0,0,2,0,0,0,0],[0,0,1,0,0,0,0],[0,0,1,0,1,0,0],[0,0,0,0,0,1,0]]);
statetemp := []:
for i from 1 to op(1,initstate)[1] do
if add(x, x in initstate[i]) = 0 then
statetemp := [op(statetemp), [0,0,0,0,0,0,0]]
else
statetemp := [op(statetemp), initstate[i]/add(x, x in initstate[i])]
end if:
od:
state := Matrix(7,7);
for i from 1 to op(1,initstate)[1] do
for j from 1 to op(1,initstate)[2] do
state(i,j) := statetemp[i][j];
od;
od;
state := MatrixMatrixMultiply(state, state);
member(max(state[3]),state,'pos');
member(max(state),state,'pos');
pos;
music := [A4,B4,C4,D4,E4,F4,G4];
song := [C4,C4,D4,C4,F4,E4,C4,C4,D4,C4,G4,F4,C4,C4]:
song := [op(song), music[pos[1]]];
for i from 1 to 10 do
state := MatrixMatrixMultiply(state, state);
member(max(state[pos[2]]),state,'pos');
member(max(state),state,'pos');
song := [op(song), music[pos[1]]];
od:
HappyBirthday:=Create(Vector(map(MK,song))):
Write("C:\\Users\\secret\\Documents\\HappyBirthday.wav",HappyBirthday);
 

Hello

I know that anames( 'user' ) recovers the contents of the workspace but I wonder whether it would be possible to recover only the contents of each m-file read within a worksheet.  

Many thanks

Ed

Dear All,

I have the following differential equation:

(diff(lambda(t), t) = lambda(t)*(rho - A*beta*h^(1 - beta + a)/(A*h^(1 - beta + a)*t*(1 - beta) + k^(1 - beta))))     (1)

where rho, A, beta, a, h and k, are known parameters. Note, that all of them are positive and parameters rho, beta and a are less than one.

When solving this equation by dsolve, I get:

lambda(t) = _C1*exp(rho*t)*(A*h^(1 - beta + a)*t*beta - A*h^(1 - beta + a)*t - k^(1 - beta))^(beta/(-1 + beta))    (2)

which when substituting, for instance, beta=0.3 into the obtained solution, transforms into:

lambda(t) = _C1*exp(rho*t)/(-0.7*A*h^(0.7 + a)*t - k^0.7)^0.4285714286                                                                (3)

However, if I already have put beta=0.3 into the original differntial equation (see (1)) and solve it, I receive:

lambda(t) = _C1*exp(rho*t)/(10*k^(7/10) + 7*A*h^(7/10 + a)*t)^(3/7)                                                                       (4)

Clearly, the resultant solutions in (3) and (4) are quite different since the solution of (3) is any complex number and in (4) it is a real number. Note, that when solving the differential function (1) numerically, there exists a solution for any value of beta.

I also have tried to rename beta (supposing that maybe beta is any protected name) but the result remains the same. Also note, that if such manipulation is done with any other parameters, the solutions (3) and (4) are the same.

I will be thankful for any clarification regarding this issue. 

Thanks in advance,

Dmitry

My maplet

mpt := Maplet(Window("Точка выше или ниже прямой",
[[Plotter[f]()], ["1st vector"],
[TextField[v1](), TextField[v2](), TextField[v5]()],
["2nd vector"],
[TextField[v3](), TextField[v4](), TextField[v6]()], ["Скаляр"],
[TextField[scalar]()], [Button("Add", Evaluate(f = 'work(1)'))],
[[ToggleButton['d2']("2D", 'value' = true, 'group' = 'BG1')],
[ToggleButton['d3']("3D", 'value' = false, 'group' = 'BG1')]]]),
ButtonGroup['BG1']())

Now I choose via toogle button what to plot 2D or 3D addition but I need to plot it in two different plotter, is there a way to do it?

Procedure 

with(Maplets);
with(Elements);
with(Tools);
[AddAttribute, AddContent, Get, ListBoxSplit, Parse, Print, Set, 

  SetTimeout, StartEngine, StopEngine]
with(plots);
with(LinearAlgebra);
with(Student[LinearAlgebra]);

with(Maplets[Elements]);

work:=proc(g)

     local x1, y1, x2, y2, scalar;z1, z2,d2,d3,u,v,U,V,r1,r2,r3,R1,R2,R3,V1:
 x1 := Get(v1::algebraic);
 y1 := Get(v2::algebraic);  

z1 := Get(v5::algebraic);  

x2 := Get(v3::algebraic);  

y2 := Get(v4::algebraic);  

z2 := Get(v6::algebraic);  

scalar := Get(scalar::algebraic);

 d2:= Get(d2);  

d3:= Get(d3);    

 if d2 = true then  u, v:=<x1, y1>, <x2, y2>;  

elif d3 = true then   u, v:=<x1, y1,z1>, <x2, y2,z2>;  

end if;  U:=arrow(u, color=red):  

V:=arrow(u, v, color=black):    

   r1:=u+v: r2:=u-v: r3:=ScalarMultiply(u,scalar):r4:=ScalarMultiply(v,scalar):  

R1:=arrow(r1, color=red):

 R2:=arrow(v, r2, color=black):  

R3:=arrow(r3, color=black):

 R4:=arrow(r4, color=black):  

V1:=arrow(v, color=blue):    

if g = 1 then return VectorSumPlot(u,v, show=1, caption="Треугольник", vectorcolors=[red,blue]);    

end if;  

end proc:

Repeat once again what do I need: I want to plot addition of vectors in 2D and 3D in 2 different plotters at once

https://www.maplesoft.com/applications/view.aspx?SID=3938&view=html

i find this example

how to write piano note A and B frequency music to wave file?

how to convert a function to transition matrix?

how to convert a differential system to transition matrix?

Hi everybody and thank you all in advanced.

I tried this:

1. convert([a,b],`*`) # this worked
2. convert([a,b],`/`) # this gives an error
3. map(`/`,a,b) # this worked

Questions:

1. How can I make code for point (2) work? Apparently `/` is not an option for convert.
2. I was looking for a function that divide two symbols, something like divide(a,b) and outputs a/b for some simple
3. programing purposes but find no one. (a,b) could be numbers, could be trig functions or polynomials.
   
   Of course, you can just input a/b or lis[i]/lis[i+n], but a knowledgeable answer to this will cast light over Maple’s ways and my wrong ways.

Let be the sequence (an) so that a1 = 1, a(n+1) = 5*a(n) + sqrt(k*a(n)^2 -8), for all n >=1. I know that, the answer is k = 24.  How can I choice number k so that the sequence is an integer sequence?

I am trying to export the data list to .dat but i received that error. New with coding and Maple.

restart;
a := sin(x^(1/3))+sin(x)^(1/3);
                      / (1/3)\         (1/3)
                   sin\x     / + sin(x)     
b := diff(a, x);
                     / (1/3)\                
                  cos\x     /      cos(x)    
                  ----------- + -------------
                      (2/3)             (2/3)
                   3 x          3 sin(x)     
output := fopen("Ospino.dat", WRITE);
for x by .5 to 50 do fprintf(output, "%f\t %f\t \n", x, b) end do;
Error, (in fprintf) number expected for floating point format
fclose(output);
 

add(q[k],k=1..5)

I'd like to be able to display a general result such as

q_1 + q_2 + ... + q_n

ideally maple would do all the special stuff for me, e.g., somehow I could do

add(q[k],k=1..2..n)

but I can't even seem to get a nice output manually that doesn't involve a lot of printing code.

Hi

Can maple compute the following integral

Integral ( z^(-0.5), z in M) where  M is the lower half-unit circle from +1 to -1

Many thanks

How I can apply the condition or assuming (a>0) in ?

I want to have a plot according to the attached figure.

Download plot 

Normally the wheel mouse will scroll a file like it does in just about every other program... but when the mouse is over a code edit region it does nothing ;/ As one is scrolling with the wheel and a code edit region pops up this breaks the flow and requires hoops, I don't like hoops!

Seems most mouse funcitons do not work well in the region. Just trying to select and drag doesn't scroll the region.

I was writing in a maple document when suddently all the text and equations collapsed into a single line at the top, making it unreadable and unusuable. I was not able to restore a previous backup of the file, however i did manage to salvage a smaller part of the document. Both files are included. 
corrupted.mw 
what_i_managed_to_save.mw

ModularSquareRoot can't handle a really big n because of the factorization problem. But I would like to know if there is any other function than can use as parameters the factorization of n.

I already have an algorithm that does that for n=pq with p, q primers, and could generalize it, but I was wondering whether there was a built-in function that does that.