ИЭ / 4 семестр / Курсач от Кочетковой / Графики сфера в матлабе / UE23

clc, close all
R1=0.01;
R121=0.025;
R2=0.05;
r11=R1:0.0001:R121;
r21=R121:0.0001:R2;
U11=-23540489.64*r11.^2-317.84./r11+44138.42; %введи своё U1
U21=835.59./r21-16711.86; %введи своё U2
E11=47080979.28*r11 -317.84./((r11).^2); %введи своё E1
E21=835.59./((r21).^2); %введи своё E2


R122=0.012;
r12=R1:0.0001:R122;
r22=R122:0.0001:R2;
U12=-23540489.64*r12.^2 -9.35./r12+13289.53; %введи своё U1
U22=144.002./r22-2880.04; %введи своё U2
E12=47080979.28*r12 -9.35./((r12).^2); %введи своё E1
E22=144.002./((r22).^2); %введи своё E2

R123=0.04;
r13=R1:0.0001:R123;
r23=R123:0.0001:R2;
U13=-23540489.64*r13.^2-652.27./r13+77580.59; %введи своё U1
U23=4721.83./r23-94436.69; %введи своё U2
E13=47080979.28*r13 -652.27./((r13).^2); %введи своё E1
E23=4721.83./((r23).^2); %введи своё E2


figure;
subplot(2,1,1);
hold on;
set(gca,'box','on');
grid on;
plot(r11,U11,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot(r21,U21,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot([R121, R121], [0, U11(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R121+0.001, 8000, 'R_{12}=0.025 м', 'FontSize',14);

plot(r12,U12,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot(r22,U22,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot([R122, R122], [0, U12(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R122+0.0005, 2000, 'R_{12}=0.012 м', 'FontSize',14);

plot(r13,U13,'--','color',[185 185 185]/255,'LineWidth',1.3);
plot(r23,U23,'--','color',[185 185 185]/255,'LineWidth',1.3);
plot([R123, R123], [0, U13(end)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R123+0.0005, 10000, 'R_{12}=0.04 м', 'FontSize',14);

xlim([0.01 0.05])
ylim([0 40000]);
xlabel('r, м','fontsize',14);
ylabel('U, В','fontsize',14,'Rotation',0,'Position', [0.008 25000 0]);
subplot(2,1,2);
hold on;
set(gca,'box','on');
grid on;
plot(r11,E11,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot(r21,E21,'-','color',[0 0 0]/255,'LineWidth',1.3);
plot([R121, R121], [-6500000, E21(1)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R121+0.001, -3000000, 'R_{12}=0.025 м', 'FontSize',14);

plot(r12,E12,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot(r22,E22,'-','color',[100 100 100]/255,'LineWidth',1.3);
plot([R122, R122], [-6500000, E22(1)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R122+0.0005, -5500000, 'R_{12}=0.012 м', 'FontSize',14);

plot(r13,E13,'--','color',[185 185 185]/255,'LineWidth',1.3);
plot(r23,E23,'--','color',[185 185 185]/255,'LineWidth',1.3);
plot([R123, R123], [-6500000, E23(1)],'--','color',[0.5 0.5 0.5],'LineWidth',1.3);
text(R123+0.0005, -3000000, 'R_{12}=0.04 м', 'FontSize',14);

xlim([0.01 0.05])
ylim([-6500000 3000000]);
xlabel('r, м','fontsize',14);
ylabel('E, В/м','fontsize',14,'Rotation',0,'Position', [0.008 0 0]);