How to Find Pareto Optimal Solutions Using Matlab?

function YY = multi_objective_function(X)

for i = 1:2
    S1(i) = -10*exp(-0.2*sqrt(X(i).^2+X(i+1).^2));
end
F1 = sum(S1); % objective function 1

for i = 1:3
    S2(i) = abs(X(i)).^0.8 + 5*sin(X(i).^3);
end
F2 = sum(S2); % objective function 2

YY = [F1 F2];

function [x,fval,exitflag,output,population,score] = solving(nvars,lb,ub,PopulationSize_Data,MaxGenerations_Data)
%% This is an auto generated MATLAB file from Optimization Tool.

%% Start with the default options
options = optimoptions('gamultiobj');
%% Modify options setting
options = optimoptions(options,'PopulationSize', PopulationSize_Data);
options = optimoptions(options,'MaxGenerations', MaxGenerations_Data);
options = optimoptions(options,'CrossoverFcn', {  @crossoverintermediate [] });
options = optimoptions(options,'Display', 'off');
options = optimoptions(options,'PlotFcn', { @gaplotpareto });
[x,fval,exitflag,output,population,score] = ...
gamultiobj(@multi_objective_function,nvars,[],[],[],[],lb,ub,[],options);

clc
clear all
close all

nvars = 3;
lb = [-5 -5 -5];
ub = [5 5 5];
PopulationSize_Data = 200;
MaxGenerations_Data = 50;

[x,fval,exitflag,output,population,score] = solving(nvars,lb,ub,PopulationSize_Data,MaxGenerations_Data);

optimal_solution = x