cw_vr.m

function phi_vr = cw_vr(t,n)
    
    phi_vr = [3*n*sin(n*t), 0, 0
                6*n*(cos(n*t)-1), 0, 0
                0, 0, -n*sin(n*t)];
    
end