If I understand what you want to do, this code can help you.
df1 = data.frame(replicate(30,rnorm(40,20,5)))
colnames(df1) = paste0("column_",1:30)
df1
y = 1:40
# attach(df1); lm(y ~ column_1 + column_2 + .... + column_30)
lm(y ~ ., data = cbind(y,df1))
df2 = matrix(NA,40,30)
for(i in 1:30){df2[,i] = df1[,i]*seq(1.1,4,0.1)[i]}
colnames(df2) = paste0("new_column_",1:30)
df = cbind(df1,df2)
df[1:5,c(1:4,56:60)]
       column_1 column_2 column_3  column_4 new_column_26 new_column_27 new_column_28 new_column_29 new_column_30
1 21.497199 12.40726 14.02247 22.410877      83.30869      81.12982      64.32517      36.74374      55.54547
2 19.976892 28.61573 26.81687  9.821782      86.66725      72.59295     104.11971      83.24529      91.68181
3  9.333362 15.95690 15.14520 11.023078      32.61563      57.83005      89.69093      68.25862      63.83636
4 16.442331 13.04171 15.54477 18.369973      85.30285     100.67551      95.59633      87.21039      89.42617
5 17.757539 19.24030 25.93277 20.191389      71.99046      68.09845      91.93291     105.11022      65.71628