求助fmincon函数的问题

function f=trelay(x)
f=x(1)+x(2)+x(3);
function[c,ceq]=mycon(x,a)
R=1;N0=1;
g(R)=[(2^(3*R)-1)*N0]^3;
ceq=g(R)/x(1)*(1/(5*x(1))+1/(2*x(2)))*(1/(12*x(1))+1/(24*x(3))))-a;
c=[];
clc
clear all
a=[0:0.01:0.05];
A=[];%线性不等式约束左边矩阵
b=[];%线性不等式约束右边向量
Aeq=[];%线性等式约束左边矩阵
beq=[];%线性等式约束右边向量
lb=zeros(1,3);
ub=[35;35;35];%自变量上限
x0=[1;1;1];%初始值
options=optimset('LargeScale','off','display','iter');
[x,fval]=fmincon(@trelay,x0,[],[],[],[],lb,ub,@mycon,options);
plot(a,x(1),'gv-');
求解答,我将约束条件中加了个变量a,最后要画图，将a作为X轴。不知怎么做