Given a hypersurface $D = h^{-1}(0)$ for some polynomial $h \in \mathbb{C} [x,y,z]$ I want to be able to use Macaulay 2 to tell if it's a free divisor or not.
What I've got so far;
Let $h_{p}$ be the reduced equation for $h$. For $D$ to be a free divisor I need $Der(-logD) := \{ \delta : \delta(h_{p}) \in (h_{p})\}$ to be a locally free $\mathcal{O}_{\mathbb{C} ^{n}}$-module, where $\delta$ is a logarithmic vector field. So for a vector field $\delta = a \frac{\partial}{\partial{x}} + b \frac{\partial}{\partial{y}} + c \frac{\partial}{\partial{z}}$ I need $a \frac{\partial{h}}{\partial{x}} + b \frac{\partial{h}}{\partial{y}} + c \frac{\partial{h}}{\partial{z}} = r h_{p}$ for some $r$.
So I think If I remove the bottom row of ones that I get when when using Macaulay 2 to calculate the kernel of the map given by the matrix $\lbrack \frac{\partial{h}}{\partial{x}} , \frac{\partial{h}}{\partial{y}} , \frac{\partial{h}}{\partial{z}} , -h_{p} \rbrack$ I'll get the module $Der(-logD)$.
But I'm not sure how to then check if its a free module, I am slightly confused too as I've been told to use the resolution function in Macaulay 2 but when ever I try I can't seem to get it to work.
I've only just started using Macaulay 2 and I am fairly new to free divisors too so sorry for any errors and thanks in advance for any help.