Dirk Laurie on Wed, 07 Aug 2013 16:26:09 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Variable checker |
I have written a routine that for a given expression checks whether it contains any variables that are not in a given list. It works. Only problem is, I don't understand why. { foreign(expr,vars) = if (type(expr)=="t_RFRAC", return(foreign(numerator(expr)) || foreign(denominator(expr)))); my(var=trap(,,variable(expr))); print(expr", "var); if(!var, return(0)); \\ expr is scalar if(subst(vars,var,0)==vars, return(var)); \\ var is not present in vars foreign(polcoeff(expr,0),vars) } addhelp(foreign,"foreign(expr,vars): 0 if expr can be written in terms of vars, otherwise an offending variable") ? foreign(x^2+a*x+b,[x,b]) x^2 + a*x + b, x b, a %2 = a The value of 'expr' on the recursive call is 'b'. There's no 'a' in it. Yet variable(expr) returns 'a'. I'm grateful it does. But it baffles me.