matlab符号运算的一个问题

starbinbin_csu

如下代码，eq表示一个过定点(-13.5,6),(-7.5,0)的圆系,c是未知数，然后我求出了该圆系关于水平对称轴下半轴 y=f(x)的关系式，但是却发现f(-13,5)不等于6，f(-7.5)不等于0,请问是怎么一回事呢？

1. syms x y c
   
2. eq=(y-0).*(y-6)+(x+7.5).*(x+13.5)+c.*((x+7.5).*(6-0)-(y-0).*(-13.5+7.5));
   
3. result=solve(eq,'y');
   
4. fun1=char(simplify(result(1)));%函数表达式
   
5. fun=eval(['@(c) @(x)' vectorize(fun1)]);

复制代码