clf

Ngrid = 35;

sigma1 = 3.0;
sigma2 = 1.25;
omega  = 0.05;

p = [sigma1; sigma2; omega];

v1grid = linspace(-1.5,1.5,Ngrid);
v2grid = linspace(-1.5,1.5,Ngrid);
[V1,V2] = meshgrid(v1grid,v2grid);
Q2 = cpldosc_fhat_q2(V1,V2,p);
h = mesh(V1,V2,Q2);
set(h,'EdgeColor','k')

xlabel('v_1')
ylabel('v_2')
zlabel('q_2')

hold on

%
% Load data from the file 'slow.dat'
%
x = load('slow.dat','-ascii');
v1 = x(:,2);
v2 = x(:,3);
q1 = x(:,4);
q2 = x(:,5);
plot3(v1,v2,q2,'b','linewidth',3)

%
% Load data from the file 'fast.dat'
%
x = load('fast.dat','-ascii');
v1 = x(:,2);
v2 = x(:,3);
q1 = x(:,4);
q2 = x(:,5);
plot3(v1,v2,q2,'r','linewidth',3)

view(125,25)
alpha(.1)