subs效率特别低怎么破

我有一个syms的(含变量有560个)矩阵560560，矩阵是通过matlab微积分分求出来，带有变量560个，这560个变量分布在矩阵的每个位置上。
c = subs(J,q,p0); —J是生成的syms矩阵560560，q是含有560个变量组成的向量，p0是我要替换的数值组成的向量。subs处理这类问题效率极低，看过很多方法：
1、matlabFunction，效率更慢；
2、str2func，变量太多，@变量的时候需要一个一个敲，太麻烦；(若有人能解决也可以)
3、symfun,与2同理，变量太多，变量可以向量代替，但是输入不能用向量代替，需要一个一个敲。