% experiment.m
% an experimental script for radial model data
% generation and analysis.
%
% coordinates of event window
A = 0;
B = 10;
C = 0;
E = 10;
% coordinates of center
X = 5;
Y = 5;
% model parameters .. user must set Nsoln to correct function
M = 1000;
D = 1;
a = 0.0;
% discretization parameters
dr = 0.25;
T = 16;
dt = 0.25;
% event detection probability
P = 1;
% profile time
t0 = 2;
% event placement
RadialDist = 1;
% repetitions for average model
count = 25;
% for predictability...unrem the next line
rand('seed',2);
% generate the events
events = RadSim(A,B,C,E,X,Y,M,D,a,dr,T,dt,P,RadialDist);
disp('Size of event matrix:');
disp(size(events));
% do radial plot of events with fronts
plotevents(A,B,C,E,X,Y,D,a,dr,T,dt,events);
% plot radial profiles of data and exact solution at select times
plotprofiles(A,B,C,E,X,Y,D,a,dr,T,dt,events);
% get profiles of averages of data at time t0 and selected probabilities
plotdataprofiles(A,B,C,E,X,Y,M,D,a,dr,T,dt,RadialDist,t0,count)
% compute error norms for averages of data and 
errors = errornorms(A,B,C,E,X,Y,M,D,a,dr,T,dt,RadialDist,count)