Question: Evalb and verify

Hi,

I am trying to proof that the right hand side is equal to the left hand side. When i do step by step, it clearly shows that my "c" and "d" are exactly the same thing. But when i use "evalb" and "verify", it comes out with "false". May I kwow what went wrong please?

Thanks.

NULL

NULLNULL

restart

with(LinearAlgebra):

``NULL

u(xi):

a := diff(u(xi), `$`(xi, 2)):

zip(`*`, E2, u(xi)) := E0*xi^2*a+E0*xi*zip(diff, u(xi), xi):

Q := E0*xi*(diff(u(xi), xi))+E1*u(xi):

X := Vector(2, {(1) = u(xi), (2) = Q})``

b := xi*zip(diff, X, xi):

c := Vector(2, {(1) = xi*(diff(u(xi), xi)), (2) = E2*u(xi)+E1*xi*(diff(u(xi), xi))})

Vector[column]([[xi*(diff(u(xi), xi))], [E2*u(xi)+E1*xi*(diff(u(xi), xi))]])

(1)

Z := Matrix(2, 2, {(1, 1) = E1/E0, (1, 2) = -1/E0, (2, 1) = -E2+E1^2/E0, (2, 2) = -E1/E0})``

d := simplify(-Typesetting:-delayDotProduct(Z, X))

Vector[column]([[xi*(diff(u(xi), xi))], [E2*u(xi)+E1*xi*(diff(u(xi), xi))]])

(2)

evalb(c = d);

false

(3)

zip(`-`, c, d);``

Vector[column]([[0], [0]])

(4)

verify(c, d, {'equal'})``

false

(5)

``

``



Download static8.1_part1.mw

Please Wait...