% hw02cbd.m  -  Derivative for fcn of M604 HW02 C(b)
%
% Input parameter xy must be a column vector with two components.
% Output is a row vector of length 2.

function z = hw02cbd(xy)

z = zeros(1,2);

x = xy(1);
y = xy(2);

z = [4*(x-2)^3+2*(x-2)*y^2,  2*y*(x-2)^2+2*y+2];
