U radu je prikazan algoritam za dvodimenzionalnu (2D) metodu rubnih elemenata (BEM) analizu potencijalnih Laplaceovih problema s heterogenom domenom. Tehnikom dvostrukog i višestrukog globalnog čvora omogućena je točna numerička aproksimacija normalne gustoće toka u točkama s njezinim fizičkim diskontinuitetom. Dobiveni sustav linearnih jednadžbi proširen je primjenom dodatnog skupa linearnih jednadžbi. Ove se jednadžbe dobivaju izjednačavanjem potencijala u istodobnim globalnim čvorovima. Dodatno je povećana točnost postupka novim izrazima dobivenim analitičkom integracijom duž linearnih rubnih elemenata. Jedan primjer s analitičkim rješenjem ilustrira razvijeni algoritam.