Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clear
- for k=0:5
- n=2^k;
- h = 1/n;
- h_1(k+1)=h;
- for i=1:n+1
- w(i)=fun(i*h);
- end
- %
- h = 1/n; % h stellt die Schrittweite dar, u den Anfangswert; beide Werte sind durch die Aufgabe vorgegeben
- u = 0;
- t=[0:h:n*h]; %Vektor erstellen mit Schrittweite h, der bei t_0 beginnt und bei t_(f-1) endet
- v(1)=u; %X_0 wird festgelegt, dieser Wert entspricht dem Anfangswert
- for i=2:(n+1) %Beginn bei X_1 da X_0 bereits festgelegt, danach X_2 etc. bis die festgelegte Schrittanzahl erreicht ist (ist immer verschoben, weil Vektoren keinen 0ten Eintrag haben)
- v(i)= v(i-1)+h*DGL(t(i-1),v(i-1)); %X_i wird festgelegt nach dem Eulerverfahren; DGL wird in einem extra File festgelegt und dann hier abgerufen, man kann jedoch alle möglichen DGL in dem File einsetzen
- end
- %
- for i=1:n+1
- e_1(i)=(v(i)-w(i))^2;
- end
- e(k+1)=sqrt((sum(e_1)));
- end
- hold off
- title('Fehler');
- loglog(h_1,e,'y')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement