这样的solve有什么问题？
>> sym ;>> =solve('A+B+C+D+E=0','A*(x0-2*h)+B*(x0-h)+C*x0+D*(x0+h)+E*(x0+2*h)=0','A*(x0-2*h)^2+B*(x0-h)^2+C*x0^2+D*(x0+h)^2+E*(x0-2*h)^2=2','A*(x0-2*h)^3+B*(x0-h)^3+C*x0^3+D*(x0+h)^3+E*(x0+2*h)^3=6*x0','A*(x0-2*h)^4+B*(x0-h)^4+C*x0^4+D*(x0-h)^4+E*(x0+2*h)^4=12*x0^2','A','B','C','D','E')
??? Error using ==> mupadengine.mupadengine>mupadengine.feval at 156
Error: illegal unknown
Error in ==> solve at 157
R = feval(symengine,'mlfsolve',eqns,vars);
把x0删掉··可以运行···
=solve('A+B+C+D+E=0','A*(-2*h)+B*(-h)+C*0+D*(h)+E*(2*h)=0','A*(-2*h)^2+B*(-h)^2+C*0^2+D*(h)^2+E*(2*h)^2=2','A*(-2*h)^3+B*(-h)^3+C*0^3+D*(h)^3+E*(2*h)^3=0','A*(-2*h)^4+B*(-h)^4+C*0^4+D*(h)^4+E*(2*h)^4=0^2','A','B','C','D','E')
Warning: Explicit solution could not be found.
> In solve at 170
ANS:
A =[ empty sym ]
B = []
C = []
D = []
E = []
但在只有A、B、C、x0的时候可是完全正常的··
>> =solve('A+B+C=0','A*(x0-h)+B*x0+C*(x0+h)=1','A*(x0-h)^2+B*x0^2+C*(x0+h)^2=2','A','B','C')
ANS:
A =-(h + 2*x0 - 2)/(2*h^2)
B =(2*x0 - 2)/h^2
C =(h - 2*x0 + 2)/(2*h^2)
>> =solve('A+B+C=0','A*(x0-h)+B*x0+C*(x0+h)=0','A*(x0-h)^2+B*x0^2+C*(x0+h)^2=2','A','B','C')
ANS:
A =1/h^2
B =-2/h^2
C =1/h^2
[ 本帖最后由 ChaChing 于 2009-4-7 07:57 编辑 ] 试了下, 没报错!?
>> sym ;
>> =solve('A+B+C+D+E=0','A*(x0-2*h)+B*(x0-h)+C*x0+D*(x0+h)+E*(x0+2*h)=0','A*(x0-2*h)^2+B*(x0-h)^2+C*x0^2+D*(x0+h)^2+E*(x0-2*h)^2=2','A*(x0-2*h)^3+B*(x0-h)^3+C*x0^3+D*(x0+h)^3+E*(x0+2*h)^3=6*x0','A*(x0-2*h)^4+B*(x0-h)^4+C*x0^4+D*(x0-h)^4+E*(x0+2*h)^4=12*x0^2','A','B','C','D','E')
A =
1/4*(-4*x0^2*h^3+4*x0*h^4+20*x0^3*h^2-h^5+16*x0^5)/(3*h^5+5*x0*h^4+2*x0^2*h^3+2*x0^3*h^2-4*x0^4*h+4*x0^5)/h^2
B =
-(-9*x0^2*h^3+14*x0^3*h^2-4*h^5+8*x0^5)/(3*h^5+5*x0*h^4+2*x0^2*h^3+2*x0^3*h^2-4*x0^4*h+4*x0^5)/h^2
C =
1/2*(-30*x0^2*h^3-4*x0*h^4+28*x0^3*h^2-15*h^5+8*x0^5)/(3*h^5+5*x0*h^4+2*x0^2*h^3+2*x0^3*h^2-4*x0^4*h+4*x0^5)/h^2
D =
-(6*x0^3-4*h^3-7*x0^2*h)/(3*h^5+5*x0*h^4+2*x0^2*h^3+2*x0^3*h^2-4*x0^4*h+4*x0^5)
E =
1/4*(4*x0^3+4*x0*h^2-h^3)/(3*h^5+5*x0*h^4+2*x0^2*h^3+2*x0^3*h^2-4*x0^4*h+4*x0^5) 这就奇怪了。·。·。
页:
[1]