include("eff_Drop.jl")

P = 5; Δt = 0.001; D_i = 0.2; D_o = 0.2; BP=0.5; x0=0.5; 
T=2; perc=0.5; L=0; R=2; k="none";

## Solve diffusion equation in phase-separated media
A, u, x, b, Nb, t_perc, αi, αo= @time solve_D("3D","FRAP","Y","N",zeros(9),
P,D_i,D_o,0,Δt,L,R,BP,x0,T,perc,zeros(9),zeros(9), k)

##
k = 0.001
A, u1, x, b, Nb, t_perc, αi, αo= @time solve_D("3D","FRAP","Y","N",zeros(9),
P,D_i,D_o,0,Δt,L,R,BP,x0,T,perc,zeros(9),zeros(9), k)
##
# Write results to file
# writedlm( "x_solution.csv",  vcat(x', u), ',')

plot_model(x, u, Nb, T, Δt)
plot_model(x, u1, Nb, T, Δt)