!rm ho
diary ho
chdir ../chap3
n=7;a=rand(n,n);t=triu(a);b=rand(n,1);
y=backsub(t,b);y1=backsubcopy(t,b);y-y1
chdir ../chap4
a=rand(n,n);b=rand(n,1);
x=gauss(a,b);x1=gauss(a,b);x-x1
a1=a; a=housmulp(a,b);a1=housmulpcopy(a1,b);a-a1
a=rand(n,n);t=triu(a);
t1=t; t=invuptr(t);t1=invuptrcopy(t1);t-t1
a=rand(n,n);a1=a;
a=phousmul(a,b);a1=phousmulcopy(a1,b);a-a1
chdir ../chap5
a=rand(n,n);a1=a;
[a,u,m,q]=compiv(a);[a1,u1,m1,q1]=compivcopy(a1);a-a1,u-u1,m-m1,q-q1
a=rand(n,n);a1=a;
[h,p,a,v]=houshess(a);[h1,p1,a1,v1]=houshesscopy(a1);h-h1,p-p1,a-a1,v-v1
a=rand(n,n);a1=a;
[q,r,a,v]=housqr(a);[q1,r1,a1,v1]=housqrcopy(a1);q-q1,r-r1,a-a1,v-v1
a=rand(n,n);a1=a;
[q,r,a,v]=housqrn(a);[q1,r1,a1,v1]=housqrncopy(a1);q-q1,r-r1,a-a1,v-v1
a=rand(n,n);a1=a;
[a,l,u]=lugsel(a);[a1,l1,u1]=lugselcopy(a1);
a=rand(n,n);a1=a;
[a,u,m]=parpiv(a);[a1,u1,m1]=parpivcopy(a1);
chdir ../chap6
a=rand(n,n);a=a*a'+eye(n,n);a1=a;
[a,h]=choles(a);[a1,h1]=cholescopy(a1);a-a1,h-h1
a=rand(n,n); l=tril(a);b=rand(n,1);
y=forelm(l,b);y1=forelmcopy(l,b);y-y1
a=rand(n,n);x0=rand(n,1);b=rand(n,1);ep=0.0001;num=20;
[x,it]=gaused(a,x0,b,ep,num);[x1,it1]=gausedcopy(a,x0,b,ep,num);
x-x1,it-it1
a=rand(n,n);b=rand(b,1);
x=gausswf(a,b);x1=gausswfcopy(a,b);x-x1
a=rand(n,n);x0=rand(n,1);b=rand(n,1);ep=0.0001;num=20;
[x,it]=sucov(a,x0,b,w,ep,num);[x1,it1]=sucovcopy(a,x0,b,w,ep,num);
x-x1,it-it1
chdir ../chap7
a=rand(n,n);
[q,r]=clgrsch(a);[q1,r1]=clgrschcopy(a);q-q1,r-r1
a=rand(n,n);
[q,r]=mdgrsch(a);[q1,r1]=mdgrschcopy(a);q-q1,r-r1
a=rand(n,n);
x=varcovar(a);x1=varcovarcopy(a);x-x1
diary off