%% Material parameter and input Matrix
P_source(material==1)=Material_mat(1,2);
P_source(material==2)=Material_mat(2,2);
P_source(material==3)=Material_mat(3,2);
P_source(material==4)=Material_mat(4,2);
P_source(material==5)=Material_mat(5,2);
P_source(material==6)=Material_mat(6,2);
P_source(material==7)=Material_mat(7,2);

Dip_angle=random('Normal',45,8,nx,ny); 

kc0(material==1)=Material_mat(1,4);sig_star(material==1)=Material_mat(1,5);
kc0(material==2)=Material_mat(2,4);sig_star(material==2)=Material_mat(2,5);
kc0(material==3)=Material_mat(3,4);sig_star(material==3)=Material_mat(3,5);
kc0(material==4)=Material_mat(4,4);sig_star(material==4)=Material_mat(4,5);
kc0(material==5)=Material_mat(5,4);sig_star(material==5)=Material_mat(5,5);
kc0(material==6)=Material_mat(6,4);sig_star(material==6)=Material_mat(6,5);
kc0(material==7)=Material_mat(7,4);sig_star(material==7)=Material_mat(7,5);