clear;
clear all;
clc;
A=rand(3000,3000);
B=rand(3000,1);

tic;
C=A*B;
time1=toc;
fprintf('CPU Time = %3.3f sec\n', time1);

Ag=gpuArray(A);
tic;
Bg=gpuArray(B);
Cg=Ag*Bg;
C=gather(Cg);
time2=toc;

fprintf('GPU Time = %3.3f sec\n', time2);