Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clc;
- clear all;
- f = @(x,y) x+y;
- yk=1;
- x0=0;
- eps=0.01;
- h=0.1;
- k=0;
- xk=x0;
- while(xk=1+h)
- yk1 = yk + hf(xk,yk);
- yk2 = yk + 2hf(xk,yk);
- if(abs(yk1-yk2)eps)
- fprintf('k=%d x=%f y=%f h=%fn',k,xk,double(yk),h);
- xk=xk+h;
- yk=yk1;
- k=k+1;
- else
- fprintf('ZMIANA h abs = %f n',double(abs(yk1-yk2)));
- h=h2;
- end
- end
- % k=359 x=1.000000 y=3.426912 h=0.001563
Advertisement
Add Comment
Please, Sign In to add comment