function [t,y] = impeuler_sajat(diffegy, t0, y0, h, N)
t=zeros(N+1,1);
y = zeros(N+1,1);
t(1) = t0;
y(1) = y0;
for i=1:N
t(i+1) = t(i) + h;
g = @(x)(x-y(i)-h*diffegy(t(i+1),x));
%g_der = @(x)(1+h);   %Kezzel kiszamitottuk elotte a g masodik valtozo szerinti derivaltjat
%Newton-modszert alkalmazunk a cikluson belul
g_der = @(x)(1-h*(1-t(i+1)^2));  % y'=y(1-t^2)+2  diffegyenlethez
x(1) = y(i)+h*diffegy(t(i),y(i));   % kezdeti iteracio az explicit Euler-megoldas
for k=1:4
    x(k+1) = x(k) - g(x(k))/g_der(x(k));
end
y(i+1)=x(k+1);
end

