Advanced Search

Maple 2016 Questions and Posts

These are Posts and Questions associated with the product, Maple 2016
View unanswered posts
Maple 2016 Questions and Posts Feed

Error, (in unknown) invalid subscript selector...

Asked by: illuminates 65
February 08 2017
0  2

Hello everyone, I use physics package and I input:

with(Physics)

Setup(mathematicalnotation = true)

Coordinates(X)

Setup(tensors = A[mu](X))

PDEtools:-declare(A(X))

PDEtools:-declare(p0(X), pm(X), pp(X), pt(X), Un(X), Ucn(X))

pt := proc (X) options operator, arrow; Matrix(2, 2, [[p0(X), sqrt(2)*pp(X)], [sqrt(2)*pm(X), -p0(X)]]) end proc

U := proc (X) options operator, arrow; Matrix(2, 2, [[1, 0], [0, 1]])+I*pt(X)*(1/f)-pt(X)*pt(X)*((1/2)/f^2) end proc

Ucn := proc (X) options operator, arrow; Matrix(2, 2, [[1, 0], [0, 1]])-I*pt(X)*(1/f)-pt(X) . pt(X)*((1/2)/f^2) end proc

And after that I write:

Trace(-LeviCivita[mu, nu, rho, sigma] . Ucn(X) . ((1/6)*e*A[nu](X)+M5*KroneckerDelta[nu, 0]) . Matrix(2, 2, [[1, 0], [0, 1]]) . Uсn(X) . d_[rho](Un(X)) . Uсn(X) . d_[sigma](Un(X)))

But I get an error:

Error, (in unknown) invalid subscript selector

Thanks for the help.

 

How to differentiate a plot containing numerical s...

Asked by: gihan 5
differential-equation odeplot
February 07 2017
1  4

Hi,

I have a Maple code which produces an output plot for a first order differential equation,

nde := evalf(subs(npar, de));
nds := dsolve({nde, sigma(0) = -1e-8}, sigma(t), type = numeric);
acc_nds := (sig0, ae, re) -> dsolve({nde, sigma(0) = sig0}, sigma(t), type = numeric,
method = lsode[backfull], abserr = ae, relerr = re, maxfun = 0, ctrl=Ctrl);

odeplot(acc_nds(-0e-7, 1e-13, 1e-13), [t, sigma(t)], t = 0..2);

This produces the outputplot that I need for sigma(t). I need to produce a outut for d(sigmat)/dt, and how can this be done? what is the command I should use?

Additionally, how can i get the data set out of the polt in to a excel file or a text file?

I am quite new to maple, so i expect your kind support

Thanks

Matrices in Physics package...

Asked by: illuminates 65
matrix physics
February 07 2017
0  1

Hello! I try to make friends physics package and matrixes. But I am faced with difficulties. To be more specific, to consider a code:

with(Physics):
Setup(mathematicalnotation = true);
Coordinates(X);
Setup(tensors = A[mu](X));
PDEtools:-declare(A(X), p0(X), pm(X), pp(X), pt(X), U(X), m5(X))

And then I create a matrix:

pt := proc (X) options operator, arrow; Matrix(2, 2, [[p0(X), sqrt(2)*pp(X)], [sqrt(2)*pm(X), -p0(X)]]) end proc;

It displays as:

Why do p0, pp, pm become function only one variable x1?

Interestingly that maple understands the next matrix:

Nevertheless even for function only one variable derivative works fine:

And I can construct more complicated functions:

But why is only one variable?

Lazy caterer's animation...

Asked by: brian bovril 889
animation
February 07 2017
1  8

Edited.

Fellas

I'm wondering how to produce the graphic in the

https://en.wikipedia.org/wiki/Lazy_caterer's_sequence.

my attempt:
caterers.mw

How do i strip my solve() result for the values on...

Asked by: pizzaluc 20
February 06 2017
1  2

I will try to explain with an example, lets say i have to solve 2 equations with 2 unknowns

P := solve({3*y+x, 2*x+3+y})
solve({3*y+x, 2*x+3+y})

I would then like to do P[1] and P[2] to get the values, but instead it gives me "x=-9/5" which make me unable to use it for things like

P[1] + P[2]

Hope you all understand, and thanks for taking your time to read my problem 

Is it possible to import questions from Maple TA i...

Asked by: mathstutor 35
export import
February 06 2017
0  4

Hi, I'd like to know if is it possible to import questions from Maple TA into Maple and use them in quizes, for example. Your help will be highly appreciated. Thank you.

I can't get solve to work...

Asked by: TheGamingEngineer 0
solve
February 06 2017
0  2

i have these commands:

restart;
with(plots);
with(VectorCalculus);
with(LinearAlgebra);
f := A-B+x^y/y;

g := x^y-y;
Solve({f = 0, g = 0}, {x, y});
              / /                     y    \         \
              | | y                  x     |         |
         Solve|< x  - y = 0, A - B + -- = 0 >, {x, y}|
              | |                    y     |         |
              \ \                          /         /

restart; with(plots); with(VectorCalculus); with(LinearAlgebra); f := A-B+x^y/y; g := x^y-y; Solve({f = 0, g = 0}, {x, y});

 

