MATLAB: Where is the mistake with the code

ode45

it keeps the initial values when it runs instead of solving, where is the problem?
B=9020;
ko=5*(10^10);
[T,x]=ode45(@diaforikes_eksiswseis(t,x,ko,B),[0  60],[10 20 0  -5651950])
function ddt=diaforikes_eksiswseis (t,x,xs,ko,B)
ddt(1)=(-ko*exp(-(B/(x(4)-xs(4)))))*(x(1)-xs(1))
ddt(2)=0
ddt(3)=ko*exp(-(B/(x(4)-xs(4))))*(x(1)-xs(1))
ddt(4)=(375*ko*exp(-(B/(x(4)-xs(4))))*(x(1)-xs(1))*(x(4)-xs(4)-25)+(331380*ko*exp(-(B/(x(4)-xs(4))))*(x(1)-xs(1)))+(90000*(60-x(4)+xs(4))))/((125*(x(1)-xs(1)))+(250*(x(3)-xs(3)))+(75*(x(2)-xs(2))))
ddt=ddt';
end

Best Answer

Try this:
[T,x]=ode45(@(t,x) diaforikes_eksiswseis(t,x,ko,B),[0  60],[10 20 0  -5651950]);
Related Question