btstrp_4.10 = function(original_data){

fit_original = lm(Emp~GNP+Def+Unemp+Mil+Pop, original_data);
print("Fitted summary of original trail");
print(summary(fit_original));

coef1 = NULL;

for (i in 1:1000){
	GNP = original_data$GNP + runif(16, -500,500);
	Pop = original_data$Pop + runif(16, -500,500);
	Def_1 = original_data$Def[1:7]+runif(7,-0.5,0.5);
	Def_2 = original_data$Def[8:16]+runif(9,-5,5);
	Def = c(Def_1, Def_2);
	Unemp = original_data$Unemp+runif(16,-5,5);
	Mil = original_data$Mil+runif(16,-5,5);
	Emp = original_data$Emp+runif(16,-50,50);

	new_Data = data.frame(cbind(GNP, Pop, Def, Unemp, Mil, Emp));
	temp_fit = lm(Emp~GNP+Def+Unemp+Mil+Pop, new_Data);
	coef1 = cbind(coef1, temp_fit$coef);
	}

print("Standard deviation of new data:");
apply(coef1, 1,sd);
}