Hi everyone, I have a problem in the code solving coupled partial differential equations. I could not find out the solution. Please help me out with this. Find the code in the attachment.

Dear maple users,
Greetings.
I am solving an ode problem with an analytical solution.
programming running properly, but my plot not exact with the already existing article plot. 
how to get the exact plot.

Thanking you.

Code:JVB.mw
 

Download JVB.mw

Analytical solution approach:

Explore the values of km_digit(n,m) using km_ for all m,0 < or equal to m and less than or equal to 8. 

Look at the output until you can make a conjecture that concerns the pattern obtained.

Hint use km_list(m,6,20) when m is not a multiple of 4 and km_list(m,6,50) when m is a multiple of 4

I am studying the motion of a beam coupled to piezoelectric strips. This continuous system is modelled by two DE:

YI*diff(w(x,t), x$4)-N[0]*cos(2*omega*t)*diff(w(x,t), x$2)+c*diff(w(x,t), t)+`ρA`*diff(w(x,t), t$2)+C[em,m]*v(t) = 0;

And:

C[p]*diff(v(t), t)+1/R[l]*v(t) = C[em,e]*(D[1,2](w)(0,t)-D[1,2](w)(ell,t));
 

where "w(x,t)" stands for the beam's vibration and "v(t)" means the electric voltage, which is constant throught the beam. I would like to numerically solve both DE simultaneosly, but maple will not let me do it. I would like to know why. I am getting the following error:

Error, (in pdsolve/numeric/process_PDEs) number of dependent variables and number of PDE must be the same
 

I suppose it is because "w(x,t)" depends on "x" and "t", while "v(t)" depends solely on time, but I am not sure. Could someone help me out? Here is my current code:

restart:
with(PDEtools):
declare(w(x,t), v(t)):

YI*diff(w(x,t), x$4)-N[0]*cos(2*omega*t)*diff(w(x,t), x$2)+c*diff(w(x,t), t)+`ρA`*diff(w(x,t), t$2)+C[em,m]*v(t) = 0;
pde1:= subs([YI = 1e4, N[0] = 5e3, c = 300, omega = 3.2233993, C[em,m] = 1], %):
ibc1:= w(0,t) = 0, D[1,1](w)(0,t) = 0, w(ell,t) = 0, D[1,1](w)(ell,t) = 0, D[2](w)(x,0) = 0, w(x,0) = sin(Pi*x/ell):

C[p]*diff(v(t), t)+1/R[l]*v(t) = C[em,e]*(D[1,2](w)(0,t)-D[1,2](w)(ell,t));
pde2:= subs([C[p] = 10, R[l] = 1000, C[em,e] = 1, ell = 5], %):
ibc2:= v(0) = 0:

pdsolve({pde1, pde2}, {ibc1, ibc2}, numeric);

Thanks.

Hello everyone!

I'm interesting in "zcoloring" funciton in colorscheme option.

I wrote simple programm which compares two results: spectrogram of signal drawn with "colormap" list and spectrogram which was plotted with zcoloring function. I use red, green, blue functions to construct JET-colormap: list and expressions in "zcoloring".

My result:

As I understand, when I use:

colorscheme = ["zcoloring", [z-> Red color function(z), z-> Green color function(z), z-> Blue color function(z)], colorspace = "RGB"]

Maple plots z-value with color RGB color coordinates defined from "zcoloring". For example, if "zcoloring" is

colorscheme = ["zcoloring", [z-> 5*z, z-> 3*z, z-> 2*z], colorspace = "RGB"]

and z value is 10, then 10 value will correspond [50,30,20]-RGB color.

My test program:

Spectrogram_zcoloring.mw

Spectrogram of my test signal:

list_test.txt

Hello,

a follow up question.

I am solving some overdetermined system of ODEs in cylindrical coordinates r,phi,Z. I obtain some equations of the following type:

(diff(_F1(phi, Z), phi)*r + diff(diff(s_r(phi, Z), phi), phi))/r = -s_r(phi, Z)/r

As can be seen, the differentiated functions do not depend on r, which is an independent variable. Thus, the correct solution is to separate the equation and have 

_F1(phi,Z)=_F1(Z), s_r(phi,Z)=s_r(Z).

By using dsolve, I always obtain a solution containing r.

