%solution for 3rd order ODE, reduced to canonical form
trange=[0 10];
x0 = [0 0.5 0];
[t x] = ode23(@myproblemmm,trange,x0)
plot(t,x(:,1)) %plots only first column since it is the solution for the 3rd order ODE