Как построить график следующей функции

clc;
clear all;
%дискретная случайная величина лямбда
lambda1 = 1;
lambda2 = 2;
lambda3 = 3;
lambda4 = 4;
lambda = [lambda1 lambda2 lambda3 lambda4];
%вероятности появления лямбда
p1 = 0.2;
p2 = 0.3;
p3 = 0.15;
p4 = 0.35;
p = [p1 p2 p3 p4];
%белый шум с 0 мат.ожиданием mu и ско sigma
mu=0;
sigma=0.1;
N = 10;
shum = normrnd(mu,sigma,1,N);
%график вероятностей лямбда
stem(lambda,p);
 
%генератор дискретной величины с заданой вероятностью
x = rand(1);
if x <=p(1) & x>0
    l = lambda(1);
elseif x <=p(2)+p(1) & x>p(1)
    l = lambda(2);
elseif x <=p(3)+p(2)+p(1) & x>p(2)+p(1)
    l = lambda(3);
elseif x<=p(1)+p(2)+p(3)+p(4) & x>p(1)+p(2)+p(3)
    l = lambda(4);
end
 
%наблюдение
z = l + shum;
t = 0:N-1;
figure(2);
plot(t,z);
 
%WZn
W = zeros(1,N); 
for n = 1:1:N 
for i = 1:1:4 
W(n) = W(n) + p(i)*((1/(sqrt(2*pi)*sigma^2))^N)*exp(-(1/(2*sigma^2))*((z(N)-lambda(i)).^2)); 
end 
end

clc, clear all, close all
 
p = [0.25 0.25 0.25 0.25];
lambda_i = 1:4;
lambda = 5;
delta = 0.1;
 
W = @(lambda)(sum(p.*delta.*(lambda - lambda_i)));

>> W(lambda)
 
ans =
 
    0.2500

clc, clear all, close all
 
p = [0.25 0.25 0.25 0.25];
lambda_i = 1:4;
n = 1:4;
sigma = 10;
 
W = @(z_n)(sum(p.*(1/(sigma.*sqrt(2*pi))).^n.*exp((-1)/(2.*sigma.^2)*sum(z_n - lambda_i).^2)));

clc, clear all, close all
 
p = [0.25 0.25 0.25 0.25];
lambda_i = 1:4;
n = 1:4;
sigma = 10;
 
z_n = 3;
 
for i = 1:4
    W(i) = (sum(p(i)*(1/(sigma.*sqrt(2*pi))).^n(i).*exp((-1)/(2.*sigma.^2)*sum(z_n - lambda_i(i)).^2)));
end
 
WW = sum(W);
 
disp(WW)

clc, clear all, close all
 
p = [0.25 0.25 0.25 0.25];
lambda_i = 1:4;
n = 1:4;
sigma = 10;
 
z_n = rand(1,4)';
 
for i = 1:length(z_n)
    W(i,:) = (sum(p.*(1/(sigma.*sqrt(2*pi))).^n.*exp((-1)/(2.*sigma.^2)*sum(z_n(i) - lambda_i).^2)));
end
 
disp([z_n, W])