

function out = wdot(t,w)

    out = [u(w(1), w(2)); v(w(1), w(2))]; 
