function f = myobj(x)
f = 6 * x(1) + 4 * x(2) + 5 * x(3) + 3 * x(4);

function [g,ceq] = myconst(x)
pR1 = 0.85; pR2 = 0.75; pR3 = 0.80; pR4 = 0.70;
pF1 = 1 - pR1; pF2 = 1 - pR2; pF3 = 1 - pR3; pF4 = 1 - pR4;
g = 0.99 - (1 - pF1 ^ x(1)) * (1 - pF2 ^ x(2)) * (1 - pF3 ^ x(3)) * (1 - pF4 ^ x(4));
ceq = [];

clear
x0 = [1;1;1;1];
a = [-1 0 0 0;0 -1 0 0;0 0 -1 0;0 0 0 -1];
b = [0 0 0 0];
lb = [1;1;1;1];ub = [19;19;19;19];
[x,fval] = fmincon(@myobj,x0,a,b,[],[],lb,ub,@myconst);

例如这种是直接在一个julia脚本里就能实现么