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