MaplePrimes Questions

I am running MAPLE 2016.1

I am able to set the color option to a function in fieldplot, e.g.

plots[fieldplot]([y, -sin(x)-(1/10)*y], x = -10 .. 10, y = -10 .. 10, arrows = SLIM, color = sin(x)*sin(y))

works as expected.  

However, for plots[fieldplot3d] I cannot set color = f(x,y,z) . E.g.

plots[fieldplot3d]([2*x, 2*y, 2*z], x = -1 .. 1, y = -1 .. 1, z = -1 .. 1, grid = [5, 5, 5], arrows = SLIM, color = sin(x)*sin(y)*sin(z));

Error, (in plot/color) invalid color specification: sin(x)*sin(y)*sin(z)

Please can anyone help?

MRB  

First Question: How to define nx1 matrix Y:=(y1,y2,...,yn) ? (n is a Natural number while it is ungiven)

Second Question: How to derivative of matrix Y with respect to nx1 matrix X:=(x1,x2,...,xn) ?

 


My efforts for the first question:  (I know to define 6x1 matrix etc. , but I dont know to define nx1 matrix

restart: Matrix(1..6,1,symbol=y) 

 

My efforts for the second question:

restart; with(VectorCalculus);
Matrix([x^2, x*y, x*z]);
Jacobian([x^2, x*y, x*z], [x, y, z]);


Can you help me? 

Hello

I am trying to write my Measuring Data into several Arrays with lengt of rowdim

I was able to extract them from the Datafile and Refering them to separate Vectors

My Code is in the attachemend and the Datafile too.

Insert the txt Path in the Messdaten in definition of A to load it.

I was able to convert the Time from Seconds to Minutes. What i realy want is an x-Axis wich is like hh:mm:ss

and Plot for Example T__1,T__2,T__Scheibe together in one Plot as y Axis and t__Minuten as x Axis

But the

plot([seq([`t__Sekunden`[i],`T__1`[i]],i=1..rowdim,LinearAlgebra[rowdim](`t__Sekunden`))

doesnt work and i also tried others

Thank you for your Help

 

Hallo,

im currently using Mathcad 15 and i want to change to a newer and better software with more possibilities.

But up to now i have not found a better software for calculating. One big advantage with mathcad is the possibilitie of symbolic formula input and calculation with units.

Now my question: Is it possible with Maple to write symbolic formulas (2D Structure of big formulas)

I dont write a formula in one row. Its nearly impossible ...

And can i calculate with units?

Thx Stefan

 

I wrote an equation and intended to iterate P(r) till that value of r when P(r) becomes 0.what is the value of r for which P(r)=0, and also  obtained its convergent decreasing graph which show convergence of P(r) to 0.
 

restart

NULL

G := 6.6743*10^(-8);

0.6674300000e-7

 

0.2997924580e11

 

0.7321400000e36

(1)

``

sys2 := diff(P(r), r) = -G*(epsilon+P(r))*((8.98*10^14*(1/3))*Pi*r^3+4*Pi*r^3*P(r)/c^2)/(c^2*(r^2-(2*G*r*(8.98*10^14*(1/3))*Pi*(r^3))/c^2)), P(0.1e-9) = 8.5561*10^34

diff(P(r), r) = -0.7426160269e-28*(0.7321400000e36+P(r))*(0.2993333333e15*Pi*r^3+0.4450600224e-20*Pi*r^3*P(r))/(r^2-0.4445794614e-13*r^4*Pi), P(0.1e-9) = 0.8556100000e35

(2)

``

``


 

Download decreasing_graph.mw

Error, More than two dependent variables, please indicate the desired scene

Trace an ellipse that has 3 or more foci using Maple.

As a geotechnical engineer, much of what I calculate relates to depth or settlement.  Which, by tradition, is signed as positive in the vertically down direction.  Hence, I want to be able to produce graphs with positive y-values (e.g. depth, settlement, etc.) plotted in the downward direction (effectively, all I want to do is simply reverse the direction of the y-axis).  One workaround is to simply multiply all y values by -1 before plotting.  But this isn't ideal (as the values show as negative).  I've had a look through the documentation and haven't found a clean way of doing this.  Is anyone able to confirm whether its possible?  Or am I wasting my time looking?  Thanks in anticipation, Ian

I can't cut and paste from my maple worksheet (because syntax somehow gets changed in the process), so perhaps the following attempt at carefully worded  prose will do (alternately, please someone tell me how to attach screen shots).  Here goes...

I execute

The optimization succeeds and the assignment takes place.  The result is what I'ld like to cut and paste so youall can see...  it looks as follows. 'r:=' (but no quotes)  followed by open square bracket, then a number, then a comma, then another open square bracket, then a vertically arranged list of stuff, then two close brackets.  I presume the vertically arranged stuff in the innermost matching square brackets somehow represents the point in the search space where the minimum takes place, and the minimum is the number that follows the first square bracket and is before the comma.  My goal is to write to disk the the point in the search space where the minimum takes place (I suppose that vector is in some sense r[2]) .  My attempt to write it to disk is

ExportMatrix("maple.txt", r[2], target = auto)

but maple complains

ExportMatrix("maple.txt", r, target = auto);
Error, invalid input: ExportMatrix expects its 2nd argument, M, to be of type {Matrix, list(Matrix)}, but received [293148935266.099426, Vector(28, {(1) = -233674.34818878854, (2) = 12.62636610207248, (3) = 0.3337043409389546e-2, (4) = 0.13303915611938133e-1, (5) = .0, (6) = .2589485687872699, (7) = 0.9341600279762804e-1, (8) = 0.3293862590239189e-1, (9) = 0.12892649003200331e-1, (10) = 0.606528143324688e-1, (11) = 0.17524827222051658e-1, (12) = -0.16059611469561555e-17, (13) = 0.621465506367271e-1, (14) = .29826235341855, (15) = -0.16288776533795546e-17, (16) = 0.15301368661023716e-1, (17) = 0.1529933711903137e-18, (18) = 0.4375115591744606...

How can I write the result of optimization to disk as a plain ascii text file (tsv format)?

Thanks

 

 

Hey there
 
I'm Trying to install maple on my Linux laptop, i have it install and starting fine but, i cant activate it. I get the massage " This Purchase code is only valid for Maple 2016 " But I know that it is the right Activation code, On a Side note I'm a student so the code is form the my school.
 

hope you can understand it.

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

1 2 3 4 5 6 7 Last Page 3 of 1394