We can use the following equation to move from 3D world to 2D camera coordinates:
And the camera center can be calculated using the following equation:
The results are shown in the following images:
We can solve for the Fundamental matrix using the following equation:
The results (drawing epilolar line) are shown below:
With the max iteration number set to 1000, the outlier ratio set to 0.99, and the error threshold set to 0.005, using the RANSAC algorithm, we get the following results.