simpler code: Q computation following Sarkka (removed abasset confusing matrix)