It would be a good idea to post the code instead of a textual description. But most likely "f = ncs1_dae(x)" is the problem already, because it accepts 1 input only, while the integrator provides at least the current position and the time, so you need 2 inputs.
The standard method to reveal such problems is using the debugger:
dbstop if error
Then Matlab stops when the error occurres and you can inspect the current variables and the calling hierarchy.
