PJ = [3*2; 5*3; 5*3; 9*5];
Gx = [...
      2.8750,0.8484;...
      3.6025,0.8057;...
      3.3209,0.3283;...
      3.1712,0.8761;...
      ];
Qxe= [...
      42.3161,6.4356;...
      23.9335,5.4995;...
      22.9832,2.3751;...
      24.4796,8.1372;...
      ];
Qxi= [...
      54.7350,12.8821;...
      64.2177,15.8709;...
      60.5471,7.2580;...
      38.5331,10.8445;...
      ];



figure;
errorbar(PJ,Qxe(:,1),Qxe(:,2)); hold on;
errorbar(PJ,Qxi(:,1),Qxi(:,2));
errorbar(PJ,Gx(:,1),Gx(:,2));
set(gca,'Xscale','log')
legend({'e heat flux';'i heat flux';'particle flux'})
xlim([5 100])
xticks([6 15 45]);
xticklabels({'(2,1)','(4,2)','(8,4)'});