Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clear;
- N = 10;
- disp(N,"N=");
- a = -2.0;
- disp(a,"a=");
- b=0.0;
- disp(b,"b=");
- h=0.001;
- function f=f(x);
- f=x*exp(x)-(x^2);
- endfunction;
- function fp=fp(x,h);
- fp=(f(x+h)-f(x))/h;
- endfunction;
- function fpp=fpp(x,h);
- fpp=(fp(x+h)-fp(x))/h;
- endfunction;
- function g=g(x);
- g=(f(x+f(x))-f(x))/f(x);
- endfunction;
- wi=0;
- if f(a)*f(b)<0 & fp(a)*fp(b)>0 & fpp(a)*fpp(b)>0 then
- if f(a)*fpp(a)>0 then
- wi=wi+1;
- xp(wi)=a;
- disp(xp(wi));
- else
- wi=wi+1;
- xp(wi)=b;
- disp(xp(wi));
- end;
- for i=(wi+1):N
- xp(i)=xp(i-1)-f(xp(i-1))/g(xp(i-1));
- disp(xp(i));
- end;
- else
- disp("Założenia nie są spełnione");
- end;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement