Substitute all solutions of a give system of equations.

I would like to substitute all the solutions of a system of equations in the original system. For example I would like to do something like this in Mupad:

eqs:=[x^2-1,y^2-1]

solutions:=solve(eqs)

and then create a matrix which contains all the couples of solutions; I tried with evalAt and subs combined with map and op but had no luck.

Any suggestions?

[Radu Trimbitas]

The idea is to consider the solution as a list and to evaluate the equation for the elements given in that list. The Mupad code is

eqs:=[x^2-1,y^2-1]:
solve(eqs,{x,y}):
zz:=op(%):
f:=l->(eqs)|l:
map([zz],f)

Answer 1

[Azzi Abdelmalek]

syms x y
eqs=[x^2-1,y^2-1]
sol=solve(eqs)
M=[sol.x sol.y]

[Azzi Abdelmalek]

[M(:,1).^2-1,M(:,2).^2-1]

Answer 2

[Alessio Nava]

And what about replacing x and y in eqs with all the couples contained in the vector M?

[Alessio Nava]

Anyone?