function [yp] = solutions_Final_p1a(t, y)

    if t>0
        yp = y*(-2*t+1/t);
    else
        yp = 1;
    end

end

