syms x y z;

z = x^2;

dz = diff(z,x);

iz = int(z,x);

iz2 = int(z,x,0,2);

eq = 'x^2 - 3*x + 2 = 4';

xx = solve(eq);

xxnum = vpa(xx);

eq1 = 'y - x^2 = 1';
eq2 = 'y - 3*x = 4';
[x,y] = solve(eq1,eq2);
vpa(x)
vpa(y)
pretty(y)