对于初学matlab的人来说,如何编写数学函数可能有点问题,因为matlab是以矩阵为支撑的,所以写matlab 函数时必须有矩阵的思想
比如下面的1,2,4阶牛顿—科特斯公式
clear
a=1;
b=3;
x=[a,(3*a+b)/4,(a+b)/2,(a+3*b)/4,b]
f=x.*x.*x-2.*x.*x+7.*x-5;
error1=(b-a)/2*(f(1)+f(5))-20-2/3
error2=(b-a)/6*(f(1)+4*f(3)+f(5))-20-2/3
error4=(b-a)/90*(7*f(1)+32*f(2)+12*f(3)+32*f(4)+7*f(5))-20-2/3
简单来说,对于一个数学函数,要想求其函数值,就必须先把自变量用矩阵的形式写出来,然后函数值也是一个矩阵的形式,想用哪个直接用下标访问就行。