Jeremy Knight Jacobian Approximation
<Text-field style="Heading 1" layout="Heading 1">Jacobian(n,f,x)</Text-field> Jacobian:=proc(n,f::Array,x::Array) #Approximates the value of the Jacobian at x #F is an array of functions, x is an array of values local i,j,e,h,J: J:=Array(1..n,1..n): e:=Array(1..n): h:=.0001: for i from 1 to n do for j from 1 to n do e[j]:=h: J[i,j]:=(f[i](x+e)-f[i](x-e))/(2*h): e[j]:=0: end do: end do: return(J): end proc: