function [z] = fkt (x,y,c,n)


h = 0.01;

m = (hang (y(1)+h) - hang(y(1)-h)/(2*h));

alpha = atan (m);

z = [y(2);
c(1) * cos((pi/2)-alpha)*cos(alpha)];