P=[0 1 -6 -4 5;0 2 1 0 0]
M=bezier(P,1000)

scf(0)
clf()
plot2d(M(1,:),M(2,:))
plot2d(P(1,:),P(2,:),-3)

M=bezierOP(P,1000)

scf(1)
clf()
plot2d(M(1,:),M(2,:))
plot2d(P(1,:),P(2,:),-3)
