This is done using Student:-NumericalAnalysis:-Euler as follows
ode:=diff(y(x),x)=2*x*y(x)+x; IC:=y(0)=1; ending_x:=2; Student:-NumericalAnalysis:-Euler(ode,IC,x=ending_x, 'numsteps'=41, 'output'='plot', 'plotoptions'=['gridlines']);
To obtain table of values instead do this
ode:=diff(y(x),x)=2*x*y(x)+x; IC:=y(0)=1; ending_x:=2; Student:-NumericalAnalysis:-Euler(ode,IC,x=ending_x, 'numsteps'=41, 'output'='information');