Question: Supernovae mesh plot

Hey Guys, 

I have plotted the following graph, its the position of an elementary volume before a supernovae takes place.

x₀=r₀sinθ0cosΦ₀

y0=r₀sinθ₀sinΦ₀

z₀=r₀cosθ₀

Created a mesh of 10 mesh points, divided the ranges  

0<r<1

0<θ<∏

0<Φ<2∏

up into 10 equal parts, and taken r₀, θ₀, Φ₀ as iterations through these values.

I then used the above euations to find x₀, y₀,z₀, and I then made a 3D plot of these values

the code used was

plots[pointplot3d]([seq(seq(seq([r,i*Pi,j*Pi],r=0..1,0.1),i=0..2,1/5),j0..1,1/10)],coords=spherical,symbol=diamond,shading=zhue);

**********************************************************************************

Now, I need to obtain 

x=rsinθ₀cosΦ₀

y=rsinθ₀sinΦ₀

z=rcosθ₀

Where I am to use r instead of r₀

r is a function of time let

r=4t

I am trying to simulate how the supernova expands over time. 

much appreciated