3.

> MODEL := diff( y(t), t ) = t / ( y(t) - t^2 * y(t) );

> IC := y(0) = 4;

(a)

(b)

(c)

(d)