A limit representation of a homogeneous symmetrical stable process with independent increments is obtained in the form of a sum of independent Wiener and purely step generalized Poisson processes. Simple and fast algorithms are developed for simulating symmetrical stable random variables and processes. The algorithms are suitable for software implementation.