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

rTSCMS1D.m
backEulerMS1D.m
getE1D.m
tridag.m
Lambda.mat
%}

M = ; %Number of grid spaces
psi0 = ; %Angle of director with respect to x-axis
Rout = ; %System size 
kappa2 = ; %First anisotropy parameter
a = 3.4049; %Temperature parameter
Sn = 0.4281; %Linked to temperature parameter, see phase diagram
err = 1e-7;
numits = ; %Max number of iterations
%Load in lambdas
load('Lambda.mat');

%-------------------------------------------------------------------
[eta,mu,nu,E] = backEulerMS1D(M,Rout,a,Sn,psi0,kappa2,err,numits,Lam1,Lam2,Lam3);
filename = ['interface_psi0',num2str(psi0),'_kappa2',num2str(kappa2),'.mat'];
%filename = strcat('MSSol_a=',num2str(a),'_k2=',num2str(kappa2),'_k3=',num2str(kappa3),'.mat');
save(['MSSols/',filename],'eta','mu','nu','E');