c=[0 0 0 1 1 0 0; 0 0 1 1 1 0 0; 1 1 1 1 1 0 0; 0 1 1 1 1 0 0; 0 0 0 0 1 0 0];
spy(c)

clc, clear all
a=0.2e-9:0.2e-9:4e-9;
b=0.2e-9:0.2e-9:4e-9;
for l1=1:length(a)
     for l2=1:length(b)
    c1=a(l1);
    c2=b(l2);
    sim('ModelAnton')
    if t.time(end)<1;
        q(l1,l2)=1;
    else
        q(l1,l2)=0;
    end
     end
end
[x y]=meshgrid(a,b);
plot3(x,y,q,'x'),grid on
axis([a(1)-1e-9 a(end)+1e-9 b(1)-1e-9 b(end)+1e-9 -1 2])
figure;
hA=axes;
spy(q)
set(hA, 'XTicklabel', [0:0.2e-9:4e-9])
set(hA, 'YTicklabel', [0:0.2e-9:4e-9])

function reshenie
clear all;clc;close all;
an=2.2e-9; ak=6e-9; ah=0.2e-9;
bn=0.2e-9; bk=4e-9; bh=0.2e-9;
a=an:ah:ak;
b=bn:bh:bk;
for l1=1:length(a)
     for l2=1:length(b)
    c1=a(l1);
    c2=b(l2);
    %sim('ModelAnton')
    z=rand(1);
    if z>0.5;
        q(l1,l2)=1;
    else
        q(l1,l2)=0;
    end
     end
end
%[x y]=meshgrid(a,b);
%plot3(x,y,q,'x'),grid on
%axis([a(1)-1e-9 a(end)+1e-9 b(1)-1e-9 b(end)+1e-9 -1 2])
%figure;
hA=axes;
spy(q)
 
[nY,nX]=size(q);
XL=get(hA,'XTicklabel');
nXL=length(XL);
kX=round(nX/(nXL-1));
set(hA, 'XTickMode', 'manual', 'XTicklabel', bn-bh:kX*bh:bk)
YL=get(hA,'YTicklabel');
nYL=length(YL);
kY=round(nY/(nYL-1));
set(hA, 'YTickMode', 'manual', 'YTicklabel', an-ah:kY*ah:ak)

set(hA, 'YDir', 'normal');