Advanced Search

panke

145 Reputation

7 Badges

8 years, 355 days
Contact panke

MaplePrimes Activity


These are questions asked by panke

Way I get this Error it's work perfect o...

panke 145 Maple
programming linearalgebra
September 07 2019
2 2


 

> 

restart
> 

with(LinearAlgebra):
> 

G := Matrix(4, 4, {(1, 1) = 1, (1, 2) = 2, (1, 3) = 0, (1, 4) = 0, (2, 1) = 2, (2, 2) = -1, (2, 3) = 0, (2, 4) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = -1, (3, 4) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = -1})
> 

``
> 

TRJH(G)

Error, (in TRJH) invalid input: LinearAlgebra:-JordanForm expects value for keyword parameter output to be of type {list, identical(J), identical(Q)}, but received Q

  
> 

``


 

Download temp.mw

dsolve error ode system...

panke 145 Maple
ode dsolve numeric
April 14 2019
1 0


 

> 

restart
> 

odeg := diff(x(t), t)-p[x](t), diff(y(t), t)-p[y](t), diff(z(t), t)-p[z](t), diff(p[x](t), t)-5, diff(p[y](t), t)-3, diff(p[z](t), t)+10, diff(L[x](t), t)-4, diff(L[y](t), t)-20, diff(L[z](t), t)-8, diff(R[xx](t), t)-(-(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[x](t)-(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[y](t)-(R[zx](t)^2+(1/2)*R[zy](t)^2+R[zz](t)^2)*L[z](t))*R[yx](t)-((R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[x](t)+(R[yx](t)^2+(1/2)*R[yy](t)^2+R[yz](t)^2)*L[y](t)+(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[z](t))*R[zx](t), diff(R[yx](t), t)-((R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[x](t)+(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[y](t)+(R[zx](t)^2+(1/2)*R[zy](t)^2+R[zz](t)^2)*L[z](t))*R[xx](t)-(-(R[xx](t)^2+(1/2)*R[xy](t)^2+R[xz](t)^2)*L[x](t)-(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[y](t)-(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[z](t))*R[zx](t), diff(R[zx](t), t)-(-(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[x](t)-(R[yx](t)^2+(1/2)*R[yy](t)^2+R[yz](t)^2)*L[y](t)-(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[z](t))*R[xx](t)-((R[xx](t)^2+(1/2)*R[xy](t)^2+R[xz](t)^2)*L[x](t)+(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[y](t)+(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[z](t))*R[yx](t), diff(R[xy](t), t)-(-(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[x](t)-(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[y](t)-(R[zx](t)^2+(1/2)*R[zy](t)^2+R[zz](t)^2)*L[z](t))*R[yy](t)-((R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[x](t)+(R[yx](t)^2+(1/2)*R[yy](t)^2+R[yz](t)^2)*L[y](t)+(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[z](t))*R[zy](t), diff(R[yy](t), t)-((R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[x](t)+(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[y](t)+(R[zx](t)^2+(1/2)*R[zy](t)^2+R[zz](t)^2)*L[z](t))*R[xy](t)-(-(R[xx](t)^2+(1/2)*R[xy](t)^2+R[xz](t)^2)*L[x](t)-(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[y](t)-(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[z](t))*R[zy](t), diff(R[zy](t), t)-(-(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[x](t)-(R[yx](t)^2+(1/2)*R[yy](t)^2+R[yz](t)^2)*L[y](t)-(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[z](t))*R[xy](t)-((R[xx](t)^2+(1/2)*R[xy](t)^2+R[xz](t)^2)*L[x](t)+(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[y](t)+(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[z](t))*R[yy](t), diff(R[xz](t), t)-(-(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[x](t)-(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[y](t)-(R[zx](t)^2+(1/2)*R[zy](t)^2+R[zz](t)^2)*L[z](t))*R[yz](t)-((R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[x](t)+(R[yx](t)^2+(1/2)*R[yy](t)^2+R[yz](t)^2)*L[y](t)+(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[z](t))*R[zz](t), diff(R[yz](t), t)-((R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[x](t)+(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[y](t)+(R[zx](t)^2+(1/2)*R[zy](t)^2+R[zz](t)^2)*L[z](t))*R[xz](t)-(-(R[xx](t)^2+(1/2)*R[xy](t)^2+R[xz](t)^2)*L[x](t)-(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[y](t)-(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[z](t))*R[zz](t), diff(R[zz](t), t)-(-(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[x](t)-(R[yx](t)^2+(1/2)*R[yy](t)^2+R[yz](t)^2)*L[y](t)-(R[yx](t)*R[zx](t)+(1/2)*R[yy](t)*R[zy](t)+R[yz](t)*R[zz](t))*L[z](t))*R[xz](t)-((R[xx](t)^2+(1/2)*R[xy](t)^2+R[xz](t)^2)*L[x](t)+(R[xx](t)*R[yx](t)+(1/2)*R[xy](t)*R[yy](t)+R[xz](t)*R[yz](t))*L[y](t)+(R[xx](t)*R[zx](t)+(1/2)*R[xy](t)*R[zy](t)+R[xz](t)*R[zz](t))*L[z](t))*R[yz](t):
> 

Icg := x(0) = 0, y(0) = 0, z(0) = 0, p[x](0) = 0, p[y](0) = 0, p[z](0) = 0, L[x](0) = 0, L[y](0) = 0, L[z](0), R[xx](0) = 1, R[yx](0) = 0, R[zx](0) = 0, R[xy](0) = 0, R[yy](0) = 1, R[zy](0) = 0, R[xz](0) = 0, R[yz](0) = 0, R[zz](0) = 1:
> 

ld := dsolve([odeg, Icg], numeric, range = 0 .. 1)

Error, (in dsolve/numeric/make_proc) invalid input: `DEtools/convertsys` expects its 2nd argument, inits, to be of type {`=`, set(`=`), list(`=`)}, but received {L[z](0), x(0) = 0, y(0) = 0, z(0) = 0, L[x](0) = 0, L[y](0) = 0, R[xx](0) = 1, R[xy](0) = 0, R[xz](0) = 0, R[yx](0) = 0, R[yy](0) = 1, R[yz](0) = 0, R[zx](0) = 0, R[zy](0) = 0, R[zz](0) = 1, p[x](0) = 0, p[y](0) = 0, p[z](0) = 0}

  
> 

``


 

Download dsolve_error.mw

 

 

 

plot error (with proc )...

panke 145 Maple 13
plot
March 10 2019
2 1

Why this plot dosen't work ?


 

``

> 

restart; with(linalg); with(plots)
> 

``
> 

FS := proc (`&xi;h`, m0, `&xi;wp`, `&mu;v`) local f; if abs(xi) <= 1 then f := proc (xi) options operator, arrow; m0*xi*(1-abs(xi))^2+xi*abs(xi)*(3-2*abs(xi)) end proc else f := proc (xi) options operator, arrow; signum(xi)*((`&xi;wp`-1)^2+`&mu;v`*(abs(xi)-1)^2)/((`&xi;wp`-1)^2+(abs(xi)-1)^2) end proc end if end proc:
> 

``
> 

> 

plot(FS(xi, 1.5, 2, .8), xi = 0 .. 1.5)

Error, (in FS) cannot determine if this expression is true or false: abs(xi) <= 1

  
> 

> 

``


 

Download plottest.mw
 

``

> 

restart; with(linalg); with(plots)
> 

``
> 

FS := proc (`&xi;h`, m0, `&xi;wp`, `&mu;v`) local f; if abs(xi) <= 1 then f := proc (xi) options operator, arrow; m0*xi*(1-abs(xi))^2+xi*abs(xi)*(3-2*abs(xi)) end proc else f := proc (xi) options operator, arrow; signum(xi)*((`&xi;wp`-1)^2+`&mu;v`*(abs(xi)-1)^2)/((`&xi;wp`-1)^2+(abs(xi)-1)^2) end proc end if end proc:
> 

``
> 

> 

plot(FS(xi, 1.5, 2, .8), xi = 0 .. 1.5)

Error, (in FS) cannot determine if this expression is true or false: abs(xi) <= 1

  
> 

> 

``


 

Download plottest.mw

 

Error with Procedure output...

panke 145 Maple 13
matrix 2dmath proc
September 16 2018
3 0


 

> 

restart
> 

" S3:=proc(psi) [[[cos(psi),-sin(psi),0],[sin(psi),cos(psi),0],[0,0,1]]]: end proc ; S1:=proc(`&varphi;`) [[[1,0,0],[0,cos(`&varphi;`),-sin(`&varphi;`)],[0,sin(`&varphi;`),cos(`&varphi;`)]]]: end proc;"

proc (psi) Matrix(3, 3, {(1, 1) = cos(psi), (1, 2) = -sin(psi), (1, 3) = 0, (2, 1) = sin(psi), (2, 2) = cos(psi), (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1}) end proc

  

proc (`&varphi;`) Matrix(3, 3, {(1, 1) = 1, (1, 2) = 0, (1, 3) = 0, (2, 1) = 0, (2, 2) = cos(psi), (2, 3) = -sin(psi), (3, 1) = 0, (3, 2) = sin(psi), (3, 3) = cos(psi)}) end proc

 (1)
> 

> 

``
> 

S3(uu)

Matrix([[cos(uu), -sin(uu), 0], [sin(uu), cos(uu), 0], [0, 0, 1]])

 (2)
> 

S1(uu)

Matrix([[1, 0, 0], [0, cos(psi), -sin(psi)], [0, sin(psi), cos(psi)]])

 (3)
> 

``


 

Download test3.mw

Way i get the wrong output with procedure S1(uu)
 

Arrow plot color error...

panke 145 Maple 13
arrow plotting
September 03 2018
6 3

arrow([u1, u2, u3], shape = cylindrical_arrow, length = sqrt(LinearAlgebra:-Transpose(f0).f0), color = [red, black, green]):

Why i get this error

Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.Error: the number of colours in the COLOUR structure is not sufficient.

3 4 5 6 7 8 9 Page 5 of 9