MATLAB中采用ODE45求解方程组

@CSDN专家-Matlab_Fans 您好，我想请问如果我用MATLAB画出一个图是这样的，通过改变中心波长，是可以改变在X轴为1的地方Y轴的值的，但是我想把这个点的值取出来，然后写一个关于中心波长和耦合效率的图，请问应该怎么做呢？

img

%观察中心波长变化时，效率的变化，ODE45解方程组画图
clc;
clear;
L=1;%声光作用长度1m

w=1550.*10.^(-9);%中心波长 单位米
s1=0.0017032967;
%s2=7.561.*10.^(-4);
s2=6.2753.*10.^(-4);
E1=9.1.*10.^(-4);%
%n2=2.05.*10.^(-3);
E2=2.47.*10.^(-3);%A02
a1=((pi.*E1)./w)-(pi./s1);%相位失配因子带入方程
a2=((pi.*E2)./w)-(pi./s2);

K0=pi/(2.*L);%A21(02)
%K0=1.11;%A11
K1=sqrt(2).*K0;
K2=K1; %（可获得的A21(02)转换效率最高）
%K2=0; %(可获得的A11转换效率最高)K2=0
f = @(z,x)[i.*K1.*x(2).*exp(-2i.*a1.*z);i.*K1.*x(1).*exp(2i.*a1.*z)+i.*K2.*x(3).*exp(-2i.*a2.*z);i.*K2.*x(2).*exp(-2i.*a2.*z)];
[z,y] = ode45(f,[0,1],[1;0;0]);%A11
%[z,y] = ode45(f,[0,1],[1;0;0]);
A=abs(y(:,1));
B=abs(y(:,2));
C=abs(y(:,3));
plot(z,A.^2,'r');%A01 红色曲线
hold on;
plot(z,B.^2,'b');%A11 蓝色曲线
hold on;
plot(z,C.^2,'k--');%A21 (02)虚线
title('耦合效率随长度变化的情况');
xlabel('长度 L ');
ylabel('耦合效率');
hold off;