Thanks, you are right the dimensions of c_2 are wrong.
I think in my function n the term m*g/c_1 (=constant) got absorbed by c_2 and in function f I forgot to do this. I’ll try to fix this and report back.
Edit:
I corrected my f-function, but unfortunately the solution is still a vector of NaNs