# \[Help\] Plot/export of solutions from PDE

**URL:** https://discourse.julialang.org/t/help-plot-export-of-solutions-from-pde/50020
**Category:** New to Julia
**Tags:** question, plotting
**Created:** [November 12, 2020, 3:06am UTC](https://discourse.julialang.org/t/help-plot-export-of-solutions-from-pde/50020 "2020-11-12T03:06:29Z")
**Posts on this page:** 1
**Page:** 1

<div class="post-metadata">

### Author: ![jinkel7](https://avatars.discourse-cdn.com/v4/letter/j/e36b37/32.png) [@jinkel7](https://discourse.julialang.org/u/jinkel7)
#### Post date: [November 12, 2020, 3:06am UTC](https://discourse.julialang.org/t/help-plot-export-of-solutions-from-pde/50020/1 "2020-11-12T03:06:30Z")

</div>

Hello, I’m quite new to Julia.  
I’m having a problem with “plot” and “exporting”.  
Can anyone give help for a newbie?

Using VoronoiFVM for one-dimensional PDE, the solution is expressed as Array{Float64,1}.

I want to express it as a graph with time as e.g. (x, solution[2,:], t)  
Like a “surface” graph.

This problem also occurs for the exporting solution.  
The exported csv file only shows the solution for the final time period.

Code below without the equations

```julia
using VoronoiFVM
using ExtendableGrids
using PyPlot
using Plots
Length=0.0075;
edge_n=100; # number of edges
X1=collect(0.0:Length/edge_n:Length)
grid=ExtendableGrids.simplexgrid(X1)
ExtendableGrids.plot(grid,Plotter=PyPlot)

using DelimitedFiles
let
t=0.0; del_t=60.0
@gif for i=1:600
         t=t+del_t      
        solve!(solution,inival,system,tstep=del_t)

        writedlm("results.csv", solution[2,:])
        
        @views begin
        P1=heatmap(reshape(solution[1,:],length(X1),1),colorbar=:right,color=:viridis,yflip=true,clim=(230,260))
        P2=heatmap(reshape(solution[2,:],length(X1),1),colorbar=:right,color=:viridis,yflip=false,clim=(230,260))
        P3=heatmap(reshape(solution[3,:],length(X1),1),colorbar=:right,color=:viridis,yflip=true)
        P4=Plots.plot(solution[2,:],ylim=(230,260), yflip = false)
        P=Plots.plot(P1,P2,P3,P4,layout=(1,4))
        end
        inival.=solution; del_t*=1.00
   end
end

```

![tmp](https://global.discourse-cdn.com/julialang/original/3X/a/6/a6dc243c1dc17dbf0dcf01faa2404d11932e1e79.gif)

This is the current graph that time flows as “gif”. P2 and P4 are what I am focusing on.  
Can it be converted to (x, solution, t) graph or exported to csv file?