Вот тот, который возвращает его в матрице 4x4, используя формулу из справочной страницы OpenGL gluPerspective :
static void my_PerspectiveFOV(double fov, double aspect, double near, double far, double* mret) {
double D2R = M_PI / 180.0;
double yScale = 1.0 / tan(D2R * fov / 2);
double xScale = yScale / aspect;
double nearmfar = near - far;
double m[] = {
xScale, 0, 0, 0,
0, yScale, 0, 0,
0, 0, (far + near) / nearmfar, -1,
0, 0, 2*far*near / nearmfar, 0
};
memcpy(mret, m, sizeof(double)*16);
}