> a<- 50  
> p1<-20
> p2<-30  
> c1<-5  
> c2<-6  
> em<-3  
> en<-2  
> e1<-1  
> e2<-2
> d<-12  
> pc<-1  
> θ<-0.5   
> E<-1000  
> r<-0.01
> B<-seq(0,8000,len = 65)
> u1<-a-p1+θ*(p2-p1)
> u2<-a-p2+θ*(p1-p2)
> y11<-{(B-c2*(u2+d*((1+θ)*e2-θ*e1)))*p2*c1*((1+θ)*e1-θ*e2)+p1*c2^2*((1+θ)*e2-θ*e1)*(u1+d*((1+θ)*e1-θ*e2))}/(p1*c2^2*((1+θ)*e2-θ*e1)+p2*c1^2*((1+θ)*e1-θ*e2))
> y21<-{B-c2*(u2+d*((1+θ)*e2-θ*e1))-c1*(u1+d*((1+θ)*e1-θ*e2))}*p1*c2*((1+θ)*e2-θ*e1)/{p1*c2^2*((1+θ)*e2-θ*e1)+p2*c1^2*((1+θ)*e1-θ*e2)}
> y12<-d*((1+θ)*e1-θ*e2)*(p1-(1+r)*c1)/p1+u1
> y22<-d*((1+θ)*e2-θ*e1)*(p2-(1+r)*c2)/p2+u2
> y13<-{(B+pc*E-(c2+pc*en)*(u2+d*((1+θ)*e2-θ*e1)))*p2*(c1+pc*em)*((1+θ)*e1-θ*e2)+p1*(c2+pc*en)^2*((1+θ)*e2-θ*e1)*(u1+d*((1+θ)*e1-θ*e2))}/(p1*(c2+pc*en)^2*((1+θ)*e2-θ*e1)+p2*(c1+pc*em)^2*((1+θ)*e1-θ*e2))
> y23<-(B+pc*E-(c2+pc*en)*(u2+d*((1+θ)*e2-θ*e1))-(c1+pc*em)(u1+d*((1+θ)*e1-θ*e2))）*p1*(c2+pc*em)*((1+θ)*e2-θ*e1)/(p1*(c2+pc*en)^2*((1+θ)*e2-θ*e1)+p2*(c1+pc*em)^2*((1+θ)*e1-θ*e2))
> y14<-d*((1+θ)*e1-θ*e2)*(p1-(1+r)*(c1+pc*em))/p1+u1
> y24<-d*((1+θ)*e2-θ*e1)*(p2-(1+r)*(c2+pc*en))/p2+u2              
> plot(B,y11,type="l",col=1,lwd =2,lty=2,xlab="initial capital B",ylab="production quantities y1,y2")
> lines(B , y21,col = 1,lwd=3,lty=1)
> lines(B ,y12,col = 2 ,lwd =2,lty = 2)
Error in xy.coords(x, y) : 'x' and 'y' lengths differ
> lines(B , y22,col = 2,lwd=3,lty=1)
Error in xy.coords(x, y) : 'x' and 'y' lengths differ
> lines(B , y13,col = 3 ,lwd =2,lty = 2)
> lines(B , y23,col = 3,lwd=3,lty=1)
Error in xy.coords(x, y) : 'x' and 'y' lengths differ
> lines(B , y14,col = 6 ,lwd =2,lty = 2)
Error in xy.coords(x, y) : 'x' and 'y' lengths differ
> lines(B , y24,col = 6,lwd=3,lty=1)
Error in xy.coords(x, y) : 'x' and 'y' lengths differ
> legend(0, 300, c("y11", "y12", "y13","y14"), col = c(1,2,3, 6), lty = c(2,2,2,2), merge = TRUE, bg = "white")
> legend(6000, 100, c("y21", "y22", "y23","y24"), col = c(1,2,3,6), lty = c(1,1,1,1), merge = TRUE, bg = "white")

求大神解答，字母大小写没有错，符号用英文的也没错，到底是哪里出了问题。求帮助。