Set of Equations for BVP

Hello everybody, I would like to solve the fallowing:

dxds= cos(φ), dyds= sin(φ), dφds= κ, dκds= F * cos(φ), dFds= 0

with

x(0) = 0 y(0) = 0 κ(0) = 0 φ(0) = n*pi/16 φ(1) = -pi/2

but my own written code does not work. I used

y1 = x, y2 = y, y3 = φ, y4 = κ

Can anybody tell me where my mistake is?

function test

F = 1;
solinit = bvpinit(linspace(0,1,1000),[0,0,pi,0],F);
sol = bvp4c(@mat4ode,@mat4bc,solinit,options);

function dydx = mat4ode(x,y,F)
dydx = [  cos(y(3))
          sin(y(3))
          y(4)
          F*cos(y(3))];

function res = mat4bc(ya,yb,F)
res = [  ya(1) 
         ya(2) 
         ya(4)
         ya(3)-pi/2
         yb(3)+pi/2];

[Torsten]

You did not tell us what you mean by "my own written code does not work" .

Best wishes

Torsten.