clear;clc
syms x y
eq=4*x*y+x^2-y;
y=solve(eq,y);
dy=subs(diff(y),x,[3,4,6,8])