MATLAB: Second order differential equation, I get a wrong result and i don’t understand what I did wrong…

iiode

Report my code. I want to solve this II order diff. eq. and as results I get y as a 2×11 double with all the values equal to zero.
It shouldn't be like that. Am I using a wrong approach? Where am I missing something ?
My code:
syms y(t);
C=8.847;
alpha=0.337; 
s_=16.86;  
d=9.5; 
A=70.9;
E=41400; 
k=pi*d/(E*A);
ode = diff(y,2) == k*(1-y/s_)*C*y^alpha;
V = odeToVectorField(ode);
F = matlabFunction(V,'vars',{'t','Y'});
sol = ode45(F,[0 250],[0 0])
x = linspace(0,250,2500);
y = deval(sol,x,1);
plot(x,y)

Best Answer

If you prescribe y(0)=y'(0)=0 as initial conditions, the solution y(t)=0 for all t is correct :-)
Best wishes
Torsten.
Related Question