Advertisement
SonicDesu

Euler RR v2

Jan 13th, 2019
122
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.43 KB | None | 0 0
  1. clc;
  2. clear all;
  3.  
  4. f = @(x,y) x+y;
  5. yk=1;
  6. x0=0;
  7. eps=0.01;
  8. h=0.1;
  9.  
  10. k=0;
  11. xk=x0;
  12. while(xk=1+h)
  13.     yk1 = yk + hf(xk,yk);
  14.     yk2 = yk + 2hf(xk,yk);
  15.     if(abs(yk1-yk2)eps)
  16.             fprintf('k=%d x=%f y=%f  h=%fn',k,xk,double(yk),h);
  17.     xk=xk+h;
  18.     yk=yk1;
  19.     k=k+1;
  20.     else
  21.         fprintf('ZMIANA h abs = %f n',double(abs(yk1-yk2)));
  22.         h=h2;
  23.        
  24.     end
  25.  
  26. end
  27. % k=359 x=1.000000 y=3.426912  h=0.001563
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement