PDE toolbox: is it possible to use coordinate dependent coefficient c in parabolic function?

I wonder is it possible to use coordinate dependent coefficient c in parabolic function?

I am trying to solve single parabolic equation in two dimensions, where c is given as

[x^2 0; 0 y^2]

where x and y are coordinates of points in traingles:

x=p(1,:); y=p(2,:);

Looks like using of

c_pde = [x.^2 y.^2]; c_pdeCenter = pdeintrp(p,t,c_pde');

does not work correctly.

Help declares variety of methods to define c, a and f coefficients, but there is a lack of examples.