Как найти максимум функции, у которой слагаемые F - это функции распределений случайной величины

clc
clear all
format short
 
T=10;
c_0=input("Введите значение c_0= ")%;
c_1=input("Введите значение c_1= ")%;
c_2=input("Введите значение c_2= ")%;
la=input("Введите значение la= ")%;
si=input("Введите значение si= ")%;
et=input("Введите значение et= ")%;
beta=input("Введите значение beta= ")%;
 
f = @(t)(1-(1-exp(-la*t)));
F = integral(f,0,100)
f_1 = @(t)(1-(1-exp(1)*(-t.^2/2*si.^2)));
F_1 = integral(f_1,0,100)
f_2 = @(t)(1-(1-exp(-(t/et).^beta)));
F_2 = integral(f_2,0,100)
u=0;  %минимальное значение
%#Задание Максимального значения функционала для сравнения на первом шаге#%'
Max_u=0;
MaxS=(c_0*F-c_1*(1-(u))*F_1-c_2*(u)*F_2)/(F+(1-(u))*F_1+(u)*F_2);
M=rand(101);
i=1;
un1=rand(1);
for u=0:0.1:T
    S=(c_0*F-c_1*(1-(u))*F_1-c_2*(u)*F_2)/(F+(1-(u))*F_1+(u)*F_2);
    M(i)=S;
    if (S>MaxS) %Нахождение максимального значения%
        MaxS=S;
        Max_u=u;
    endif
    i=i+1;
end
M;
printf("Максимум функционала = %d\n", MaxS);
printf("Значение u при максимуме функционала= %d\n",Max_u)
k=1;
for u=0:0.1:T
    
    un1(k)=u1;
    k=k+1;
end
 
 
un1;
surf(un1, M)