but I can't get it to actually solve the equations. what do I do  ?

solve ode srm motor...

Asked by: lucaud 35
February 06 2017
0  2

How I can solve it for P?

P=i^(2)r+(&DifferentialD;)/(&DifferentialD; t) (1/(2)l(tetha)i^(2))+1/(2)i^(2)(&DifferentialD;l(tetha))/(&DifferentialD; theta) w

attach i(t) "corrente", l(t) "induttanza", theta "angular", w "rotary speed"graph of funtions

thanks

Symbolic cut-off of high-order terms...

Asked by: illuminates 65
polynomial
February 05 2017
0  4

Hi! I have a variable polynomial expression and I want to cut off all the terms of order 3 and more, for every product of variables.

For example consider the polynomial P = x + y + a*x^2 + b*x^3 + c*x*y + d*x^2*y + e*x*y^2 + y^2, I want an operation who returns me a*x^2 + c*x*y + y^2

(terms for which the sum of exponents in x and y exceeds or less two are being deleted)

Is it possible?

 

how i can write this line ''for'' in maple?...

Asked by: 9009134 110
February 04 2017
0  2

hi..

how i can write this line ''for'' in maple?

thanks..

for.mw
 

TopEdge = "C";
BottomEdge = "C";
LeftEdge = "C";
RightEdge = "C";

"For[i = 1, i < m + 2, i++, w[i, 1] = 0; w[i, n + 1] = 0]; For[j = 1, j < n + 2, j++, w[1, j] = 0; w[m + 1, j] = 0]; If[TopEdge == "C", For[i = 1, i < m + 2, i++, w[i, 0] = w[i, 2]], If[TopEdge == "S", For[i = 1, i < m + 2, i++, w[i, 0] = -w[i, 2]], "Invalid Input Data!"]] If[BottomEdge == "C", For[i = 1, i < m + 2, i++, w[i, n + 2] = w[i, n]], If[BottomEdge == "S", For[i = 1, i < m + 2, i++, w[i, n + 2] = -w[i, n]], "Invalid Input Data!"]] If[LeftEdge == "C", For[j = 1, j < n + 2, j++, w[0, j] = w[2, j]], If[LeftEdge == "S", For[j = 1, j < n + 2, j++, w[0, j] = -w[2, j]], "Invalid Input Data!"]] If[RightEdge == "C", For[j = 1, j < n + 2, j++, w[m + 2, j] = w[m, j]], If[RightEdge == "S", For[j = 1, j < n + 2, j++, w[m + 2, j] = -w[m, j]], "Invalid Input Data!"]] w[0, 0] = w[2, 2]; w[0, n + 2] = w[2, n]; w[m + 2, 0] = w[m, 2]; w[m + 2, n + 2] = w[n, m];"

``


 

Download for.mw

 

Error, (in eval/diff) invalid input: diff received...

Asked by: 9009134 110
diff
February 04 2017
0  8

hi .in substuting in ode equation i encounter with error

Error, (in eval/diff) invalid input: diff received (f1[i+1, i]-2*f1[i, i]+f1[i-1, i])/h^2, which is not valid for its 2nd argument

please help me

thanksh_2.5_eq4.mw
 

 

restart; with(DEtools); with(plots); interface(rtablesize = 25); ode1 := -P1*(diff(f1(x, y), x, x, x, x, x, x)+diff(f1(x, y), y, y, y, y, y, y)+3*(diff(f1(x, y), x, x, x, x, y, y))+3*(diff(f1(x, y), x, x, y, y, y, y)))+P2*(diff(f1(x, y), x, x, x, x)+diff(f1(x, y), y, y, y, y)+2*(diff(f1(x, y), x, x, y, y)))-N_x*(diff(f1(x, y), x, x))-N_y*(diff(f1(x, y), y, y)); F1xx := i, proc (j) options operator, arrow; (f1[i+1, j]-2*f1[i, j]+f1[i-1, j])/h^2 end proc; F1yy := i, proc (j) options operator, arrow; (f1[i, j+1]-2*f1[i, j]+f1[i, j-1])/hy^2 end proc; F1xxxx := i, proc (j) options operator, arrow; (f1[i+2, j]-4*f1[i+1, j]+6*f1[i, j]-4*f1[i-1, j]+f1[i-2, j])/h^4 end proc; F1yyyy := i, proc (j) options operator, arrow; (f1[i, j+2]-4*f1[i, j+1]+6*f1[i, j]-4*f1[i, j-1]+f1[i, j-2])/hy^4 end proc; F1xxxxxx := i, proc (j) options operator, arrow; (f1[i+3, j]-6*f1[i+2, j]+15*f1[i+1, j]-20*f1[i, j]+15*f1[i-1, j]-6*f1[i-2, j]+f1[i-3, j])/h^6 end proc; F1yyyyyy := i, proc (j) options operator, arrow; (f1[i, j+3]-6*f1[i, j+2]+15*f1[i, j+1]-20*f1[i, j]+15*f1[i, j-1]-6*f1[i, j-2]+f1[i, j-3])/hy^6 end proc; F1xxyyyy := i, proc (j) options operator, arrow; (f1[i-1, j-2]-4*f1[i-1, j-1]+6*f1[i-1, j]-4*f1[i-1, j+1]+f1[i-1, j+2]-2*f1[i, j-2]+8*f1[i, j-1]-12*f1[i, j]+8*f1[i, j+1]-2*f1[i, j+2]+f1[i+1, j-2]-4*f1[i+1, j-1]+6*f1[i+1, j]-4*f1[i+1, j+1]+f1[i+1, j+2])/(hy^4*h^2) end proc; F1yyxxxx := i, proc (j) options operator, arrow; (f1[i-2, j-1]-2*f1[i-2, j]+f1[i-2, j+1]-4*f1[i-1, j-1]+8*f1[i-1, j]-4*f1[i-1, j+1]+6*f1[i, j-1]-12*f1[i, j]+6*f1[i, j+1]-4*f1[i+1, j-1]+8*f1[i+1, j]-4*f1[i+1, j+1]+f1[i+2, j-1]-2*f1[i+2, j]+f1[i+2, j+1])/(h^4*hy^2) end proc; eq := simplify(eval(ode1, {diff(f1(x, y), x, x) = F1xx(i, j), diff(f1(x, y), y, y) = F1yy(i, j), diff(f1(x, y), x, x, x, x) = F1xxxx(i, j), diff(f1(x, y), y, y, y, y) = F1yyyy(i, j), diff(f1(x, y), x, x, x, x, x, x) = F1xxxxxx(i, j), diff(f1(x, y), x, x, x, x, y, y) = F1yyxxxx(i, j), diff(f1(x, y), x, x, y, y, y, y) = F1xxyyyy(i, j), diff(f1(x, y), y, y, y, y, y, y) = F1yyyyyy(i, j)}))

