how to plot 3d pareto front for multiobjective linear programming problem

clear
clc
close all
fun = @(x) [84718*x(1)+72627.42*x(2)+41021.98*x(3)+58077.44*x(4),...
    1124.424*x(1)+1004.247*x(2)*8.40858*x(3)+10.52292*x(4),...
    51.524*x(1)+47.666*x(2)+19.136*x(3)+19.134*x(4)];
%options = optimoptions('gamultiobj',"PlotFcn","gaplotpareto","PopulationSize",200);
options = optimoptions('paretosearch','UseVectorized',true,'ParetoSetSize',200,...
    'plotfcn','psplotparetof');
Aineq = [1 1 1 1; 220.086 210.356 92.672 135.994; 464.75 203.194 261.548 206.638];
Bineq = [49.9; 660; 1136];
lb = [0.135277; 0.075732; 0; 0];
ub = [];
rng default
[x,f] = paretosearch(fun,4,[],[],lb,ub,[],options);
%opts.PopulationSiz = 500
%[xg,fg] = gamultiobj(fun,4,[],[],Aineq, Bineq,lb,ub,[],options);
view(-47,50)
disp(xg(4,:))
%what is wrong in this code kindly explain it