Constant E=2.0*10^11,I=22,l=1,d=22;
loopconstant r1=[1000000:200000:200000001],r2=[1000000:200000:200000001];
function t1-(1/r1);
t2-(1/r2);
SubDivision;
function (t1^2*(k^2*E*I)^2+(k^2*E*I)^2*t2^2-2*(k^2*E*I)^2*t1*t2+k^4*l^2+2*(k^2*E*I)^2*k^2*l^2*t1*t2+t1^2*t2*2*(k^2*E*I)^3*l+2*k^2*E*I*k^2*l*t1-2*(k^2*E*I)^3*l*t1*t2^2-2*k^2*E*I*k^2*l*t2+t2^2*(k^2*E*I)^2*k^2*l^2+t1^2*(k^2*E*I)^2*k^2*l^2-2*(k^2*E*I)^2*k^2*t1*t2*l^2+4*k^2+4*k^2*k^2*E*I*l*t2-4*k^2*k^2*E*I*l*t1+(k^2*E*I)^4*l^2*t1^2*t2^2)^(1/2)*sin(k*l+arctan((t2*k*k^2*E*I*l-t1*k^2*E*I*k*l+2*k)/(k^2*E*I*t1-k^2*E*I*t2+(k^2*E*I)^2*l*t1*t2+k^2*l)))-2*k;
function p-k*k*E*I;
我结果是想把P算出来的，但结果却只算到了t1，t2。请问大佬们我哪里出错了啊，万分感激，结果如下：该怎么改代码呢
！