Advertisement
Guest User

Untitled

a guest
Jan 16th, 2018
84
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.90 KB | None | 0 0
  1. clear
  2. for k=0:5
  3. n=2^k;
  4. h = 1/n;
  5. h_1(k+1)=h;
  6. for i=1:n+1
  7. w(i)=fun(i*h);
  8. end
  9. %
  10. h = 1/n; % h stellt die Schrittweite dar, u den Anfangswert; beide Werte sind durch die Aufgabe vorgegeben
  11. u = 0;
  12. t=[0:h:n*h]; %Vektor erstellen mit Schrittweite h, der bei t_0 beginnt und bei t_(f-1) endet
  13. v(1)=u; %X_0 wird festgelegt, dieser Wert entspricht dem Anfangswert
  14. 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)
  15. 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
  16. end
  17. %
  18. for i=1:n+1
  19. e_1(i)=(v(i)-w(i))^2;
  20. end
  21. e(k+1)=sqrt((sum(e_1)));
  22. end
  23. hold off
  24. title('Fehler');
  25. loglog(h_1,e,'y')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement