Advertisement
Guest User

Untitled

a guest
Sep 25th, 2017
60
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.69 KB | None | 0 0
  1. public static void trapezia(){
  2.  
  3.  
  4. double yi1 = 0;
  5. double yi = y0;
  6. double xi = x0;
  7.  
  8. point2.add(new Point(xi,yi));
  9.  
  10. for(int i = 1;i < N + 1;i++){
  11.  
  12. double yiCurr = yi;
  13. double yiPrev = yi;
  14.  
  15. do{
  16. yiPrev = yiCurr;
  17. yiCurr = yiPrev - (yiPrev - yi - 0.5 * h * (f(xi,yi) + f(xi + h,yiPrev))) / (1 - 0.5 * h * f1(xi + h,yiPrev));
  18. }
  19. while(Math.abs(yiCurr - yiPrev) > 0.0000001);
  20.  
  21. yi1 = yi + 0.5 * h * (f(xi,yi) + f(xi + h,yiCurr));
  22.  
  23. yi = yi1;
  24. xi = x0 + i * h;
  25.  
  26. point2.add(new Point(xi,yi));
  27. }
  28.  
  29. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement