x = var('x')    # define a variable x
y = function('y',x)   # define y to be a function of that variable
DE = diff(y,x) == 1 + 4*y^2
print 'Differential Equation:'
show(DE)
soln = desolve(DE, [y,x], ics=(0,0))
print 'Solution:'
show(soln)

print 'Explicit solution:'
solve(soln,y)