Matlab - Porblemas con ode

Hola, tengo el siguiente sistema de ecuaciones diferenciales:

function [f]=3bodies(t,x)

%Parametros
u=0.001;

% Ecuaciones
f(1)=x(2);

f(2)=2*x(4)+ x(1) -(1-u)*(x(1)+u)/((x(1)+ u)^2 + x(3)^2)^1.5 - u*(x(1)+ u-1)/((x(1)+ u -1)^2 + x(3)^2)^1.5;

f(3)=x(4);

f(4)=x(3)-2*x(2)-(1-u)*x(3)/((x(1)+ u)^2 + x(3)^2)^1.5 - y(3)*u/((x(1)+ u-1)^2 + x(3)^2)^1.5;

f = f';

Y al ejecutar >tr=[0,10];x0=[1,1,0,0];[t,x]=ode23('3bodies',tr,x0)

me aparece:

??? Invalid function name '3bodies'.

Error in ==> funfunprivateodearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...


Ojala me puedan ayudar!


Gracias

Es un problema con facil solución, el feval salta porque no puede evaluar los parametros que le metes, como son vecores pon un punto antes de ^, / y * así operará elemento a elemento

vale gracias ya lo arregle.