也许像这样
信谊S C C t y
[T、Y] = meshgrid (10:0.5:10 10:0.5:10);
S = ((2。* y - t) / (2。* T-Y))。
dT = 1(大小(S));%方程
大叶= s /√(dT。^ 2 + s ^ 2);
dxu = dT. /√(dT。^ 2 + s ^ 2);
箭袋(T Y dxu、大叶)
包含“t”,ylabel“y”
标题(的方向场dy / dt = (2 * y - t) / (2 * t - y) ')
集(gca),“XLim”(10 0),“YLim”10 [0]);
这不是一个简单的颂歌,奇点在某些领域。可以用dsolve无法联系到,可以达到解决数值歌唱