syms x y z
f = y * sin(x) + z * cos(x);
int1 = int(f,z, z_1, z_2); % z_1, z_2 分别为z的上下限
int2 = int(int1, y, y_1, y_2); % y_1, y_2 分别为y的上下限
int_f = int(int2, x, x_1, x_2); % x_1, x_2 分别为x的上下限
你是要算解析解,还是数值解,数值解有n种方法,例如euler ,runge-kutta,