> plot(x^2-16,x=-7..7);
> deq1:=diff(x(t),t)=x(t)^2-16;
> DEtools[dfieldplot](deq1,x(t),t=0..2,x(t)=-7..7);
> plot(x^3-16*x,x=-7..7);
> deq2:=diff(x(t),t)=x^3-16*x;
> DEtools[dfieldplot](deq2,x(t),t=0..1,x(t)=-6..6);
> DEtools[dfieldplot](deq2,x(t),t=0..1,x(t)=-0.1..0.1);
>
DEtools[dfieldplot](deq2,x(t),t=0..1,x(t)=-4.1..-3.9);
> plot(x+sin(x),x=-20..20);
> deq3:=diff(x(t),t)=x(t)+sin(x(t));
> DEtools[dfieldplot](deq3,x(t),t=0..10,x(t)=-7..7);
>