stimulate.m

有限差分

syms a b c d E f g h I J T t h1 h2;
 jn2=E*(a-b)*(c*exp(a-b)-d)/(h2*(exp(a-b)-1));
 jn1=g*(a-b)*(c*exp(a-b)-d)/(h1*(exp(a-b)-1));
 R=((d*I)-1)/(T*(I+1)+t*(d+1));
jp2=E*(a-b)*(I*exp(a-b)-J)/(h2*(exp(a-b)-1));
jp1=g*(a-b)*(I*exp(a-b)-J)/(h1*(exp(a-b)-1));
h=(h1+h2)/2;
Fn=jn2-jn1-h*R
Fp=jp2-jp1+h*R