# Use assume() on find_fit

Hi everyone,

does anyone know how I can use assumptions with find_fit? A minimal example would be

```
var('a,b,x')
assume(a>0,b>0)
find_fit([(2,2),(3,2.9),(4,4.1)],(a+1)/b*x,variables=[x],initial_guess=(-2,-2))
```

and the output is

```
[a == -105.66654155674746, b == -104.30685450333799]
```

Any help or hint on where to find a solution is appriciated. Thanks in advance.

Considering that `find_fit` uses Scipy and assumptions are in Maxima, I think it's quite unlikely that this would do much. I guess we could try to run them through assumptions, though this might slow it down appreciably.