In this article, we examine a stochastic partial differential equation (SPDE) driven by a symmetric $α$-stable (S$α$S) Lévy noise, that is multiplied by a linear function $σ(u)=u$ of the solution. The solution is interpreted in the mild sense. For this models, in the case of the Gaussian noise, the solution has an explicit Wiener chaos expansion, and is studied using tools from Malliavin calculus. These tools cannot be used for an infinite-variance Lévy noise. In this article, we provide sufficient conditions for the existence of a solution, and we give an explicit series expansion of this solution. To achieve this, we use the multiple stable integrals, which were developed in Samorodnitsky and Taqqu (1990, 1991), and originate from the LePage series representation of the noise. To give a meaning to the stochastic integral which appears in the definition of solution, we embed the space-time Lévy noise into a Lévy basis, and use the stochastic integration theory (Bichteler and Jacod 1983, Bichteler 2002) with respect to this object, as in other studies of SPDEs with heavy-tailed noise: Chong (2017a), Chong (2017b), Chong, Dalang and Humeau (2019). As applications, we consider the heat and wave equations with linear multiplicative noise, also called the parabolic/hyperbolic Anderson models.