From what I can figure out,bearing in mind it is late ^^ try
vectorw.f =cy*cp*cr - sy*sp*sr vectorx.f =cy*cp*sr + sy*sp*cr vectory.f =sy*cp*cr + cy*sp*sr vectorz.f =cy*sp*cr - sy*cp*sr