Нужно пояснить строки, которые подчеркнул

clc; clear all;
N=8;
load front.mat
r=0:0.001:1;
p=0:0.005:6.25;
[P,R]= meshgrid(p,r);
Z=zeros(length(r),length(p),N);
Q=zeros(1,N);Q1=zeros(1,N);
[x,y]=meshgrid(0:0.01:1,0:0.05:2*pi);
 
global ZZZ1
G=interp2(x,y,A,R,P);
for k=1:N
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%    
      n=floor((sqrt(8*k+1)-1)/2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%    
m=2*k-n^2-2*n;
   if m==0
       E=2; 
   else E=1;
   end;
   Z=Zern(x,y,k);
   Z2=interp2(x,y,Z,R,P);
   ZZZ1=Z2.*G.*R;
   Q(k)=dblquad(@zh,0,6.25,0,1);
   Q1(k)=Q(k)*((2*n+2)/(E*pi));
end;
disp (Q1)
 
function [zh]=zh(p,r)
zh=zeros(1,length(p));
global ZZZ1
for i=1:length(p)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%    
      zh(i)=ZZZ1(round(r./0.001+1),round(p(i)./0.005+1)); 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;