怎样输出这样一个运用ODE45解的偏微分方程图像
 

方程组为：A1=x(1)+ix(2);
A2=x(3)+ix(4); *是一个复数*

m文件为：
function yx=qwe(t,x)
yx=zeros(4,1);
yx(1)=-x(1)-5.*x(2)+x(1).*x(3)+x(2).*x(4);
yx(2)=-x(2)+5.*x(1)+x(1).*x(4)-x(2).*x(3);
yx(3)=8-x(3)+x(4)-x(1).*x(1)+x(2).*x(2);
yx(4)=-x(4)-x(3)-2.*x(1).*x(2);

我的程序为：
tspan=0:0.01:500;
x0=[1 1 1 1];
[t,x]=ode45('qwe',tspan,x0);
plot(x(18000:27000,2),x(18000:27000,1))

上面这个方程输出的是以横纵坐标分别为x(2)和x(1)所输出的图像
我想请问如果想输出横坐标为x(2)^2+x(1)^2纵坐标为x(3)^2+x(4)^2的函数图像，也就是横纵坐标为A1模的平方和A2模的平方，应该做如何的改动？
请大家帮忙，谢谢大家！