The only difference of this part with partII is here I extended the original lens computation module in 2D case to 3D case so that 3D rays can be traced into the lens system and compute outgoing rays. The results are similar to partII. The new code is here.