sztuka142.txt

%%zadanie 1
% clc;clear all;close all;
% f=inline('(x-(x-exp(-x))')
% xp=0.5; eps0=10^-3; max0=100; a=0; b=1;
% [fp,it]=fixpoint(f,a,b,xp,eps0,max0,1);

% %%zadanie 2
% clc;clear all;close all;
% f=inline('x-(-4+3*x-0.5*x^2)')
% xp=4; eps0=10^-3; max0=100; a=1; b=5;
% [fp,it]=fixpoint(f,a,b,xp,eps0,max0,1);

% %%zadanie dodatkowe
clc;clear all;close all;
f=inline('r*x*(1-x)')
% xp=0.21; eps0=10^-6; max0=100; a=0; b=1;
% r=3.2;
% [fp,it]=fixpoint(f,a,b,xp,eps0,max0,1,r);
% 
% figure
% plot(fp)

%%wykresik
rz=1:0.01:4; nr=length(rz)  ; max0=1000; eps0=10^-6; p0=0.2; a=0; b=1;
%% tyle jest r-n
figure
hold on
xmin=1; xmax=4; ymin=0; ymax=1;%%granica dla osi
axis([xmin xmax ymin ymax])
for j=1:nr
    r=rz(j);
    [fp,it]=fixpoint(f,a,b,p0,eps0,max0,0,r); %nie chcemy rysowac funkcji;
    if it<max0  %%je?li max ilo?c iteracji=1000
        fw=fp(it); %%fw-to, co kre?limy
    else
        fw=fp(800:max0);
    end
    n=length(fw);
    rr=r*ones(1,n); %%ci?g punkcik?w w jednej linii
    plot(rr,fw,'b.','markersize',1)
    drawnow
   
end