Desdeaquí, Estoy tratando de resolver un sistema simbólico de ecuaciones como esta

syms x y z;
[x, y, z] = solve('z = 4*x', 'x = y', 'z = x^2 + y^2')
x =
0
2

y =
0
2

z =
0
8

excepto que mis ecuaciones se generan en diferentes puntos del archivo m y con coeficientes aleatorios. Mi pregunta es cómo puedo lograr lo siguiente ...

// Generate the first equation.
n = *random number generated here*;
E1 = (z == n*x + 2*n);   // <--- How to save this symbolic equation to use in "solve(...)" later?

// Other work, generate other eqs.
...

// Solve system of eqs.
[x, y, z] = solve( E1 , E2, E3)   // What/How to call/use the previously saved symbolic equations.

¡Gracias por la ayuda!

EDITAR

Actualizado para explicar mejor el objetivo.