diff --git a/plot4.R b/plot4.R
new file mode 100644
index 00000000000..9a3e9b68ffc
--- /dev/null
+++ b/plot4.R
@@ -0,0 +1,17 @@
+temp <- tempfile()
+download.file('http://d396qusza40orc.cloudfront.net/exdata%2Fdata%2Fhousehold_power_consumption.zip',temp)
+all_data<-read.csv(unzip(temp),sep = ";",header = T,quote = "",na.strings='?')
+unlink(temp)
+my_data<-subset(all_data,Date=='1/2/2007' | Date=='2/2/2007')
+my_data$dt<-strptime(paste(my_data$Date,my_data$Time),"%d/%m/%Y %T")
+png('plot4.png',bg = "transparent")
+par(mfrow = c(2,2) )
+plot(my_data$dt,my_data$Global_active_power,ylab='Global Active Power',xlab='',t='l')
+plot(my_data$dt,my_data$Voltage,ylab='Voltage',xlab='datetime',t='l')
+plot(my_data$dt,my_data$Sub_metering_1,ylab='Energy Submetering',xlab='',t='n')
+lines(my_data$dt,my_data$Sub_metering_1,col='black')
+lines(my_data$dt,my_data$Sub_metering_2,col='red')
+lines(my_data$dt,my_data$Sub_metering_3,col='blue')
+legend('topright',col=c("black","red","blue"),legend = c("Sub_metering_1","Sub_metering_2","Sub_metering_3"),lwd=1,bty ='n' )
+plot(my_data$dt,my_data$Global_reactive_power,ylab='Global_reactive_power',xlab='datetime',t='l')
+dev.off()
\ No newline at end of file
diff --git a/plot4.png b/plot4.png
new file mode 100644
index 00000000000..3090a948a76
Binary files /dev/null and b/plot4.png differ