Вы можете попробовать заменить вторую строку на:
f=['x+y+' num2str(k)];
Я не смог проверить это, так как у меня нет панели инструментов PDE.
Надеюсь, это поможет.
A.
PS: Если вам нужно решить для большого количества значений k, есть способы решить этот pde за один выстрел для всех значений k в данном интервале. В результате вы получите функцию трех переменных: u (x, y, k) (кто сказал, что k является параметром, а не переменной). Хорошим моментом является то, что стоимость памяти и время вычислений будут по-прежнему масштабироваться как 2D Пуассон, вы никогда не решаете / не сохраняете вещи в трехмерной сетке (где k - третье измерение).
Если вам интересны эти новые методы, которые называются: надлежащая обобщенная декомпозиция (PGD), вы можете найти несколько недавних работ по этим вопросам.