clear
clc
format long
syms x y;
%I(x,y)=p00+p10*x+p01*y+p20*x.^2+p11*x*y+p02*y.^2;(30)
I1(x,y)=-6.171e+04-1.768*x+-763.1*y+847*x.^2+0.2231*x*y+715.7*y.^2;
I11=int(int(I1(x,y),y,-3,0),x,-3,3);
I12=int(int(I1(x,y),y,0,3),x,-3,3);
detaI1=(I11-I12)./(I11+I12)
vpa(detaI1)
I2(x,y)=-6.163e+04-1.826*x-786.1*y+845.8*x.^2+0.2181*x*y+714.3*y.^2;
I21=int(int(I2(x,y),y,-3,0),x,-3,3);
I22=int(int(I2(x,y),y,0,3),x,-3,3);
detaI2=(I21-I22)./(I21+I22)
vpa(detaI2)
I3(x,y)=-6.154e+04-1.861*x-799*y+844.6*x.^2+0.2141*x*y+713.1*y.^2;
I31=int(int(I3(x,y),y,-3,0),x,-3,3);
I32=int(int(I3(x,y),y,0,3),x,-3,3);
detaI3=(I31-I32)./(I31+I32)
vpa(detaI3)
I4(x,y)=-6.146e+04-1.931*x-806.6*y+843.4*x.^2+0.2102*x*y+712*y.^2;
I41=int(int(I4(x,y),y,-3,0),x,-3,3);
I42=int(int(I4(x,y),y,0,3),x,-3,3);
detaI4=(I41-I42)./(I41+I42)
vpa(detaI4)
clc
format long
syms x y;
%I(x,y)=p00+p10*x+p01*y+p20*x.^2+p11*x*y+p02*y.^2;(30)
I1(x,y)=-6.171e+04-1.768*x+-763.1*y+847*x.^2+0.2231*x*y+715.7*y.^2;
I11=int(int(I1(x,y),y,-3,0),x,-3,3);
I12=int(int(I1(x,y),y,0,3),x,-3,3);
detaI1=(I11-I12)./(I11+I12)
vpa(detaI1)
I2(x,y)=-6.163e+04-1.826*x-786.1*y+845.8*x.^2+0.2181*x*y+714.3*y.^2;
I21=int(int(I2(x,y),y,-3,0),x,-3,3);
I22=int(int(I2(x,y),y,0,3),x,-3,3);
detaI2=(I21-I22)./(I21+I22)
vpa(detaI2)
I3(x,y)=-6.154e+04-1.861*x-799*y+844.6*x.^2+0.2141*x*y+713.1*y.^2;
I31=int(int(I3(x,y),y,-3,0),x,-3,3);
I32=int(int(I3(x,y),y,0,3),x,-3,3);
detaI3=(I31-I32)./(I31+I32)
vpa(detaI3)
I4(x,y)=-6.146e+04-1.931*x-806.6*y+843.4*x.^2+0.2102*x*y+712*y.^2;
I41=int(int(I4(x,y),y,-3,0),x,-3,3);
I42=int(int(I4(x,y),y,0,3),x,-3,3);
detaI4=(I41-I42)./(I41+I42)
vpa(detaI4)