I have a symbolic function, I want to know the points of discontinuity.

 For example the function is : 

syms t
 k=t*heaviside(5-t)+5*heaviside(t-5) ;

Answer 1

[Walter Roberson]

http://www.mathworks.com/help/symbolic/mupad_ref/discont.html

so,

syms t
k = t*heaviside(5-t)+5*heaviside(t-5);
disc = feval(symengine, 'discont', k, t);

[Kanav]

can we extend this function for a two variable function?