A similar problem that does no contain derivatives is solved by solve/identity.

Is there something similar for dsolve?

EDIT: I again put here more info and file. I solve some overdetermined system of differential equations.

[diff(s_r(r, phi, Z), r) = 0, diff(s_r(r, phi, Z), phi) = -diff(s_phi(r, phi, Z), r)*r^2,
diff(s_Z(r, phi, Z), r) = -diff(s_r(r, phi, Z), Z), diff(s_phi(r, phi, Z), phi) = -s_r(r, phi, Z)/r, 
diff(s_Z(r, phi, Z), phi) = -diff(s_phi(r, phi, Z), Z)*r^2, diff(s_Z(r, phi, Z), Z) = 0, 
diff(m(r, phi, Z), r) = s_Z(r, phi, Z)*B_phi(r, phi, Z) - s_phi(r, phi, Z)*B_Z(r, phi, Z),
 diff(m(r, phi, Z), phi) = s_r(r, phi, Z)*B_Z(r, phi, Z) - s_Z(r, phi, Z)*B_r(r, phi, Z),
 diff(m(r, phi, Z), Z) = s_phi(r, phi, Z)*B_r(r, phi, Z) - s_r(r, phi, Z)*B_phi(r, phi, Z), 
s_r(r, phi, Z)*diff(W(r, phi, Z), r) + s_phi(r, phi, Z)*diff(W(r, phi, Z), phi) + 
s_Z(r, phi, Z)*diff(W(r, phi, Z), Z) = 0]

After some time, I arrive at the equation in the original question. So the independance of the other functions on r is the consequence of the other equations.

Here is the file (shortened): mwquestion2.mw

Sorry if it's a dumb question, i'm new with Sage.

I'm doing the multiplication of two matrix and I obtain the expression 

[142√(2x−2√)+122√y] 

I want to represent this on 3d, but I realized that the expression has to be written with operators '*'

How can I add the operators to my expression?

   1 1 1 2 2 3 3  

(  1 3 3 2 2 1 2 )

how to input this permutation into permgroup ?

Who deleted www.mapleprimes.com/questions/228364-Graphics-Problem

and www.mapleprimes.com/posts/211539-Buggy-Lights

I thought I had alzheimers or something.  Who's deleting these posts and questions.  I realize they are similar but there is different information in there that was useful.  Why don't we delete all similar questions .. of course this would be horrible since some have slight differences - granted some don't have answers tagged on the end and are exact copies those are of course waranted and I would agree to be deleted becuase they were duplicates.

Please undelete those questions - there were referenses and extensive conversations in there that are not present in the saving question.

This is another problem I just found in Maple 2019.2 on windows 10. professional.

I wanted to close Maple, so did  File->Exit 

But Maple did nothing. It did not close.  Also Alt-F4 did not close Maple. I had to click on the little X on top right corner of the open window to close Maple.  

In earlier version this used to work to close Maple.

Do others see this as well?  To reproduce, simply start Maple, and do File->Exit.

Here is a movie also

This may be a bug.  In Maple 2019.0

plot3d(x^2+y^2)

Now grab the graphing window and drag to adjust it's size.  It doesn't adjust it's size until you let go of the mouse button.

If maple can auto discover derived equations, 

then most variable are unknown name which may or may not be physics variable.

so, how to guess which physics data suitable for these look like meaningless variable?

or these tools only visualize the relationship of existing known equations?

How to solve this DE with IC by using DTM.
D^m u(x,t)=u''(x,t)-u^2 (x,t), where n-1< m < n 

IC: u(x,0)=1+sin(x), and u'(x,0)=0

I have a problem for school that I need help with. 

Solve: by variation of parameters.

x³y’’’ – x²y’’ − 2xy’ + 6y = x²

Use any software, e.g. Maple, as an aid in computing roots of the auxiliary equation and the Wronskian-based determinants W, W1, W2, W3

I was able to hand calculate the roots as [m₁ = 3, m₂ = 2, m₃ = -1]
 

Download project_27.mw

I have just upgraded my laptop from Windows 7 to Windows 10.  On starting up Maple 2018,  I receive the attached message on screen.  This is after previously loading the worksheet successfully.   Today,  I am not able to do so.  I need to permanently register my firewall to allow Maple to run; can anyone help?

Melvin