wieviele Schritte sollen höchstens ausgeführt werden?
numsteps = wichtig: numew <= numstep <= 300

wählen Sie die Reorthogonalisierungstechnik

	keine Reorthogonalisierung
	selektive Orthogonalisierung
	vollständige Reorthogonalisierung

wählen sie epsaccept für das Parlett-Kriterium
epsaccept = wichtig: 1.0e-16 <= epsaccept <= 1.0e-2

geben Sie hier den Startvektor an: n Einträge. die Einträge bitte mit Kommata oder Leerzeichen trennen. Sie können mehrere Zeilen benutzen.
Falls Sie keinen eigenen Vektor eingeben wollen , tragen Sie hier nur "xxx" ein. Dann wird intern ein Zufallsvektor generiert.

x0:

geben Sie hier das Programmstück ein, das die Operation y=A*x ausführt. das Resultat muss y heissen, der Eingangsvektor muss x heissen, die Dimension ist n . Ansonsten können Sie benutzen: eine Matrix A, deren Einträge Sie dann definieren müssten, einen Hilfsvektor zeisum mit n Komponenten, die Laufvariablen i,j,k, die logischen Variablen bool1, bool2, bool3, die Konstanten pi, e1=exp(1), sqrt2=sqrt(2.0) . Eigene Variablen können Sie nicht definieren. c quindiagonalmatrix mit den eigenwerten c 16sin(i*pi/(2(n+1)))**4, i=1,...,n y(1)=5.d0*x(1)-4.d0*x(2)+x(3) y(2)=-4.d0*x(1)+6.d0*x(2)-4.d0*x(3)+x(4) y(n-1)=-4.d0*x(n)+6.d0*x(n-1)-4.d0*x(n-2)+x(n-3) y(n)=5.d0*x(n)-4.d0*x(n-1)+x(n-2) do i=3,n-2 y(i)=6.d0*x(i)-(x(i-1)+x(i+1))*4.d0+x(i-2)+x(i+2) enddo

Klicken Sie auf "Auswerten", um die Aufgabe an das Programm zu schicken.