So Im doing a curve againist data plots for angular freq,

Ive reached the stage where im trying to evaluate the function,

ive put

>fit :=eval(f(x),R[2]);

where R[2] is the the values for the constants, in the function

F/sqrt(m^2*(omega^2-omega_0^2)^2+b^2*omega^2

which is also defined as f(x)