Question: Integration with vectors and units

restart

with(VectorCalculus); SetCoordinates('cartesian[x, y, z]')

with(ScientificConstants)

"with(Units[Simple]):"

``

g := evalf(Constant(g, units))

9.80665*Units:-Unit(m/s^2)

(1)

`#mover(mi("\`v__0\`"),mo("&rarr;"))` := `<,>`(v[0]*Unit('m'/'s'), 0*Unit('m'/'s'), 0*Unit('m'/'s'))

Vector(3, {(1) = v[0]*Units:-Unit(m/s), (2) = 0, (3) = 0})

(2)

`#mover(mi("a"),mo("&rarr;"))` := `<,>`(0*Unit('m'/'s'^2), g, 0*Unit('m'/'s'^2))

Vector(3, {(1) = 0, (2) = 9.80665*Units:-Unit(m/s^2), (3) = 0})

(3)

int(`#mover(mi("a"),mo("&rarr;"))`, t)

Vector(3, {(1) = 0, (2) = 9.80665*t*Units:-Unit(m/s^2), (3) = 0})

(4)

As you can see in (4) units are wrong. I have tried to insert them in int expression but nothing has work. Moreover using Units[Standard] packecge i get

Error, (in int) wrong number (or type) of arguments: wrong type of integrand passed to indefinite integration.

Thanks

Download Units_and_Integrals.mw

Please Wait...