想实现高斯分布，请看这是怎么回事。该如何解决
 

本人是新手，想实现高斯分布的图，可是N=1时，应该会出现高斯分布，可是只出现最大值到最小值的一半的图形，望高手指教

clc,clear,close all;
global N
t0=0;tf=200;
y1=0.66e-6;%设φm=0.66，φ（0，0）=φm*e-4
y2=0.0001;%设初始A=0

%options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[T,Y]=ode45('lie',[t0,tf],[y1,y2]);

plot(T,Y(:,1),'r-')


lie.m文件的内容如下

function dy=lie(t,y)
N=2;
dy=zeros(2,1);
dy(1)=y(1)*N*exp(-y(2))-y(1);
dy(2)=y(1);