function [ yend ] = endpoint( C )
    y0 = [0, C];
    xmin = 0;
    xmax = 1;
    N = 100;
    [xout, yout]=rk4_integrator(y0, xmin, xmax, N, @f_ivp);
    yend=yout(end,1);
end

