I'm currently working on models presenting high or infinite dimensionality, missing observations and non linearity. These issues make direct use of numeric approaches infeasible. I'm trying to find regularities and approximations in a symbolic form. I need support for basic operations like marginalization or Bayes rule and exploiting the properties of probability distribution (both general (>0 and sum up to 1) and specific (gaussian, poisson..). I need it to check that I am not missing and doing any transcription error. I know matlab has a symbolic extension, and I used a lot of years ago mathematica. But I did not find online any example using the specific features for probabilistic modelling that I need. Any suggestion is welcome