  This paper is a demonstration of how the WKB approximation can be used tohelp solve the linearised 3D MHD equations. Using Charpit's Method and aRunge-Kutta numerical scheme, we have demonstrated this technique for apotential 3D magnetic null point, ${\bf{B}}=(x,\epsilon y -(\epsilon +1)z)$.Under our cold plasma assumption, we have considered two types of wavepropagation: fast magnetoacoustic and Alfv\'en waves. We find that the fastmagnetoacoustic wave experiences refraction towards the magnetic null point,and that the effect of this refraction depends upon the Alfv\'en speed profile.The wave, and thus the wave energy, accumulates at the null point. We havefound that current build up is exponential and the exponent is dependent upon$\epsilon$. Thus, for the fast wave there is preferential heating at the nullpoint. For the Alfv\'en wave, we find that the wave propagates along thefieldlines. For an Alfv\'en wave generated along the fan-plane, the waveaccumulates along the spine. For an Alfv\'en wave generated across the spine,the value of $\epsilon$ determines where the wave accumulation will occur:fan-plane ($\epsilon=1$), along the $x-$axis ($0<\epsilon <1$) or along the$y-$axis ($\epsilon>1$). We have shown analytically that currents build upexponentially, leading to preferential heating in these areas. The workdescribed here highlights the importance of understanding the magnetic topologyof the coronal magnetic field for the location of wave heating.
Using Charpit's Method and a Runge-Kutta numerical scheme, we have demonstrated this technique for a potential 3D magnetic null point, ${\bf{B}}=(x,\epsilon y -(\epsilon +1)z)$.