clc beta1=0.3; beta2=0.45; epsilon=1e-10; n=4; ntilde=n*(n+1)/2; X0=10*[1,0,0,0;0,1,0,0;0,0,1,-0.5;0,0,-0.5,1]; Y0=10*eye(n); A1=[1,0,0,0;0,1,0,0;0,0,0,0;0,0,0,0]; A2=[0,0,0,0;0,0,1,0;0,1,0,0;0,0,0,0]; A3=[0,0,0,0;0,0,1,0;0,1,1,0;0,0,0,0]; A4=[0,0,0,0;0,0,0,0;0,0,1,0;0,0,0,-1]; A5=[0,0,0,0;0,0,0,0;0,0,0,1;0,0,1,1]; A=[svec(A1)';svec(A2)';svec(A3)';svec(A4)';svec(A5)';zeros(ntilde/2,ntilde)]; B1=[1,0,0,0;0,-1,0,0;0,0,0,0;0,0,0,0]; B2=[0,1,0,0;1,0,0,0;0,0,0,0;0,0,0,0]; B3=[0,0,1,0;0,0,0,0;1,0,0,0;0,0,0,0]; B4=[0,0,0,1;0,0,0,0;0,0,0,0;1,0,0,0]; B5=[0,0,0,0;0,0,0,1;0,0,0,0;0,1,0,0]; B=[zeros(ntilde/2,ntilde);svec(B1)';svec(B2)';svec(B3)';svec(B4)';svec(B5)']; q=[1;zeros(9,1)]; tic; x=NTAlgorithm2(A,B,q,X0,Y0,beta1,beta2,epsilon); elapsedtime=toc; ratio=x(2:1+x(1)); file=fopen('testprogramoutput','w'); fprintf(file,'No of iterations = %d\n',x(1)); fprintf(file,'Elapsed Time = %f\n',elapsedtime); fprintf(file,'Ratio of duality gaps with iteration:\n'); fprintf(file,'%f\n',ratio); fclose(file); noofiterations=x(1) X=smat(x(x(1)+2:ntilde+x(1)+1)) Y=smat(x(ntilde+x(1)+2:2*ntilde+x(1)+1))