MATLAB earthquake forward program