solve function of ode45

function dzdt= speed(x,z)
M=10; m=0.1;  c=0.2;  k=0.4; F=10*x ;
dzdt(1)=z(2);
dzdt(2)=-(c*z(1) - F + k*z)/(M + m);
dzdt=dzdt';
end

i want to use this function to plot y versus x versus

Answer 1

[Star Strider]

There is an error in the code you posted. I believe it should be:

dzdt(2,:)=-(c*z(1) - F + k*z(2))/(M + m);

Otherwise, this is a straightforward application of ode45 (link). See the function documentation for details.

Answer 2

[madhan ravi]

Try this:

There is nothing wrong with your code:

dzdt(2)=-(c*z(1) - F + k*z(2))/(M + m);
%                         ^^^---missed it