There is a way to do this, but it might be the long way around. If you have the curves (graphs from the manufacturer), there is software that you can use to capture data points (x,y) from the graph itself. In the olden days, I used software that was free, but you had to click on each point you wanted and have it stored to a data set. Nowadays, there is software that will get the data points for you (Quintessa, for example).
After you have all of the data points on your curve, you can use an equation fitter (easily found online) that will give you a function that fits your curve. The last time I did this, Excel had a function that did this. It made up very nice equations for me that worked very well, with sometimes up to 4th power variables.
Avalanche and Zener diodes are discontinuous in their outputs, and I don't know if other programs will see this, but you can always use a Heaviside function to break it into parts.
So, you have two equations associated with two Heaviside functions that most decent math programs can easily deal with. That is, for V < a, you get one curve and for V >= a, you get a second. This can all be expressed in a single equation, if it's written properly.
There is a harder, but more precise, way to do this, but it will still involve two Heaviside functions. That, of course, is to write the equations according to the device physics and it's characteristics and doping. I think that most device manufacturers keep this information confidential.