m=input('tridiag test eingabe m= ');
xlow=input('eingabe xlow= ');
xup=input('eingabe xup= ');
a=input('eingabe faktor a= ');
x=linspace(xlow,xup,m);
 h=(xup-xlow)/(m-1);
 A=diag(-2*ones(m,1))+diag(ones(m-1,1),-1)+diag(ones(m-1,1),1) ...
+h*diag(1./x)*( diag(-ones(m,1))+diag(ones(m-1,1),1) ) + ...
 h^2*a*diag(1./x.^2)*diag(ones(m,1));
eig(A)