第一次练习
教学要求：熟练掌握Matlab软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab软件解决微积分、线性代数与解析几何中的计算问题。
补充命令
vpa(x,n)            显示x的n位有效数字，教材102页
fplot(‘f(x)’,[a,b])    函数作图命令，画出f(x)在区间[a,b]上的图形
在下面的题目中为你的学号的后3位（1-9班）或4位（10班以上）
1.1 计算与
程序：
syms x
limit((1001*x-sin(1001*x))/x^3,x,0)
结果：
程序：
syms x
limit((1001*x-sin(1001*x))/x^3,x,inf)
结果：
0
1.2 ，求
程序：
syms x
diff(exp(x)*cos(1001*x/1000),2)
结果：
南京有什么大学-2001/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/1000*x)
1.3 计算
程序：
dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)
结果：
1.4 计算
程序：
syms x
int(x^4/(1000^2+4*x^2))
结果：
1.5
程序：
syms x
diff(exp(x)*cos(1000*x),10)
结果：
1.6 给出在的泰勒展式(最高次幂为4).
程序：
syms x
taylor(sqrt(1001/1000+x),5)
结果：
1.7 Fibonacci数列的定义是，
用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。
程序：
x=[1,1];
for n=3:20
    x(n)=x(n-1)+x(n-2);
end
x
结果：
Columns 1 through 10
          1          1          2          3          5          8          13          21          34          55
  Columns 11 through 20
          89        144        233        377        610        987        1597        2584        4181        6765
1.8 对矩阵，求该矩阵的逆矩阵，特征值，特征向量，行列式，计算，并求矩阵（是对角矩阵），使得。
程序与结果：
a=[-2,1,1;0,2,0;-4,1,1001/1000];
inv(a)
eig(a)
  2.00000000000000
[p,d]=eig(a)
p =
  0.3355 - 0.2957i  0.3355 + 0.2957i  0.2425         
        0                  0            0.9701         
  0.8944            0.8944            0.0000         
注：p的列向量为特征向量
d =
  -0.4995 + 1.3223i        0                  0         
        0            -0.4995 - 1.3223i        0         
        0                  0            2.0000       
a^6
  11.9680  13.0080  -4.9910
        0  64.0000        0
  19.9640  -4.9910  -3.0100 
1.9 作出如下函数的图形（注：先用M文件定义函数，再用fplot进行函数作图）：
函数文件f.m:
function y=f(x)
if 0<=x&x<=1/2