function [J] = jacobian(x,t)
%
%  Jacobian of [H3.14]
%
A = [ 1250 -25113 -60050 -42647 -23999
       500 -10068 -24057 -17092  -9613
       250  -5060 -12079  -8586  -4826
      -750  15101  36086  25637  14420
       250  -4963 -11896  -8438  -4756 ];
%
J = A;
return