testbidiff.cpp itpp/itbase.h ../bdm/stat/libBM.h pmsm.h work/git/mixpp/pmsm/testbidiff.cpp ../bdm/stat/libBM.h pmsm.h stat/libFN.h itpp/itbase.h int int main () main #include<itpp/itbase.h> #include"../bdm/stat/libBM.h" #include"pmsm.h" usingnamespaceitpp; //Theselinesareneededforuseofcoutandendl usingstd::cout; usingstd::endl; intmain() { vecx0="1001201002"; vecu0="11"; doubleh=1e-6; IMpmsmI; I.set_parameters(0.28,0.003465,h,0.1989,1.5,4.0,0.04,0.0); vecx(x0);x(0)+=h; cout<<(I.eval(x,u0)-I.eval(x0,u0))/h<<endl; x=x0;x(1)+=h; cout<<(I.eval(x,u0)-I.eval(x0,u0))/h<<endl; x=x0;x(2)+=h; cout<<(I.eval(x,u0)-I.eval(x0,u0))/h<<endl; x=x0;x(3)+=h; cout<<(I.eval(x,u0)-I.eval(x0,u0))/h<<endl; cout<<endl; matA(4,4); I.dfdx_cond(x0,u0,A); cout<<A.T()<<endl; //Exitprogram: return0; }