solving Coupled nonlinear differential equation?

theta(s)=Const(c1*cos(theta_initialvalue-alpha)+(yfinal-y(s))*cos(alpha)-(xfinal-x(s)*sin(alpha))

diff(x(s))=cos(theta(s))

diff(y(s))=sin(theta(s))

tan(aplha)=(yfinal-c1*(1-cos(theta_initialvalue))/(xfinal-c1*sin(theta_initialvalue))

% initial conditions

theta(s=0)=0;

x(s=0)=0;

y(s=0)=0;

c1=0.002;

Const=F/EI

EI=3.06*10^3;

F=13.3;,

xfinal.yfinal,theta_initialvalue and alpha=0;

initially x(s),y(s),theta(s)

how to get analytical and numerical solution of these coupled non linear differential equations?