%{
Setup and execution script for tactoidSolnCubeMS
Edit and run this script to run tactoidSolnCubeMS
Make sure the following files are located in this
directory:

rTSCMS.m
equationsMS.m
jacobian2MS.m
getTacEnergyMS.m
tactoidSolnCubeMS.m
Lambda.mat
%}

M = ; %Number of grid spaces
m = ; %Topological charge
Rin = ; %Initial tactoid/defect size
Rout = ; %System size 
kappa2 = ; %First anisotropy parameter
kappa3 = ; %Second anisotropy parameter
a = ; %Temperature parameter
Sn = ; %Linked to temperature parameter, see phase diagram
dt = 1; %SOR parameter
err = 1e-7; %Convergence parameter
numits = ; %Max number of iterations
%Load in lambdas
load('Lambda.mat');

%-------------------------------------------------------------------
[u,E] = tactoidSolnCubeMS(M-1,m,Rin,Rout,kappa2,kappa3,a,Sn,dt,err,numits,Lam1,Lam2,Lam3);
%filename
filename = ['tacSol_m',num2str(m),'_kappa2',num2str(kappa2),'.mat'];
save(['MSSols/',filename],'u','E');