I have a question about Pymc3 fit the equation:

S = A* f**(-a)+C

and I have known : log10(A) obeys Uniform(-10,10)

log10(c) obeys Uniform (-20,10)

a abeys Uniform (1,4)

And likelihood just a normal distribution,

How can I use Pymc3 determine the A ,a ,C (and I have the data about S, and f , like S=[1,2,3…], f=[1,2,3…])

Thanks!