How to resolve a pde equation using pdetool extracted data

Hi

My question is how to resolve programmaticaly this equation (*)using pdetool data "p,e,t,b,,, etc" ,"Direclet and Neumann boundary's included" then using Jiles-Atherton method to plot hysteresis loop.

(*) [ (1/delta(t)) * [K] + [M] ] * {A}^(t+delta(t)) = (1/delta(t)) * [K] * {A}^t +{G}^t + {F}^t

with

Js = J(max)*sin(omega*t).

{F} = (delta(e)/3)*Js

delta(e) =1/2[(x1-x2)/(y1-y2)]

{G}^t = mu0*delta(e)[b1*My-C1Mx;b2*My-C2Mx;b3*My-C3Mx]

A : potential vector.

thank you