1 | initial version |

One option is to use `find_fit`

.

Here is some sample code.

```
# example data
f(x,y) = x^2+y^2
data=[]
for i in range(0,10):
for j in range(0,10):
data.append([i,j,f(i,j)])
# do the fitting
var('a b c')
model(x,y)=a*x^2+b*x*y+c*y^2
find_fit(data,model)
```

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.