求大神指导一下,应该怎么修改代码呢(红色字体标明)?万分感激!!!
#vz
F2vcz<-function(AB,vcz,lch)
{alpha<-theta(AB,e3)
F2vczRD<-abs(1/2*r*d2*(vcz-lch)^2*((sin(alpha))^3*cn90+picf)*l0)
RL<-abs(1/2*r*d2*(vcz-lch)^2*(sin(alpha))^2*cos(alpha)*cn90*l0)
if (det(rbind(AB,e1,e2))==0)
{if(vcz>=lch)
{F2vczRD=F2vczRD
F2vczRLx=0
F2vczRLy=0
}
else
{F2vczRD=-F2vczRD
F2vczRLx=0
F2vczRLy=0
}
}
else if (det(rbind(AB,e1,e3))==0)
{if (vcz>=lch)
{if (alpha>=0&&alpha<=pi/2)
{F2vczRD=F2vczRD
F2vczRLx=RL
F2vczRLy=0
}
else (alpha>pi/2&&alpha<=pi)
{F2vczRD=F2vczRD
F2vczRLx=-RL
F2vczRLy=0
}
}
else
{if (alpha>=0&&alpha<=pi/2)
{F2vczRD=-F2vczRD
F2vczRLx=-RL
F2vczRLy=0
}
else (alpha>pi/2&&alpha<=pi)
{F2vczRD=-F2vczRD
F2vczRLx=RL
F2vczRLy=0
}
}
}
Error in if (vcz>=lch) {: 需要TRUE/FALSE值的地方不可以用缺少值
|