各位大神，请问如何绘制以下曲线，本人是新手。求指教。下面是我做的程序。
alfa=[5 10 15 25 27 30];
i=1;
L=62;
Rd=73.5;
Rc=69.5;
i=1;
j=1;
for i=1:1:6
for j=1:1:720
theta[i][j]=pi*j/360;
end
end
for i=1:1:6
for j=1:1:180
Rdc[i][j]=Rd*sqrt(sin(theta[i][j]^2+(cos(theta(j))*cos(alfa(i)))^2);
delta[i][j]=asin(Rd*sin(theta[i][j])/Rdc-theta[i][j]);
fi[i][j]=asin(Rdc[i][j]*sin(delta[i][j])/sqrt((Rc-Rd*cos(delta[i][j]))^2+(Rdc[i][j]*sin(delta[i][j]))^2));
beta[i][j]=asin(sqrt(Rc^2+Rdc[i][j]^2-2*Rc*Rdc[i][j]*cos(delta[i][j]))/L);
end
for j=181:1:360
Rdc[i][j]=Rd*sqrt(sin(theta[i][j])^2+(cos(theta[i][j])*cos(alfa(i)))^2);
delta[i][j]=-asin(Rd*sin(theta[i][j])/Rdc[i][j]-theta[i][j])-theta[i][j]+pi;
fi[i][j]=asin(Rdc[i][j]*sin(delta[i][j])/sqrt((Rc-Rd*cos(delta[i][j]))^2+(Rdc[i][j]*sin(delta[i][j]))^2));
beta[i][j]=asin(sqrt(Rc^2+Rdc[i][j]^2-2*Rc*Rdc[i][j]*cos(delta[i][j]))/L);
end
for j=361:1:540
Rdc[i][j]=Rd*sqrt(sin(theta[i][j])^2+(cos(theta[i][j])*cos(alfa(i)))^2);
delta[i][j]=-asin(Rd*sin(theta[i][j])/Rdc[i][j]-theta[i][j])-theta[i][j]+pi;
fi[i][j]=asin(Rdc[i][j]*sin(delta[i][j])/sqrt((Rc-Rd*cos(delta[i][j]))^2+(Rdc[i][j]*sin(delta[i][j]))^2));
beta[i][j]=asin(sqrt(Rc^2+Rdc[i][j]^2-2*Rc*Rdc[i][j]*cos(delta[i][j]))/L);
end
for j=541:1:720
Rdc[i][j]=Rd*sqrt(sin(theta[i][j])^2+(cos(theta[i][j])*cos(alfa(i)))^2);
delta[i][j]=asin(Rd*sin(theta[i][j])/Rdc-theta[i][j])-theta[i][j]+2*pi;
fi[i][j]=asin(Rdc[i][j]*sin(delta[i][j])/sqrt((Rc-Rd*cos(delta[i][j]))^2+(Rdc[i][j]*sin(delta[i][j]))^2));
beta[i][j]=asin(sqrt(Rc^2+Rdc[i][j]^2-2*Rc*Rdc[i][j]*cos(delta[i][j]))/L);
end
end
for i=1:1:6
for j=1:1:720
plot(theta[i][j]*360/(2*pi),delta[i][j],theta[i][j]*360/(2*pi),fi[i][j],theta[i][j]*360/(2*pi),beta[i][j]);
hold on;
end
end