Error, (in eval/diff) invalid input: diff received (f1[i+1, i]-2*f1[i, i]+f1[i-1, i])/h^2, which is not valid for its 2nd argument

  

NULL

 

expand(%)

i

 (1)

simplify(4*h^4*f3[i]*f4[i]*%)

4*h^4*f3[i]*f4[i]*i

 (2)

diff(f1(x, y), x, x, y, y, y, y)

diff(diff(diff(diff(diff(diff(f1(x, y), x), x), y), y), y), y)

 (3)

``

``

NULL


 

Download h_2.5_eq4.mw

 

nonlinear fit gives the right fit, but wrong param...

Asked by: shesia 20
parameters
February 03 2017
1  5

hello im currently working on a project and im fitting some non-linear data with the follwing model:

model:= 1.*10^5*exp(-t*k-*t/B)

where B and k are the unkown parameters.

i have the following data:

X=<0,2,4,6,8,10>

Y := <100000, 86089.76983, 74114.4849, 63804.98946, 54929.56828>

if i do the fit like this: fit(model,X,Y,t) i get a fit that looks like this:

100000*exp(-0.520664970792415e-1*t) which i alright. however when i try to get the parameters i get something like this

fit(model,X,Y,t,output=[parametervalues])

[B = .999563650781966, k = -.948370042622550] - when i define k=0.05 however i get [B = 483.910474528817], which is the right parameter.

insereting the two different sets yields:

1) k=0.05, B = 483.910474528817 -> 100000*exp(-0.520664970792415e-1*t) (good aprox to the fit)

2) B = .999563650781966 and k = -.948370042622550 -> 100000*exp(-0.5206649794e-1*S) (way off)

how is this possible anyone able to help, been sitting with this for days now?

 

 

Maple 2016 not starting...

Asked by: FadiBunni 10
February 02 2017
2  5

Hallo,
I got problem with starting maple 2016, when I click on the app files, it thinks for a second then nothing happens, I have tried to uninstall and install again, but no defference. I am not getting any error at any point, that is why it is hard to fix. maybe it has to do with java, jre or PATH?

Thanks for your help.  

Corrupt Maple document (I know that this exact pro...

Asked by: andreasamdisen 10
worksheet corruption
February 02 2017
2  9

https://www.dropbox.com/s/zpj7kudgn3gikfa/INS%20-%20Bygninger%20energibehov.mw?dl=0

https://www.dropbox.com/s/vj0t0p3r9stoy3c/INS%20-%20F%C3%B8ringsveje%20og%20indeklima.mw?dl=0

"There were problems during the loading process. Your worksheet may be incomplete."

I have seen threads where this has been solved, but i do not know how to do it.

It happened to both files, after restarting my computer.

Are anyone willing to help me?

complex tensor expression with Pauli matrices...

Asked by: illuminates 65
physics tensor
February 02 2017
0  0

I attempt to understand is it possible whether to calculate the next expression by means of mathematica:

where

{,}- anti-communicator;

 a=0,1,2,3;

\tau^{0}-unity matrix; \tau^{i} - Pauli matrix;

and:
 

 

 


 
finally - Levi-Civita symbol; -Hermitian conjugation.

 

Thank you for your kind replies .

First 32 33 34 35 36 37 38 Last Page 34 of 60