prob5_5.m

来自「有关的矩阵处理的matlab源程序 较多」· M 代码 · 共 8 行

r=0.5;
dx=r/10;
x=-r:dx:r;
y=sqrt(0.25-x.^2);
z=x+0.5;
f=2*y.*z;
sum(f-0.5*(f(1)+f(end)))*dx