function y=F(p)

% This function computes the function to be integrated by the
% Runge-Kutta routine.
%
% Translated into MATLAB by Ronald Kumon, 21 Oct 1995

% Declare global variables
global a n nhar

% Function to integrate
y=-a.*p + j/4*n.*PconvP(p,nhar);