# dsolve2 is intended to solve differential equations with
# initial conditions and return the resulting function.
#
# Example: 
#
#   > de := {diff(y(x),x) = 1,y(0) = 1}:
#
#   > f := dsolve2( de, y(x) ):
#
#   > f(2);
#     3;
#
# Note that you have to tell dsolve what the independent and dependent
# variables are by telling it to solve with respect to y(x).
#
# Definition of dsolve2:

dsolve2 := (de, fexpr) -> unapply( rhs( dsolve( de, fexpr )), op(1,fexpr)):