% trisys, Hur effektiv är tridia
  n=25;
  disp('          n      op. tridia   op. A\b ')
  for sys=1:4
    dia=4*ones(n,1);
    sup=ones(n-1,1);  sub=sup;
    b=(1:n)';
    flops(0) ,x=tridia(dia,sup,sub,b);
    tridiaflops=flops;
    A=4*eye(n)+diag(sup,1)+diag(sub,-1);
    flops(0),  x=A\b;
    gaussflops=flops;
    disp([n tridiaflops gaussflops])
    n=2*n;
  end