Advanced Search

Yappa

5 Reputation

One Badge

7 years, 133 days
Contact Yappa

MaplePrimes Activity


These are questions asked by Yappa

How do i plot an array of unit scaled po...

Yappa 5 Maple
plot units plotting
February 25 2018
1 0

Hi!

If i do something like this:

my_points := [[1, 2], [2, .9]];
plot(my_points, style = point)

It works ok. But then i try to take units in place:

my_scaled_points := [[1*Unit('kg'), 2*Unit('hour')], [21*Unit('kg'), 7*Unit('hour')]];
plot(my_scaled_points, style = point, useunits = [Unit('kg'), Unit('hour')])

..i got empty graph. How it can be fixed?

Thank you!

 

How do i can simplify polynomial with un...

Yappa 5 Maple 2017
simplify units
February 23 2018
1 3
 

Download test.mw
Udp: document attached.

> 

Q__gr := Typesetting[delayDotProduct](Vector[row](8, {(1) = 0, (2) = 5, (3) = 10, (4) = 15, (5) = 20, (6) = 25, (7) = 30, (8) = 35}), Unit('m'^3/'day'), true)

Vector[row](%id = 18446746861178193550)

 (1)
> 

`η__gr` := Vector[row](8, {(1) = 0, (2) = 9.324, (3) = 17, (4) = 23.232, (5) = 27.192, (6) = 29.6, (7) = 29.516, (8) = 24.592})

Vector[row](%id = 18446746861214593262)

 (2)
> 

`pointsη` := [seq([Q__gr[i], `η__gr`[i]], i = 1 .. 8)]

[[0, 0], [5*Units:-Unit(m^3/d), 9.324], [10*Units:-Unit(m^3/d), 17], [15*Units:-Unit(m^3/d), 23.232], [20*Units:-Unit(m^3/d), 27.192], [25*Units:-Unit(m^3/d), 29.6], [30*Units:-Unit(m^3/d), 29.516], [35*Units:-Unit(m^3/d), 24.592]]

 (3)
> 

CurveFitting[PolynomialInterpolation](`pointsη`, q)

0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d)

 (4)
> 

eta := proc (q) options operator, arrow; 0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d) end proc

proc (q) options operator, arrow; 0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d) end proc

 (5)
> 

eta(12*Unit('m'^3/'day'))

.8415811058*Units:-Unit(m^3/d)^7/Units:-Unit(m^3/d)^7-8.900090268*Units:-Unit(m^3/d)^6/Units:-Unit(m^3/d)^6+36.88243199*Units:-Unit(m^3/d)^5/Units:-Unit(m^3/d)^5-75.97578241*Units:-Unit(m^3/d)^4/Units:-Unit(m^3/d)^4+80.63761922*Units:-Unit(m^3/d)^3/Units:-Unit(m^3/d)^3-45.58150401*Units:-Unit(m^3/d)^2/Units:-Unit(m^3/d)^2+31.86194286*Units:-Unit(m^3/d)/Units:-Unit(m^3/d)

 (6)
> 

simplify(19.76619849)

Error, (in Units:-Standard:-+) the units `m^21/s^7` and `m^18/s^6` have incompatible dimensions

  
> 

``


 

Download test.mw

 

Hi!

I got a function (from CurveFitting), that produce a polynomial with the some units inside:

 

As i can see, it can be easy simplified, but i get:


 

Is where any trick to do it, without stripping units with convert(unit_free) or something?

Thank you!

 

Page 1 of 1