data b(drop=x_1 x_2);
retain x_1 x_2;
set a;
output;
if  n in (4,7) then
  do;
    n=n-1;
    x=(x+x_2)/2;
    output;
  end;
x_1=x;
x_2=x_1;
run;
proc sort;
by n;
run;
proc print;
run;

楼上高人，学习

Thank a lot! What is the purpose of creating x_2? Why can't I simply use x=(x+x_1)/2?

pls run the codes below, then you can find your results.

data b;
retain x_1 x_2;
set a;
output;
if  n in (4,7) then
  do;
    n=n-1;
    x=(x+x_2)/2;
    output;
  end;
x_1=x;
x_2=x_1;
run;

proc print;
run;

I use the same method for data with different gaps, but the result is not exactly the same as what I expect.
For example,

data a;                                                                                                                                   
input   n   x;                                                                                                                           
cards;                                                                                                                                    
1   10                                                                                                                                    
2   11.5                                                                                                                                 
4   12                                                                                                                                    
5   10.6                                                                                                                                 
7   9.8                                                                                                                                   
10   10.1                                                                                                                                 
11  12                                                                                                                                    
;                                                                                                                                         
run;                                                                                                                                      

data b;                                                                                                                                   
retain x_1 x_2;                                                                                                                           
set a;                                                                                                                                    
output;                                                                                                                                   
if  _n_ in (3,5) then                                                                                                                     
   do;                                                                                                                                    
    n=n-1;                                                                                                                                
    x=(x+x_2)/2;                                                                                                                          
    output;                                                                                                                              
  end;                                                                                                                                    
if  _n_ in (6) then                                                                                                                       
   do j=1 to 2;                                                                                                                           
    n=n-1;                                                                                                                                
    x=(x_2)*j/3 + x*(3-j)/3;                                                                                                                     
    output;                                                                                                                              
  end;                                                                                                                                    

x_1=x;                                                                                                                                    
x_2=x_1;                                                                                                                                 
run;                                                                                                                                      

proc sort data=b;                                                                                                                        
by n;                                                                                                                                    
run;                                                                                                                                      
proc print;                                                                                                                              
run;

The result is as follow:


It is different from the expected result for n= 8, 9

n      x
8      9.9
9      10