%plot light curves for ZTF data

    n=0;
    infile=['lc_' num2str(n) '.dat'];
    A=readmatrix(infile,'CommentStyle','#','Delimiter',' ');
    HJD=A(:,1);
    flux=A(:,2);
    eflux=A(:,3);
    clear A
    plot(HJD,flux,'.')
    ylim([0 2]*1E-4)