In this talk, we present three new classes recently introduced in YUIMA package. These classes allow the user to manage three different problems: • Construction of a multidimensional stochastic differential equation driven by a general multivariate Lévy process. In particular we show how to define and then simulate a SDE driven by a multivariate Variance Gamma process. • Definition and simulation of a functional of a general SDE. • Definition and simulation of the integral of an object from the class yuima.model. In particular, we are able to evaluate Riemann Stieltjes integrals, deterministic integrals with random integrand and stochastic integrals. Numerical examples are given in order to explain the new methods and classes.