MaplePrimes Questions

I'm looking for the Maplets package for Maple 7 that was available on the old mapleprimes.  Does anyone have it?

I did not realize Maple 2016, since Maple 2015 actually, had stopped Windows XP support. 

hi

how i can remove this error?

1) ''Error, invalid input: lhs received diff(diff(diff(diff(w(x, y, t), x), x), x), x)+.12*(diff(diff(diff(diff(w(x, y, t), x), x), y), y))+0.11e-2*(diff(diff(diff(diff(w(x, y, t), y), y), y), y))+10-.4*(diff(diff(w(x, y, t), x), x))-0.4e-2*(diff(diff(w(x, y, t), y), y))+diff(diff(w(x, y, t), t), t)-w(x, y, t)-tau, which is not valid for its 1st argument, expr'''''

2)'' Error, incorrect number of extra arguments in select''

thanks

scale.mw
 

 

NULL

restart:Digits := 15: beta := 1: alpha := 0.1: Upsilon := .1: xi := 10: eta := .1: N_X := .4: N_Y := .4: psi:NULL=4.73:tau = 10

tau = 10

(1)

A5 := 1; A6 := 2*alpha^2; A7 := alpha^4; A8 := xi; A9 := xi*alpha^2; A10 := xi*alpha^4; A11 := -N_X; A12 := -N_Y*alpha^2; A13 := 1; A14 := -beta; A15 := -tau

-tau

(2)

EOM := A5*(diff(w(x, y, t), x, x, x, x))+A6*(diff(w(x, y, t), x, x, y, y))+A7*(diff(w(x, y, t), y, y, y, y))+A8(diff(w(x, y, t), x, x, x, x))+A9*(diff(w(x, y, t), x, x, y, y))+A10*(diff(w(x, y, t), y, y, y, y))+A11*(diff(w(x, y, t), x, x))+A12*(diff(w(x, y, t), y, y))+A13*(diff(w(x, y, t), t, t))+A14*w(x, y, t)+A15

diff(diff(diff(diff(w(x, y, t), x), x), x), x)+.12*(diff(diff(diff(diff(w(x, y, t), x), x), y), y))+0.11e-2*(diff(diff(diff(diff(w(x, y, t), y), y), y), y))+10-.4*(diff(diff(w(x, y, t), x), x))-0.4e-2*(diff(diff(w(x, y, t), y), y))+diff(diff(w(x, y, t), t), t)-w(x, y, t)-tau

(3)

BC := w(0, y, t) = 0, (D[1, 1](w))(0, y, t) = 0, w(a, y, t) = 0, (D[1, 1](w))(a, y, t) = 0, w(x, 0, t) = 0, (D[2, 2](w))(x, 0, t) = 0, w(x, b, t) = 0, (D[2, 2](w))(x, b, t) = 0

w(0, y, t) = 0, (D[1, 1](w))(0, y, t) = 0, w(a, y, t) = 0, (D[1, 1](w))(a, y, t) = 0, w(x, 0, t) = 0, (D[2, 2](w))(x, 0, t) = 0, w(x, b, t) = 0, (D[2, 2](w))(x, b, t) = 0

(4)

 

 

BC1:=subs(a=1,b=1,[BC])

[w(0, y, t) = 0, (D[1, 1](w))(0, y, t) = 0, w(1, y, t) = 0, (D[1, 1](w))(1, y, t) = 0, w(x, 0, t) = 0, (D[2, 2](w))(x, 0, t) = 0, w(x, 1, t) = 0, (D[2, 2](w))(x, 1, t) = 0]

(5)

We use the method of multiple scales to directly attack EOM1 and BC1. To transform the time derivatives in EOM1 in terms of the scales  and , we let

timeScales:=T[0],T[1]

T[0], T[1]

(6)

dt[1]:=expr->add(epsilon^i*diff(expr,timeScales[i+1]),i=0..1)

proc (expr) options operator, arrow; add(epsilon^i*(diff(expr, timeScales[i+1])), i = 0 .. 1) end proc

(7)

msForm:={w(x,y,t)=w(x,y,timeScales),seq(diff(w(x,y,t),t$i)=dt[i](w(x,y,timeScales)),i=1..2)}

{diff(diff(w(x, y, t), t), t) = dt[2](w(x, y, T[0], T[1])), diff(w(x, y, t), t) = diff(w(x, y, T[0], T[1]), T[0])+epsilon*(diff(w(x, y, T[0], T[1]), T[1])), w(x, y, t) = w(x, y, T[0], T[1])}

(8)

multiScales:=`union`(map(s->subs(w=s,msForm),[w])[])

{diff(diff(w(x, y, t), t), t) = dt[2](w(x, y, T[0], T[1])), diff(w(x, y, t), t) = diff(w(x, y, T[0], T[1]), T[0])+epsilon*(diff(w(x, y, T[0], T[1]), T[1])), w(x, y, t) = w(x, y, T[0], T[1])}

(9)

Then, we seek a second-order approximate solution in the form

solRule:=w(x,y,timeScales)=add(epsilon^j*w[j](x,y,timeScales),j=1..3)

w(x, y, T[0], T[1]) = epsilon*w[1](x, y, T[0], T[1])+epsilon^2*w[2](x, y, T[0], T[1])+epsilon^3*w[3](x, y, T[0], T[1])

(10)

We introduce the detuning  to describe the nearness of the edge lengths  and , consider the case of primary resonance, and hence let

tau=epsilon^3*tau(x,y)*convert(cos(Omega*T[0]),exp)

tau = epsilon^3*tau(x, y)*((1/2)*exp(I*Omega*T[0])+(1/2)*exp(-I*Omega*T[0]))

(11)

lhs(EOM)

Error, invalid input: lhs received diff(diff(diff(diff(w(x, y, t), x), x), x), x)+.12*(diff(diff(diff(diff(w(x, y, t), x), x), y), y))+0.11e-2*(diff(diff(diff(diff(w(x, y, t), y), y), y), y))+10-.4*(diff(diff(w(x, y, t), x), x))-0.4e-2*(diff(diff(w(x, y, t), y), y))+diff(diff(w(x, y, t), t), t)-w(x, y, t)-tau, which is not valid for its 1st argument, expr

 

Substituting multiScales, solRule, and scaleRule2 into EOM1, expanding the result for small , and discarding terms of order higher than , we obtain

expr83a:=seq(convert(series(value(subs(multiScales,solRule,lhs(EOM))),epsilon,4),polynom),i=1)

Error, invalid input: lhs received diff(diff(diff(diff(w(x, y, t), x), x), x), x)+.12*(diff(diff(diff(diff(w(x, y, t), x), x), y), y))+0.11e-2*(diff(diff(diff(diff(w(x, y, t), y), y), y), y))+10-.4*(diff(diff(w(x, y, t), x), x))-0.4e-2*(diff(diff(w(x, y, t), y), y))+diff(diff(w(x, y, t), t), t)-w(x, y, t)-tau, which is not valid for its 1st argument, expr

 

NULL

NULL

t0 := time(); ST1_like := map(combine, select(has, select(has, ST, sin(m*Pi*x)), sin(n*Pi*y))); time_taken := time()-t0

.187

 

Error, incorrect number of extra arguments in select

 

0.

(12)

``


 

Download scale.mw

 

Hello,


Let matrix A be nxn matrix. How to find modal matrix M, the modal matrix for the matrix A is the n × n matrix formed with the eigenvectors of A as columns in M? (D=M^-1*A*M, where D is dioganal matrix with the eigenvalues of A on the main diagonal of D)

For example;

restart;
with(LinearAlgebra):
A := Matrix(2, 2, [-4, -2, 3, 1])

How to find M modal matrix?

Goodmorning,

 

Has the command DGinfo(T,"CoefficientSet") an equivalent in the Physics package ? is there another way to extract the non-zero elements from the Einstein tensor in a format acceptabele to pdsolve?

 

Many thanks,

 

Luc Van Gompel

Hi.

In some cases when dealing with vectofields an such the are integral has to be expressed in terms of r(t).

the general form for r is r^2=(r*cos(t)-a)^2+(r*sin(t)-b)^2, When I solve this in maple it seems like I get the inverse of the desired result.

If I knew that was always the case I could just inverse my result to get the right expression for r, but im not sure if it only applies for this particular cas or all cases.

I would be happy if anyone took a quick look and suggested a way to obtain the desired solution for any center (a,b) for the circle.

 


 

Expression for radius, circle centred at (a,b)

RA := r^2 = (r*cos(t)-a)^2+(r*sin(t)-b)^2

r^2 = (r*cos(t)-a)^2+(r*sin(t)-b)^2

(1)

isolate(r^2 = (r*cos(t)-a)^2+(r*sin(t)-b)^2, r)

r = (1/2)*(a^2+b^2)/(cos(t)*a+b*sin(t))

(2)

eval(%, [a = -1, b = 0])

r = -(1/2)/cos(t)

(3)

plot3d([-2*x, x^2+y^2], y = -sqrt(-x^2-2*x) .. sqrt(-x^2-2*x), x = -2 .. 0, color = [green, red], orientation = [0, 0, 0])

 

(1/2)*Pi <= t and t <= 3*Pi*(1/2)

(1/2)*Pi <= t and t <= (3/2)*Pi

(4)

0 <= r and r <= -2*cos(t)

0 <= r and r <= -2*cos(t)

(5)

Area_off_center = int(r, [r = 0 .. -2*cos(t), t = (1/2)*Pi .. 3*Pi*(1/2)]); 1; Area_at_center = int(r, [r = 0 .. 1, t = 0 .. 2*Pi])

Area_off_center = Pi

 

Area_at_center = Pi

(6)

``

``


 

Download parametrization_of_r_not_centred_at_orgin.mw
 

I would also happily like to know how I can solve for the range r can take, obviously in the example i´m working with here r starts at 0, but that is not always the case i guess.

 

Thank you, your help is much apperciated

Hello

I changed the color of the background of a worksheet in black (I introduced a Table for that reason) because the white color hurts my eyes. I still have a couple of things that I'd like to modify:

1. The blinking cursor is black and is now invisible. How can I change its color?

2. The default input color is red ; i'd like to change it in white. How can I do?

Thanks Nicola

How do I genereate random negative numbers?

 

when i put a minus sign in the rand-function i get an error.

 

I need to for example generate a random number between -r and r, where r is the input

 

 

Is there a difference between scoping rules for regular procedures and for arrow operators?

      Inspired by the theme
http://www.mapleprimes.com/questions/219995-Finding-A-Convinient-Parametrization-Of-Surfaces
Examples in the Mathematica did Alexander Bannikov.
It is equidistant radius 0.1 to the surface

   (x1 ^ 2 + x2 ^ 2-0.4) ^ 2 + (x3 + sin (x1 * x2 + x3)) ^ 4-0.1 = 0;

https://vk.com/doc7819263_439405418?hash=af46d61d8aad95f70b&dl=9f245f5b6b68b47075

and an example of parameterization the same surface

https://vk.com/doc7819263_439432143?hash=36cf31d52c97e2e373&dl=7e4fa17a771dffb331

As I have understood from the words of Alexander Bannikov, parameterization was performed using the functions: RegionFunction, ContourPlot3D, ClippingPlanes.

It turns out that Maple functions inferior?

Hi everybody, 

I am working with Maple 2015.2 ; system = OS X El Capitan
I would like to compare the outputs of two files by using  the Worksheet:-Comparator command

OriginalFile.mw file contains this single instruction
> a:=1;

ModifiedFile.mw contains this single instruction
> a := 2;

Both us are located on my desktop, which is also my current directory in the worksheet from which I am going to run the comparison.

MyDir := cat(currentdir(), "/Desktop/");
OriginalFile := cat(MyDir, "OriginalFile.mw"):
ModifiedFile := cat(MyDir, "ModifiedFile.mw"):

# to verify that all is going well I use these two commands 
# Worksheet:-DisplayFile(OriginalFile):
# Worksheet:-DisplayFile(ModifiedFile):

Running the  command  Worksheet:-Comparator(OriginalFile, ModifiedFile) returns me error messages :
The first one is  (after having selected  "compare" and next "run" in the Comparator Maplet) 

and the second one (after having cliked ont the "OK" button" on this alert window)

Note that the same thing happens if I use Worksheet:-Comparator() and select the two files with the browsers
Could you please help me to understand what happens ?

Thanks in advance

 

How do I solve the problem that attached? Please give me the correction. thank you

Hello! Hope everything fine with you. I am try to solve the three-point differential by numerical method in attached file but failed. Please see the attachement and solve my problem. I am very thankful your kind effort. Please take care.

three-point.mw

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

What code to use in a maplet that would open the open file window? 

ie when I push a button the open file window opens allowing the user to choose a file.


 

refresh

refresh

(1)

G := 6.6743*10^(-8); 1; c := 2.99792458*10^10; 1; pi := 3.143; 1; rho := 5.3808*10^14

0.6674300000e-7

 

0.2997924580e11

 

3.143

 

0.5380800000e15

(2)

diff(P(r), r) = -G*(rho*c^2+P(r))*((4*pi*r^3*(1/3))*rho+4*Pi*r^3*P(r)/c^2)/(c^2*(r^2-2*G*r*(4*pi*r^3*(1/3))*rho/c^2)), diff(v(r), r) = 1.485232054*10^(-28)*((4*pi*r^3*(1/3))*rho+4.450600224*10^(-21)*Pi*r^3*P(r))/(r^2-1.485232054*10^(-28)*r*(4*pi*r^3*(1/3))*rho)

diff(P(r), r) = -0.7426160269e-28*(0.4836021866e36+P(r))*(0.2254913920e16*r^3+0.4450600224e-20*Pi*r^3*P(r))/(r^2-0.3349070432e-12*r^4), diff(v(r), r) = 0.1485232054e-27*(0.2254913920e16*r^3+0.4450600224e-20*Pi*r^3*P(r))/(r^2-0.3349070432e-12*r^4)

(3)

condition; -1; P(0) = 0, v(1014030) = -.4283

P(0) = 0, v(1014030) = -.4283

(4)

``


 

Download maple_soft.mw

I found the solution of P(r) at P(0)=0, but could obtain the result of v(r) at v(1014030)=-0.4283, v(r) may have a graph such that i can goes from -0.4283 to 0.

First 918 919 920 921 922 923 924 Last Page 920 of 2308