ROTATE MATRICE

vector N = 0;
vector scale = 1; //No difference in scale
vector rot = 0; //Additonal rot
vector pivot = point(2, "P", 0);
vector axis = normalize(point(1, "P", 0));
vector axis_new = cross(axis, {0, 1, 0}); // Get perpendicular angle for rotation

float degrees = radians(point(0, "time", @ptnum));

matrix3 m = ident();

rotate(m, degrees, axis_new);
vector4 orient = quaternion(m);
matrix m2 = instance(pivot, N, scale, rot, orient, pivot);

//Apply
@P *= m2;