请教数学物理系的学生二道难题

Tuexy

我碰到了两个问题，已经多天查找，没有结果。请大家帮我解决这两个难题。
第一个问题是这样：
根据图一的方程（其中T0 = 0.1, epsilon = 0.1 a = 0.1 0 <= b <= 4.0; 0 <= c <= 100），求出图三中的相图，它的坐标为b, c, 对应的稳定点区域和非稳定的相图区域边界。I：单稳极限环，II单稳，III：双稳
我的想法是先求jacobian，其行列式为０，　然后求含有b，　c　参数的方程　（左边微分用０代替），最后将三个方程联合，得到bc的一个函数，然后ezplot就可以得到最后的图。
但是我不会用solve解三个方程，四个参数的非线性方程，方程中不包含x(1), x(2). 下面是我的代码，但是在solve的地方我不会了。

1. %Condition:
   
2. %(1) F1 = F2 = 0 (2) Jacobian = 0
   
3. %
   
4. 
   
5. a = 0.1; epslion = 0.1; T0 = 5.0;
   
6. %use the symbolic calcution
   
7. syms x1 x2 y b c
   
8. 
   
9. F = [((epslion ^ 2 + x1 ^ 2 ) / ( 1 + x1 ^ 2) ) / (1 + x2) - a * x1;
   
10.     b / T0 - x2 / ( T0 * (1 + c * x1 ^ 2) )];
    
11. V = [x1, x2];
    
12. 
    
13. %calculate the jacobian matrix
    
14. J = jacobian(F, V);
    
15. J_det = det(J);
    
16. 
    
17. %function fc = g(x)
    
18. y(1) = ((epslion ^ 2 + x1 ^ 2 ) / ( 1 + x1 ^ 2) ) / (1 + x2) - a * x1;
    
19. y(2) = b / T0 - x2 / ( T0 * (1 + c * x1 ^ 2) );
    
20. y = [J_det y(1) y(2)];
    
21. 
    
22. result = solve(y, 'b', 'c');

复制代码

第二个问题，是图２的偏微分。其中Ｄ＝　０，　右边只有前两项。matlab中没有这种类型的方程，它只能解elliptic, parobalic 之类的，然后PDE tool双只能解两阶，我这个是一阶三参数方程（t, x(1), x(2)）这三个变量。
它的边值条件是反射边值条件。J(x, t) = FP − D· 偏P／偏x　＝　０　（t在任一时候都是）。这个好像第三类边界条件。