1function plotresultsPEBBL(nx,ny)
2
3control = load([int2str(nx),'x',int2str(ny),'.sol.txt']).';
4X       = load(['X_',int2str(nx),'_',int2str(ny),'.txt']);
5Y       = load(['Y_',int2str(nx),'_',int2str(ny),'.txt']);
6tmp = control;
7for i = 1:nx
8  for j = 1:ny
9    control(j+(i-1)*ny) = tmp(i+(j-1)*nx);
10  end
11end
12figure,
13scatter3(2*X/nx,Y/ny,control,100,control,'filled')
14view(2)
15axis equal;
16axis tight;
17box on;
18xlabel('x');
19ylabel('y');
20title('Computed controls');
21colorbar
22print('-depsc2',['control_',int2str(nx),'x',int2str(ny),'.eps']);
23
24end
25