function [] = vexconv(vex, r) r = (1/r + 1/250)^-1;