clc
clear all
dsolve('Dy=1+y^2','y(0)=1','x')