# MATLAB: Can someone help with this

differential equations

Hi I need help plotting the resopnse to this differential equation with 10^4 points. thanks if you can help.
!x!+ cx! + kx + εx3 = Bcost0 ≤ c ≤ 0.5, −1 ≤ k ≤1, 0 ≤ ε ≤1, B = 8.5 ignore the exclamation points.

(1 + c + k) * x + ε*x^3 = B*cost
roots([epsilon, 0, (1 + c + k), -B*cost])
N = 100;[EPSILON, CK1] = ndgrid(linspace(0,1,N), linspace(1 + 0 + -1, 1 + 1/2 + 1, N));X = roots([EPSILON, 0, CK1, -B*cost]);surf(EPSILON, CK1, X(:,:,1));  same for (:,:,2) and (:,:,3) for the three different X values