Plotting a 3-axis pareto front

I am doing PSO for 3 objective functions with the code:

function PlotCosts(pop,rep)

    pop_costs=[pop.Cost];
    sc=size(pop_costs,2);
    pop_costs=reshape(pop_costs,3,sc/3);
    plot(pop_costs(1,:),pop_costs(2,:),'ko');
    hold on;

    rep_costs=[rep.Cost];
    sc=size(rep_costs,2);
    rep_costs=reshape(rep_costs,3,sc/3);
    plot(rep_costs(1,:),rep_costs(2,:),'r*');

    xlabel('Energetic efficiency (%)');
    ylabel('Life cycle cost (RM)');

    grid on;

    hold off;
end;

May I know how should i modify the above code to plot a 3-d pareto front?