two questions: A) list length of solve result, B) 0^0
Subject: two questions: A) list length of solve result, B) 0^0
From: Robert Dodier
Date: Sat, 7 Jun 2008 19:00:20 -0600
On 6/7/08, Stavros Macrakis <macrakis at alum.mit.edu> wrote:
> I suggest you call algsys directly rather than via solve for uniform
> result format.
Or linsolve, if the equations in question are linear.
The inconsistent return value of solve seems like a bug to me.
> If youb want to force 0^0 to simplify to 1, I believe you can do this
> using tellsimp, though at some cost in efficiency.
That seems to work as expected.
simp : false;
tellsimp (0^0, 1);
simp : true;
0^0;
=> 1
my_bernstein(n,k,x):=binomial(n,k)*(x+1)^k*(1-x)^(n-k);
my_bernstein(0,0,1);
=> 1
best
Robert