%% constitutive model Drucker-Prager屈服函数 %% = （A+B(1+C*EXP(nε））*（1+b1*η） %% constitutive model Modified-Drucker-Prager屈服函数 ... %% = （A+B(1+C*EXP(nε））*（a1+a2*exp(a3η）) %% constitutive model parameter A=1; B=2; C=3; n=1; %% Drucker-Prager Parameter b1 = 0.1; %% Modified-Drucker-Prager Parameter a1 = 1; a2 = 2; a3 = 3; %% plastic strain ss = 1; %% stress state tria = sqrt(3)/3; %% Drucker-Prager屈服函数 f1 = (A+B(1+C*EXP(n*ss)))*(1+b1*tria); %% Modified-Drucker-Prager屈服函数 f2 = (A+B(1+C*EXP(n*ss)))*(a1+a2*exp(a3*tria)); 将这段程序改成，输入矩阵ss，求出f1和f2的值

可以将原有的程序改成以下形式：

% constitutive model parameters
A = 1;
B = 2;
C = 3;
n = 1;
b1 = 0.1;
a1 = 1; 
a2 = 2;
a3 = 3;

% input matrix ss
ss = [1 2 3; 4 5 6; 7 8 9];

% stress state tria
tria = sqrt(3)/3;

% calculate f1 and f2
f1 = (A + B.*(1+C.*exp(n.*ss))).*(1+b1.*tria);
f2 = (A + B.*(1+C.*exp(n.*ss))).*(a1+a2.*exp(a3.*tria));

这样，输入矩阵 `ss` 就可以是任意大小的矩阵，`f1` 和 `f2` 的值也会相应地计算出来，并以大小与 `ss` 相同的矩阵形式返回。