01-10-2015, 03:52 AM

fsolve([x+y=3,2*x+2*y=6,x-y=0],[x,y],[0,0]); --> invalid dimension

fsolve([x+y=3,2*x+2*y=6,x-y=0],[x,y]); --> [3/2, 3/2]

fsolve([x+y=3,x-y=0],[x,y],[0,0]); --> [1.5, 1.5]

fsolve() seems use a naive dimension check -- if the number of equations is not equal to the number of variables, then giving initial values results in invalid dimension.

fsolve([x+y=3,2*x+2*y=6,x-y=0],[x,y]); --> [3/2, 3/2]

fsolve([x+y=3,x-y=0],[x,y],[0,0]); --> [1.5, 1.5]

fsolve() seems use a naive dimension check -- if the number of equations is not equal to the number of variables, then giving initial values results in invalid dimension.