Я решил Вашу систему в mathematica. Без сюрпризов. Если Вы не опечатались в условии.
-- 21 окт 2009, 00:02 --На фазовой плоскости несколько интереснее. Хотя тоже понятно.
Если у Вас установлена mathematica, то код простой
Код:
sol = NDSolve[{x'[t] == Log[1 + y[t] + x[t]] - x[t],
y'[t] == 1 - Sin[t x[t]] - y[t], x[0] == 0, y[0] == 0}, {x[t],
y[t]}, {t, 0, 10}]
Plot[x[t] /. sol, {t, 0, 10}]
Plot[y[t] /. sol, {t, 0, 10}]
Только с разделом Вы, ИМХО, немного ошиблись.