
function R = myfun3 (z, y)
kh_z = kh_z2 (z);
R = [y(2) ; -kh_z*y(1)];
end




