function deriv = corioliskraft3d(status,t,param)

% Dette er et miniprogram som brukes ved integrering av Newtons 
% 2. lov i mekanikken når akselerasjonen skyldes Coriolis og 
% sentrifugalkraft.
% Hensikten med denne rutinen er rett og slett å få gitt inn 
% likningsettet som benyttes som input i Runge-Kutta-rutinen rk4
% eller i rk4a (adaptiv RK).
% Input er statusvektor (posisjon og hastighet, begge i tre (!!!)
% dimensjoner), og returverdiene er den deriverte av posisjon 
% og hastighet, det vil si hastighet og akselerasjon.
% Programmet er skrevet av Arnt Inge Vistnes 1. mai 2005.

% Beregner akselerasjonen:
r = [status(1) status(2) status(3)];  % Henter ut posisjonen
v = [status(4) status(5) status(6)];  % og hastighet
omega(3) = param;                     % Henter omega fra parameterne
omega(1) = 0; omega(2) = 0;
a = -cross(omega,cross(omega,r))-2*cross(omega,v);
%a = -2*cross(omega,v);
                                      % sentrifugal- og
                                      % Coriolisakselerasjon
                                                       
% Returnerer de deriverte av statusvektor (dvs hast. og aksel.)                       
deriv = [v(1) v(2) v(3) a(1) a(2) a(3)]; 
return;