Glu-less пример рисования куба?

Предлагаю взглянуть на http://openglbook.com/ Там написано, что это OpenGL 4, но всё описано здесь ( пока только основы) напрямую относится и к OpenGL 3. Я думаю, что это лучший ресурс для быстрого начала работы с последними версиями OpenGL. Другим очень хорошим ресурсом является Learning Modern 3D Graphics Programming, который чуть более углубленно.

Вы должны заметить, что OpenGL 3 и новее игнорирует фиксированный конвейер. Вы все еще можете использовать фиксированный конвейер, используя профиль совместимости, но IMO лучше изучить программируемый конвейер. Это означает, что больше не будет непосредственного режима (например, glVertex3f и подобного) и списков отображения. Вы должны отправить свои данные вершин непосредственно на графическую карту перед рисованием. Это делается с помощью буферных объектов, которые представляют собой не что иное, как массив данных вершин. Затем вы должны отправить программу GLSL, которая будет обрабатывать эти данные, например, применять модели и проекционные матрицы (даже если вам не нужны причудливые шейдеры).

И последнее, но не менее важное: текущие версии OpenGL не реализуют матричные операции и освещение (опять же, если не используется профиль совместимости). Вам нужно создать все матрицы самостоятельно, но так как это утомительная работа, я предлагаю использовать внешнюю библиотеку.

Все это может звучать пугающе, и, честно говоря, мне было труднее изучить программируемый конвейер, чем фиксированный. Однако это дает вам гораздо больше гибкости.

В OpenTK есть свои математические классы, включая матрицу 4x4 — Matrix4. Он включает в себя как метод ортогональной проекции, так и метод перспективной проекции - Matrix4.CreateOrthographic(...); и Matrix4.CreatePerspectiveFieldOfView(...);. Также есть метод lookAt: Matrix4.LookAt(...);

Привязки OpenTK для OpenGL принимают Matrix4 для таких методов, как GL.LoadMatrix(), GL.UniformMatrix4() и т. д.

Из Ch4 книги OpenGL GLu.Perspective больше или менее эквивалентно:

98  Matrix CreateProjectionMatrix(
99      float fovy,
100     float aspect_ratio,
101     float near_plane,
102     float far_plane
103 )
104 {
105     Matrix out = { { 0 } };
106  
107     const float
108         y_scale = Cotangent(DegreesToRadians(fovy / 2)),
109         x_scale = y_scale / aspect_ratio,
110         frustum_length = far_plane - near_plane;
111  
112     out.m[0] = x_scale;
113     out.m[5] = y_scale;
114     out.m[10] = -((far_plane + near_plane) / frustum_length);
115     out.m[11] = -1;
116     out.m[14] = -((2 * near_plane * far_plane) / frustum_length);
117  
118     return out;
119 }