Different a coefficient for each element not possible anymore with solvepde?

In Matlab 2015a I was using the pde solver for a problem in which the a coefficient was different for each Element:

u = assempde(problem,p,e,t,c,a,f);

'a' was a 1xN matrix with N the number of elements.

In version 2016a I get the remark that I should not use assempde anymore, but solvepde instead. So I was transferring all my code to the new functions. But when defining the coefficients:

specifyCoefficients(problem,'m',0,'d',0,'c',c,'a',a,'f',f);

I get the error: "Size of the matrix 'a' is not consistent with the pde.PDEModel/PDESystemSize property." Is it not possible anymore to define an a coefficient which is dependent on the element?

Answer 1

[Alan Weiss]

Sorry, but to use the latest solvers, I think that you will need to update your coefficient syntax to functional form. This form has the benefit that it works with quadratic elements as well as linear elements.

Alan Weiss

MATLAB mathematical toolbox documentation