clear all

syms x; %%deklaracja zmiennej x
I=0.3;
C=0.001;
U=5;
eps=0.01;

fun=(I/C)-U*100*pi* sin(100*pi*x);
%fun= x^2 - 5;
f=inline(fun);

z=diff(f(x));
df=inline(z); %inline coby mi wynikow nie pieprzylo, chuj wie co to robi
m=diff(diff(f(x)));
ddf=inline(m);

xo=0.001;
x=xo;

x1=(xo-(f(xo)/df(xo)))

%if (f(x)*ddf(x)>0)
   %disp('Funkcja spełnia warunki')
   N=0
while(abs(f(x1))>eps)
N=N+1
x1=(x1-(f(x1)/df(x1)))
end

disp('WYNIKI I SPRAWDZENIE')
disp('PRZEDZIAŁ')
xo
x1

disp('SPRAWDZENIE')
f(xo)*f(x1)
disp('JAK JEST <0 to OKEJ')



%else disp('Funkcja nie spełnia warunków')
%end