clear, clc;
%% build matrice
m1=[0.8+0.1i,0.5i;0.3i,0.3+0.4i];
m2=[0.1+0.05i,0.2i;
0,0.4+0.9i];
S=[1 0.6+0.4i];
%% differentiate symbolic M,
syms kx ky
M = m1 + m2*exp((kx -ky)*1i) + m2'*exp(-(kx-ky)*1i);
dmdkx = diff(M,kx); % 符号函数偏微分
%% get function F
F = S*dmdkx*S';
global ki; ki=0;
I =dblquad('f', -1,1,-2,2)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% function 单独放在一个文件(和Untitle文件同目录)
function F=f(kx,ky)
global ki;
ki=ki+1;
F= exp(- kx.*1i + ky.*1i)*(- 319/500 - (97*1i)/250) + exp(kx.*1i - ky.*1i)*(- 319/500 + (97*1i)/250);