matlab设置不等分刻度,坐标轴单位怎么改成不等分的间距

最近根据别人的文献在仿真图形，别人的图是这个样子的我自己仿真的图形是这样的

想像别人那样把纵坐标单位换成10^13到10^19,请问各位应如何操作？

自己检查了下，感觉自己的图形应该没错，就是不知道怎么让图像对方那样

程序附下：

clc;clear all;

global T;

yy_1=[];

yy_2=[];

%对变量作定义

h = waitbar(0,'Please wait...');

ti = 0;

for T=0.1*10^4:0.01*10^4:3*10^4

ti = ti + 1;

syms x1 y1 z1 xt y21 zt w22;

P=101.325*10^3;

k=1.3806503*10^(-23);

Me=9.10938188*10^(-31);

h=6.62606876*10^(-34);

O_1=1313900;

N_1=1402300;

O_2=3388200;

N_2=2856100;

N_A=6.02*10^23;

%只考虑第一次电离时，x1是中性粒子数密度,y1是离子数密度，z1是电子数密度

sol=solve(P-(x1+y1+z1)*k*T,...

y1-z1,...

y1*z1/x1-2*(2*pi*k*T*Me/h^2)^(3/2)*exp(-(N_1/N_A)/(k*T))*0.78084-2*(2*pi*k*T*Me/h^2