mat4 QuaternionToMat4(vec4 quat)
{
	// normalize and create a local copy of the vector.
	vec3 axis = normalize(quat.xyz);

	// calculate angles
	float c = cos(-quat.w);
	float s = sin(-quat.w);
	float t = 1.0f - c;

	// do the conversion math once
	float tXX = t * axis.x * axis.x;
	float tXY = t * axis.x * axis.y;
	float tXZ = t * axis.x * axis.z;
	float tYY = t * axis.y * axis.y;
	float tYZ = t * axis.y * axis.z;
	float tZZ = t * axis.z * axis.z;
	float sinX = s * axis.x;
	float sinY = s * axis.y;

scaling matrix would look like:
x 0 0 0
0 y 0 0
0 0 z 0
0 0 0 1

GL.BindBuffer(BufferTarget.ArrayBuffer, mbo[material]);
GL.VertexAttribPointer(pos1, 4, VertexAttribPointerType.Float, false, stride, Vector4.SizeInBytes * 0);
GL.VertexAttribPointer(pos2, 4, VertexAttribPointerType.Float, false, stride, Vector4.SizeInBytes * 1);
GL.VertexAttribPointer(pos3, 4, VertexAttribPointerType.Float, false, stride, Vector4.SizeInBytes * 2);
GL.VertexAttribPointer(pos4, 4, VertexAttribPointerType.Float, false, stride, Vector4.SizeInBytes * 3);
GL.EnableVertexAttribArray(pos1);