This question is a follow-up to [my previous question][1]. The code I use to plot is the following: Data = NIntegrate[myrotorz1[x, y, t], {x, -5, 5}, {y, -5, 5}]; Table[Plot[Re[Data], {t, 0, 6, 0.1}], PlotRange -> All] I would like to know whether to use the `NIntegrate` within the `NDSolve` domain. I tryed with Plot[Re[NIntegrate[myrotorz1[x, y, t], {x, -5, 5}, {y, -5, 5}]], {t, 0, 6}] got this error NIntegrate::inumr: The integrand myrotorz1[x,y,0.000122571] has evaluated to non-numerical values for all sampling points in the region with boundaries {{-5,5},{-5,5}}. >> NIntegrate::inumr: The integrand myrotorz1[x,y,0.000122571] has evaluated to non-numerical values for all sampling points in the region with boundaries {{-5.,5.},{-5.,5.}}. >> NIntegrate::inumr: The integrand myrotorz1[x,y,0.000122571] has evaluated to non-numerical values for all sampling points in the region with boundaries {{-5.,5.},{-5.,5.}}. >> General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >> [1]: http://mathematica.stackexchange.com/questions/96419/how-to-do-vectorplot