I have a simple system which has 3dof. Stiffness matrix of this system depends on dof's displacement such as ki=k0*[1-0.1*sqrt(ui)]. I would like to solve this problem using ode45. But I could not manage this for MDOF systems. Any help will be greatly appreciated.