clear; clc;

syms x y;
y = cos(x)
A = [x, 1; 1, x]
z = A * [x; y]
w = subs(z, 2)
eval(w)

diff(z)
int(z)

% See https://www.maplesoft.com/products/Maple/
% See http://www.wolfram.com/?source=nav