The element to be simulated is divided into regions, and each region is further divided into a plurality of quadrature nodes. Pairs are formed for all the quadrature nodes. Green's functions are computed and stored for the pairs. Each of the pairs is allocated to either the far field or the near field for purposes of simulation in accordance with a criterion. A Gaussian quadrature is computed for the pairs allocated to the far field while a high order quadrature is computed for those allocated in the near field. The component simulation is arrived after combining information derived from the Gaussian quadrature and the high order quadrature into a matrix which is then solved to obtain the charge distribution. Summation of the charges thus obtained yields the capacitance of the element. The high order quadrature is computed using a plurality of basis functions. The basis functions, denoted .psi..sub.i.sup.k (r'), are 1,x,y,x.sup.2,xy,y.sup.2. The basis functions are used to compute a set of weights v.sub.j.sup.k. The weights are computed by solving ##EQU1## where .psi..sub.i.sup.k (r') are the basis functions, G(r,r') are the Green's functions for each of the pairs allocated to the near field separated by a distance r-r', and a.sub.ij is a matrix satisfying the relationship a.sub.ij =G(r,r.sub.j.sup.k).psi..sub.i (r.sub.j.sup.k) and where index k counts the regions in the near field T.sub.k, index i counts the number of pairs, and index j counts up to a number p of the quadrature nodes in the near field. The advantage of this approach is that a fast solution to an integral equation descriptive of the element to be simulated can be achieved in exchange for constructing quadratures.

Das simuliert zu werden Element wird in Regionen geteilt, und jede Region wird weiter in eine Mehrzahl der Quadraturnullpunkte geteilt. Paare werden für alle Quadraturnullpunkte gebildet. Grünfunktionen werden für die Paare berechnet und gespeichert. Jedes der Paare wird entweder weit auffangen zugeteilt, oder die nahen fangen zwecks der Simulation in Übereinstimmung mit einem Kriterium auf. Eine Gaußsche Quadratur wird für die Paare berechnet, die weit auffangen zugeteilt werden, während eine hoher Auftrag Quadratur für die berechnet wird, die im nahen zugeteilt werden, auffangen. Die Teilsimulation wird angekommen, nachdem man die Informationen kombiniert hat, die von der Gaußschen Quadratur und von der hoher Auftrag Quadratur in eine Matrix abgeleitet werden, die dann gelöst wird, um die Aufladung Verteilung zu erreichen. Summierung der Aufladungen folglich erhaltenen Ergebnisse die Kapazitanz des Elements. Die hoher Auftrag Quadratur wird mit einer Mehrzahl der Grundlage Funktionen berechnet. Die Grundlage Funktionen, bezeichnetes psi..sub.i.sup.k (r '), sind 1, x, y, x.sup.2, xy, y.sup.2. Die Grundlage Funktionen werden verwendet, um einen Satz Gewichte v.sub.j.sup.k zu berechnen. Die Gewichte werden berechnet, indem man ## EQU1 ## löst, in dem psi..sub.i.sup.k (r ') die Grundlage Funktionen sind, G(r, sind r ') die Grünfunktionen für jedes der Paare, die dem nahen zugeteilt werden, auffangen getrennt durch eine Abstand Eisenbahn ', und a.sub.ij ist eine Matrix, die das Verhältnis a.sub.ij = G(r, r.sub.j.sup.k).psi..sub.i (r.sub.j.sup.k) erfüllt und wo Index k zählt, fangen die Regionen im nahen T.sub.k, Zählimpulse des Index I die Zahl Paaren auf, und Index J zählt bis zu einer Anzahl p von den Quadraturnullpunkten im nahen auffangen. Der Vorteil dieser Annäherung ist, daß eine schnelle Lösung zu einer Integralgleichung, die vom simuliert zu werden Element beschreibend ist, gegen das Konstruieren von von Quadraturen erzielt werden kann.