>> c=eye(5) c = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 >> for kk=3:5 c(kk,kk-1)=kk/150.0 end ??? ,kk-1)=kk/150.0 end | Missing operator, comma, or semi-colon. >> for kk=3:5 c(kk,kk-1)=kk/150.0 end c = 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0.0200 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 c = 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0.0200 1.0000 0 0 0 0 0.0267 1.0000 0 0 0 0 0 1.0000 c = 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0.0200 1.0000 0 0 0 0 0.0267 1.0000 0 0 0 0 0.0333 1.0000 >> for kk=3:5 c(kk-1,kk))=kk/150.0 ??? kk=3:5 c(kk-1,kk)) | Missing operator, comma, or semi-colon. >> for kk=3:5 c(kk-1,kk)=kk/150.0 end c = 1.0000 0 0 0 0 0 1.0000 0.0200 0 0 0 0.0200 1.0000 0 0 0 0 0.0267 1.0000 0 0 0 0 0.0333 1.0000 c = 1.0000 0 0 0 0 0 1.0000 0.0200 0 0 0 0.0200 1.0000 0.0267 0 0 0 0.0267 1.0000 0 0 0 0 0.0333 1.0000 c = 1.0000 0 0 0 0 0 1.0000 0.0200 0 0 0 0.0200 1.0000 0.0267 0 0 0 0.0267 1.0000 0.0333 0 0 0 0.0333 1.0000 >> >> c c = 1.0000 0 0 0 0 0 1.0000 0.0200 0 0 0 0.0200 1.0000 0.0267 0 0 0 0.0267 1.0000 0.0333 0 0 0 0.0333 1.0000 >> det(c) ans = 0.9978 >> [vec, val]=eig(c) vec = 1.0000 0 0 0 0 0 -0.6708 0.6708 0.2236 0.2236 0 0.5000 0.5000 -0.5000 0.5000 0 0.2236 -0.2236 0.6708 0.6708 0 -0.5000 -0.5000 -0.5000 0.5000 val = 1.0000 0 0 0 0 0 0.9851 0 0 0 0 0 1.0149 0 0 0 0 0 0.9553 0 0 0 0 0 1.0447 >> uff=sqrt(val) uff = 1.0000 0 0 0 0 0 0.9925 0 0 0 0 0 1.0074 0 0 0 0 0 0.9774 0 0 0 0 0 1.0221 >> uff*uff ans = 1.0000 0 0 0 0 0 0.9851 0 0 0 0 0 1.0149 0 0 0 0 0 0.9553 0 0 0 0 0 1.0447 >> y ??? Undefined function or variable y. >> const y ??? const y | Missing operator, comma, or semi-colon. >> c1=[ 1 1 1 1 1 ] c1 = 1 1 1 1 1 >> >> c1=0.01*c1 c1 = 0.0100 0.0100 0.0100 0.0100 0.0100 >> c2=c1; c3=c1; c4=c1;c5=c1; >> c5 c5 = 0.0100 0.0100 0.0100 0.0100 0.0100 >> c1(1)= 1.0; >> c2(2)= 1.0; >> c3(3)= 1.0; >> c4(4)= 1.0; >> c5(5)= 1.0; >> >> c4 c4 = 0.0100 0.0100 0.0100 1.0000 0.0100 >> c4'*c4 ans = 1.0004 >> c= [c1'*c1 c1'*c2 c1'*c3 c1'*c4 c1'*c5 c2'*c1 c2'*c2 c2'*c3 c2'*c4 c2'*c5 c3'*1 c3'*c2 c3'*c3 c3'*c4 c3'*c5 c4'*c1 c4'*c2 c4'*c3 c4'*c4 c4'*c5 c5'*c1 c5'*c2 c5'*c3 c5'*c4 c5'*c5 ] ??? All rows in the bracketed expression must have the same number of columns. >> >> [c2'*c1 c2'*c2 c2'*c3 c2'*c4 c2'*c5] ans = 0.0203 1.0004 0.0203 0.0203 0.0203 >> c= [c1'*c1 c1'*c2 c1'*c3 c1'*c4 c1'*c5 c2'*c1 c2'*c2 c2'*c3 c2'*c4 c2'*c5 c3'*c1 c3'*c2 c3'*c3 c3'*c4 c3'*c5 c4'*c1 c4'*c2 c4'*c3 c4'*c4 c4'*c5 c5'*c1 c5'*c2 c5'*c3 c5'*c4 c5'*c5 ]; >> >> det(c) ans = 0.9980 >> c c = 1.0004 0.0203 0.0203 0.0203 0.0203 0.0203 1.0004 0.0203 0.0203 0.0203 0.0203 0.0203 1.0004 0.0203 0.0203 0.0203 0.0203 0.0203 1.0004 0.0203 0.0203 0.0203 0.0203 0.0203 1.0004 >> v=[c1 c2 c3 c4 c5] v = 1.0000 0.0100 0.0100 0.0100 0.0100 0.0100 1.0000 0.0100 0.0100 0.0100 0.0100 0.0100 1.0000 0.0100 0.0100 0.0100 0.0100 0.0100 1.0000 0.0100 0.0100 0.0100 0.0100 0.0100 1.0000 >> v(1) ans = 1 >> v(1,1:i) ans = [] >> v(1,1:5) ans = 1.0000 0.0100 0.0100 0.0100 0.0100 >> v(;,5) ??? v(; | Missing variable or function. >> v(:,5) ans = 0.0100 0.0100 0.0100 0.0100 1.0000 >> c5=(v:,5) ??? c5=(v: | A closing right parenthesis is missing. Check for a missing ")" or a missing operator. >> c5 c5 = 0.0100 0.0100 0.0100 0.0100 1.0000 >> v(:,5) ans = 0.0100 0.0100 0.0100 0.0100 1.0000 >> stran=[ 0 0 0 0 0 ] stran = 0 0 0 0 0 >> stran=zeros(5) stran = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 >> stran(1,1)=1.0 stran = 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 >> ny=eye(4); >> s=[0.1;0.1;0.1;0.1;] s = 0.1000 0.1000 0.1000 0.1000 >> h=[ny s] h = 1.0000 0 0 0 0.1000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 >> s=[0.1 0.1 0.1 0.1 1.0]; >> h=[h;s] h = 1.0000 0 0 0 0.1000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 0.1000 0.1000 0.1000 0.1000 1.0000 >> det(h) ans = 0.9600 >> h(5,1)=0.2 h = 1.0000 0 0 0 0.1000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 0.2000 0.1000 0.1000 0.1000 1.0000 >> h(1,5)=0.2 h = 1.0000 0 0 0 0.2000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 0.2000 0.1000 0.1000 0.1000 1.0000 >> det(h) ans = 0.9300 >> q=1.0-h(5,:)*h(:,5) q = -0.0700 >> q=1.0-h(5,1:4)*h(1:4,5) q = 0.9300 >> q=1.0-h(5,1:4)*h(1:4,5) q = 0.9300 >> >> q = ??? q = | Missing variable or function. >> >> 0.9300 ans = 0.9300 >> >> >> ??? > | Missing variable or function. >> ny=eye(4); >> s=[0.1;0.1;0.1;0.1;] s = 0.1000 0.1000 0.1000 0.1000 >> h=[ny s] h = 1.0000 0 0 0 0.1000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 >> s=[0.1 0.1 0.1 0.1 1.0]; >> h=[h;s] h = 1.0000 0 0 0 0.1000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 0.1000 0.1000 0.1000 0.1000 1.0000 >> det(h) ans = 0.9600 >> h(5,1)=0.2 h = 1.0000 0 0 0 0.1000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 0.2000 0.1000 0.1000 0.1000 1.0000 >> h(1,5)=0.2 h = 1.0000 0 0 0 0.2000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 0.2000 0.1000 0.1000 0.1000 1.0000 >> >> det(h) ans = 0.9300 >> q=1.0-h(5,1:4)*h(1:4,5) q = 0.9300 >> ny=eye(4); >> s=[0.1;0.1;0.1;0.1;]; >> h=[ny s]; >> s=[0.1 0.1 0.1 0.1 1.0]; >> h=[h;s] h = 1.0000 0 0 0 0.1000 0 1.0000 0 0 0.1000 0 0 1.0000 0 0.1000 0 0 0 1.0000 0.1000 0.1000 0.1000 0.1000 0.1000 1.0000 >> q=1.0-h(5,1:4)*h(1:4,5) q = 0.9600 >> [lvec,lval]=eig(h) lvec = 0.8333 -0.1667 -0.1667 0.3536 0.3536 -0.1667 0.8333 -0.1667 0.3536 0.3536 -0.1667 -0.1667 0.8333 0.3536 0.3536 -0.5000 -0.5000 -0.5000 0.3536 0.3536 0 0 0 0.7071 -0.7071 lval = 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 0 0 0 0 0 1.2000 0 0 0 0 0 0.8000 >> det(lval) ans = 0.9600 >> quit 2093 flops.