各个子程序如下:
function z=f_1(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) x.^(a+c-1).*y.^(b+d),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) x.^(a+c-1).*y.^(b+d),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_2(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) a*c*x.^(a+c-1).*y.^(b+d),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) a*c*x.^(a+c-1).*y.^(b+d),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_3(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) b*d*x.^(a+c+1).*y.^(b+d-2),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) b*d*x.^(a+c+1).*y.^(b+d-2),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_4(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) a*x.^(a+c-1).*y.^(b+d),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) a*x.^(a+c-1).*y.^(b+d),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_5(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) c*x.^(a+c-1).*y.^(b+d),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) c*x.^(a+c-1).*y.^(b+d),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_6(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) b*x.^(a+c).*y.^(b+d-1),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) b*x.^(a+c).*y.^(b+d-1),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_7(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) d*x.^(a+c).*y.^(b+d-1),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) d*x.^(a+c).*y.^(b+d-1),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_8(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) c*b*x.^(a+c).*y.^(b+d-1),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) c*b*x.^(a+c).*y.^(b+d-1),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_9(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) a*d*x.^(a+c).*y.^(b+d-1),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) a*d*x.^(a+c).*y.^(b+d-1),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
function z=f_10(a,b,c,d,yi,yo,r,h)
z1=quad2d(@(x,y) x.^(a+c+1).*y.^(b+d),r,1+h/2,-1,yo);
z2=quad2d(@(x,y) x.^(a+c+1).*y.^(b+d),r,1-h/2,-1,yi);
z=z1-z2;
z=vpa(z);
|