At the bottom of the vase the radius of the hole is bez[0].x but we need to find the radius at y = 2. The function bezier_line_intersection() returns a list of u-values where a given line intersects ...

These files are MATLAB and Python wrappers for the Okada DC3D0 point source and the DC3D rectangular dislocation surface fortran subroutines. The Matlab wrappers are written using the MEX functions.