syms y(t) ode = diff(y,t) == t; ini_con = y(0)==10; ysol(t)= dsolve(ode,ini_con); ezplot(ysol,[0,10]);