Intrinsic camera matrix from blender
The formula for alpha stated is incorrect.
If you set it to (focal length)*(image width)/(sensor width)
it will work.

Note that focal length and sensor width are in mm, so they can't both be in the numerator, just due to units issues.