求几道matlab编程题

clc
clear all
close all
syms x y z k p3 p2 p1
%% 第12题
jacobian((sin(x*y)+cos(y*z)+tan(x*z)),[x,y])
% ans =
%
% [ z*(tan(x*z)^2 + 1) + y*cos(x*y), x*cos(x*y) - z*sin(y*z)]
%% 第13题
integral2(@(x,y)(2*x.^2.*sin(y)+3*x.*y.^2),0,3,0,3)
% ans =
%
% 157.3199
%% 第14题
% fun=(k+1)*p3-((2*k+1)*x*p2-k*p1);
% solve(fun,p3)
% ans =
%
% (p2*x - k*p1 + 2*k*p2*x)/(k + 1)
%上面是求出p(k+1)的表达式
k=10;
p(1,k+1)=x;
p(0+1)=1;%p_0(x)=1;
p(1+1)=x;
for k=2:10
p(k+1)=(p(k)*x - k*p(k-1) + 2*k*p(k)*x)/(k + 1);
end
for k=1:11
p(k)
end
%上面这3句用于显示答案，答案太长，我这里只列出p0到p3
% ans =
%
% 1
%
%
% ans =
%
% x
%
%
% ans =
%
% (5*x^2)/3 - 2/3
%
%
% ans =
%
% (x*((5*x^2)/3 - 2/3))/4 - (3*x)/4 + (x*(10*x^2 - 4))/4