Numerically solve a (nonlinear) system of 64 differential equations

Hi all,

I'm trying to numerically solve a system of 64 equations with Matlab2014a and ode45 (see picture below for the equation, where Erf is the error function and Q is 64 by 64 matrix). But typing every single equation seems very tedious. So I was wondering (more like hoping) if there is some way to implement this formula in matlab without explicitly typing all the equations (and sum terms asf) i.e. using the notion below somehow with ode45?

Answer 1

[Walter Roberson]

N = 64;
X = sym(zeros(1, N));
lhs = sym(zeros(1,N));
rhs = sym(zeros(1,N));
for i = 1 : N
  X(i) = symfun( sprintf('x%d(t)', i), t );
  lhs(i) = diff(X(i), t);
end
for i = 1 : N
  rhs(i) = erf(sum( (Q(i,:) .* X) .* (1 + mu * X(i)) ) - lambda*X(i));
end

Now you can go for

dsolve(lhs == rhs)

However, dsolve() is not able to solve the system symbolically. dsolve() is also not going to be able to solve the system numerically unless you also create boundary conditions.

You could proceed to

syms y
odefun = matlabFunction( rhs(:), 'vars', {t y});

and then odefun would be a function handle returning a column vector of differential values.

However... each of your xi and xj is not a variable but is instead a function xi(t) xj(t), so the odefun that was generated in this way would contain calls to functions x1, x2, x3, .... x64 . And that is a problem because you do not know what the functions are.

Perhaps, then:

N = 64;
X = sym('x', [1, N]);
rhs = sym(zeros(1,N));
for i = 1 : N
  rhs(i) = erf(sum( (Q(i,:) .* X) .* (1 + mu * X(i)) ) - lambda*X(i));
end
odefun = matlabFunction( rhs(:), 'vars', {t X});

and then you could pass odefun as the first parameter to ode45

[holistic]

Thanks for the answer!

I tried to run the code:

N = 64;
X = sym('x', [1, N]);
rhs = sym(zeros(1,N));
for i = 1 : N
  rhs(i) = erf(sum( (Q(i,:) .* X) .* (1 + 0.05 * X(i)) ) - 0.15*X(i));
end
odefun = matlabFunction(rhs(:), 'vars', {t X});

but I'm getting the following error message:

Error using sym/matlabFunction>checkVars (line 157)
The second input must be a string, a cell of strings or sym arrays, or an array of symbolic
variables. 
Error in sym/matlabFunction (line 104)
vars = checkVars(funvars,opts);

Any idea, what the problem could be?

[Walter Roberson]

Which MATLAB version are you using?

[Walter Roberson]

For the above to work you would need

syms t

and you would have gotten a failure if you had a numeric t in scope.

This t corresponds to the fact that you asked for dxi/dt, implying that you want an ode function whose first parameter is named t and whose second parameter is x .

If you have a numeric t in scope, then you could use

odefun = matlabFunction(rhs(:), 'vars', {sym('t') X});

[holistic]

Thanks again!

I'm using Matlab2014a. Using symbolic t works now. But solving the system with ode45 gives me another error (sorry, no experience with the Matlab equation solving environment so far).

I ran the following code now with some given Matrix Q and some initial conditions given by init:

N = 64;
X = sym('x', [1, N]);
rhs = sym(zeros(1,N));
for i = 1 : N
  rhs(i) = erf(sum( (Q(i,:) .* X) .* (1 + 0.05 * X(i)) ) - 0.15*X(i));
end
odefun = matlabFunction(rhs(:), 'vars', {sym('t') X});

init=ones(1,64);

[T,Y] = ode45(odefun,[0 100],init);

The error I get after trying to run ode45 is:

Index exceeds matrix dimensions.

Error in
symengine>makeFhandle/@(t,in2)[-erf(in2(:,1).*(3.0./2.0e1)+in2(:,1).*(in2(:,1).*(1.0./2.0e1)+1.0).*7.99809683656279e-2+in2(:,2).*(in2(:,1).*(1.0./2.0e1)+1.0).*1.07414091919116e-1+in2(:,3).*(in2(:,1).*(1.0./2.0e1)+1.0).*3.314786330234999e-2+in2(:,4).*(in2(:,1).*(1.0./2.0e1)+1.0).*1.66198757047694e-17-in2(:,5).*(in2(:,1).*(1.0./2.0e1)+1.0).*1.580341064489109e-2-in2(:,6).*(in2(:,1).*(1.0./2.0e1)+1.0).*8.488241673305028e-2-in2(:,7).*(in2(:,1).*(1.0./2.0e1)+1.0).*2.045905882245121e-3-in2(:,8).*(in2(:,1).*(1.0./2.0e1)+1.0).*2.95904362181023e-17+in2(:,9).*(in2(:,1).*(1.0./2.0e1)+1.0).*1.419777972710162e-2+in2(:,18).*(in2(:,1).*(1.0./2.0e1)+1.0).*4.503567840706202e-3+in2(:,32).*(in2(:,1).*(1.0./2.0e1)+1.0).*1.764135287417524e-18-in2(:,36).*(in2(:,1).*(1.0./2.0e1)+1.0).*2.680116550681772e-3-in2(:,46).*(in2(:,1).*(1.0./2.0e1)+1.0).*1.907907552419566e-2-in2(:,47).*(in2(:,1).*(1.0./2.0e1)+1.0).*1.949607300943512e-2);-erf(in2(:,2).*(3.0./2.0e1)+in2(:,1).*(in2(:,2).*(1.0./2.0e1)+1.0).*3.488184486703885e-2+in2(:,2).*(in2(:,2).*(1.0./2.0e1)+1.0).*6.50657580042813e-2+in2(:,3).*(in2(:,2).*(1.0./2.0e1)+1.0).*1.265343658920942e-1+in2(:,4).*(in2(:,2).*(1.0./2.0e1)+1.0).*1.847419376107615e-17-in2(:,5).*(in2(:,2).*(1.0./2.0e1)+1.0).*8.943144198443781e-2-in2(:,6).*(in2(:,2).*(1.0./2.0e1)+1.0).*1.512080821669974e-1-in2(:,7).*(in2(:,2).*(1.0./2.0e1)+1.0).*9.397856507600125e-2-in2(:,8).*(in2(:,2).*(1.0./2.0e1)+1.0).*1.239247432784212e-17+in2(:,11).*(in2(:,2).*(1.0./2.0e1)+1.0).*1.925331758414416e-2+in2(:,15).*(in2(:,2).*(1.0./2.0e1)+1.0).*2.406361664217272e-2);-erf(in2(:,3).*(3.0./2.0e1)+in2(:,1).*(in2(:,3).*(1.0./2.0e1)+1.0).*8.81744922108718e-3+in2(:,2).*(in2(:,3).*(1.0./2.0e1)+1.0).*9.84746187318824e-2+in2(:,3).*(in2(:,3).*(1.0./2.0e1)+1.0).*3.741460468952067e-4+in2(:,4).*(in2(:,3).*(1.0./2.0e1)+1.0).*1.289789011370725e-17-in2(:,5).*(in2(:,3).*(1.0./2.0e1)+1.0).*1.378203831834125e-2-in2(:,6).*(in2(:,3).*(1.0./2.0e1)+1.0).*7.586531395896701e-2-in2(:,7).*(in2(:,3).*(1.0./2.0e1)+1.0).*9.401074831617101e-2-in2(:,8).*(in2(:,3).*(1.0./2.0e1)+1.0).*3.054444847054488e-17);-erf(in2(:,4).*(3.0./2.0e1)+in2(:,1).*(in2(:,4).*(1.0./2.0e1)+1.0).*8.909615715171684e-18+in2(:,2).*(in2(:,4).*(1.0./2.0e1)+1.0).*1.694202850300424e-18+in2(:,3).*(in2(:,4).*(1.0./2.0e1)+1.0).*1.159609179833099e-17+in2(:,4).*(in2(:,4).*(1.0./2.0e1)+1.0).*2.538586397098392e-33-in2(:,5).*(in2(:,4).*(1.0./2.0e1)+1.0).*7.949303445418413e-19-in2(:,6).*(in2(:,4).*(1.0./2.0e1)+1.0).*6.41826725793848e-18-in2(:,7).*(in2(:,4).*(1.0./2.0e1)+1.0).*1.073800535129517e-17-in2(:,8).*(in2(:,4).*(1.0./2.0e1)+1.0).*2.213387425827061e-33-in2(:,60).*(in2(:,4).*(1.0./2.0e1)+1.0).*1.551455650421417e-2);-erf(in2(:,5).*(3.0./2.0e1)-in2(:,1).*(in2(:,5).*(1.0./2.0e1)+1.0).*7.229263327801956e-2-in2(:,2).*(in2(:,5).*(1.0./2.0e1)+1.0).*6.879396588573246e-2-in2(:,3).*(in2(:,5).*(1.0./2.0e1)+1.0).*5.274693002223463e-2-in2(:,4).*(in2(:,5).*(1.0./2.0e1)+1.0).*1.416033938543907e-17+in2(:,5).*(in2(:,5).*(1.0./2.0e1)+1.0).*6.064430781703089e-2+in2(:,6).*(in2(:,5).*(1.0./2.0e1)+1.0).*8.224691476083075e-2+in2(:,7).*(in2(:,5).*(1.0./2.0e1)+1.0).*3.938410176440853e-2+in2(:,8).*(in2(:,5).*(1.0./2.0e1)+1.0).*3.19917421938282e-17-in2(:,46).*(in2(:,5).*(1.0./2.0e1)+1.0).*2.792798752717326e-2-in2(:,48).*(in2(:,5).*(1.0./2.0e1)+1.0).*1.151617405540475e-2);-erf(in2(:,6).*(3.0./2.0e1)-in2(:,1).*(in2(:,6).*(1.0./2.0e1)+1.0).*6.502896423675533e-2-in2(:,2).*(in2(:,6).*(1.0./2.0e1)+1.0).*8.6511404479342e-2-in2(:,3).*(in2(:,6).*(1.0./2.0e1)+1.0).*4.545559964353429e-2-in2(:,4).*(in2(:,6).*(1.0./2.0e1)+1.0).*1.89873977025125e-17+in2(:,5).*(in2(:,6).*(1.0./2.0e1)+1.0).*7.661426129612155e-2+in2(:,6).*(in2(:,6).*(1.0./2.0e1)+1.0).*6.665801620908114e-2+in2(:,7).*(in2(:,6).*(1.0./2.0e1)+1.0).*4.582526993230324e-2+in2(:,8).*(in2(:,6).*(1.0./2.0e1)+1.0).*1.162129508321451e-17+in2(:,14).*(in2(:,6).*(1.0./2.0e1)+1.0).*5.764268645577905e-2+in2(:,17).*(in2(:,6).*(1.0./2.0e1)+1.0).*2.481389795531033e-2-in2(:,36).*(in2(:,6).*(1.0./2.0e1)+1.0).*1.896801189431572e-2-in2(:,42).*(in2(:,6).*(1.0./2.0e1)+1.0).*7.420710624832825e-3);-erf(in2(:,7).*(3.0./2.0e1)-in2(:,1).*(in2(:,7).*(1.0./2.0e1)+1.0).*7.796265056908606e-3-in2(:,2).*(in2(:,7).*(1.0./2.0e1)+1.0).*1.17044487098766e-1-in2(:,3).*(in2(:,7).*(1.0./2.0e1)+1.0).*2.06831449925103e-3-in2(:,4).*(in2(:,7).*(1.0./2.0e1)+1.0).*1.34836975731288e-17+in2(:,5).*(in2(:,7).*(1.0./2.0e1)+1.0).*2.681747473972178e-2+in2(:,6).*(in2(:,7).*(1.0./2.0e1)+1.0).*1.8694165750749e-2+in2(:,7).*(in2(:,7).*(1.0./2.0e1)+1.0).*6.824711111573827e-3+in2(:,8).*(in2(:,7).*(1.0./2.0e1)+1.0).*1.999753374748547e-17-in2(:,38).*(in2(:,7).*(1.0./2.0e1)+1.0).*5.153560821885156e-2-in2(:,56).*(in2(:,7).*(1.0./2.0e1)+1.0).*3.688463656669122e-2-in2(:,62).*(in2(:,7).*(1.0./2.0e1)+1.0).*2.446566557116911e-2);-erf(in2(:,8).*(3.0./2.0e1)-in2(:,1).*(in2(:,8).*(1.0./2.0e1)+1.0).*4.797252442398348e-18-in2(:,2).*(in2(:,8).*(1.0./2.0e1)+1.0).*1.468406612758519e-17-in2(:,3).*(in2(:,8).*(1.0./2.0e1)+1.0).*7.964571877758049e-18-in2(:,4).*(in2(:,8).*(1.0./2.0e1)+1.0).*3.557220374340735e-33+in2(:,5).*(in2(:,8).*(1.0./2.0e1)+1.0).*2.752707904651532e-17+in2(:,6).*(in2(:,8).*(1.0./2.0e1)+1.0).*3.280068254728857e-17+in2(:,7).*(in2(:,8).*(1.0./2.0e1)+1.0).*2.554656599572182e-17+in2(:,8).*(in2(:,8).*(1.0./2.0e1)+1.0).*1.194027822280651e-32+in2(:,10).*(in2(:,8).*(1.0./2.0e1)+1.0).*9.677976619425696e-3+in2(:,29).*(in2(:,8).*(1.0./2.0e1)+1.0).*1.370423753194491e-2-in2(:,43).*(in2(:,8).*(1.0./2.0e1)+1.0).*5.499782528961678e-2-in2(:,55).*(in2(:,8).*(1.0./2.0e1)+1.0).*1.087637102369617e-1);-erf(in2(:,9).*(3.0./2.0e1)+in2(:,8).*(in2(:,9).*(1.0./2.0e1)+1.0).*2.909708459482493e-2+in2(:,9).*(in2(:,9).*(1.0./2.0e1)+1.0).*4.487343177616468e-3+in2(:,10).*(in2(:,9).*(1.0./2.0e1)+1.0).*5.768584528897285e-2+in2(:,11).*(in2(:,9).*(1.0./2.0e1)+1.0).*4.107306215416469e-2+in2(:,12).*(in2(:,9).*(1.0./2.0e1)+1.0).*7.356909107536882e-18-in2(:,13).*(in2(:,9).*(1.0./2.0e1)+1.0).*1.873751840825473e-2-in2(:,14).*(in2(:,9).*(1.0./2.0e1)+1.0).*1.425234359908031e-1-in2(:,15).*(in2(:,9).*(1.0./2.0e1)+1.0).*6.984027170493991e-2-in2(:,16).*(in2(:,9).*(1.0./2.0e1)+1.0).*3.1557422512336e-17+in2(:,20).*(in2(:,9).*(1.0./2.0e1)+1.0).*1.260683260403629e-1);-erf(in2(:,10).*(3.0./2.0e1)+in2(:,1).*(in2(:,10).*(1.0./2.0e1)+1.0).*1.126107857314477e-2+in2(:,3).*(in2(:,10).*(1.0./2.0e1)+1.0).*2.639767294478128e-2+in2(:,9).*(in2(:,10).*(1.0./2.0e1)+1.0).*2.430476481872912e-2+in2(:,10).*(in2(:,10).*(1.0./2.0e1)+1.0).*1.473669731314094e-3+in2(:,11).*(in2(:,10).*(1.0./2.0e1)+1.0).*3.711458155824415e-2+in2(:,12).*(in2(:,10).*(1.0./2.0e1)+1.0).*2.33368928318137e-17-in2(:,13).*(in2(:,10).*(1.0./2.0e1)+1.0).*7.649899469148876e-2-in2(:,14).*(in2(:,10).*(1.0./2.0e1)+1.0).*3.615764250330596e-2-in2(:,15).*(in2(:,10).*(1.0./2.0e1)+1.0).*1.134256111802032e-1-in2(:,16).*(in2(:,10).*(1.0./2.0e1)+1.0).*2.179630312407674e-17+in2(:,22).*(in2(:,10).*(1.0./2.0e1)+1.0).*1.467313539922843e-1+in2(:,30).*(in2(:,10).*(1.0./2.0e1)+1.0).*1.513649368077297e-2+in2(:,31).*(in2(:,10).*(1.0./2.0e1)+1.0).*1.288639950955502e-2);-erf(in2(:,11).*(3.0./2.0e1)+in2(:,6).*(in2(:,11).*(1.0./2.0e1)+1.0).*1.029739744720864e-2+in2(:,9).*(in2(:,11).*(1.0./2.0e1)+1.0).*3.168442596866235e-2+in2(:,10).*(in2(:,11).*(1.0./2.0e1)+1.0).*1.064925985089645e-2+in2(:,11).*(in2(:,11).*(1.0./2.0e1)+1.0).*8.508887696246523e-2+in2(:,12).*(in2(:,11).*(1.0./2.0e1)+1.0).*1.340906737649788e-17-in2(:,13).*(in2(:,11).*(1.0./2.0e1)+1.0).*7.075321819831534e-2-in2(:,14).*(in2(:,11).*(1.0./2.0e1)+1.0).*1.08321016952543e-1-in2(:,15).*(in2(:,11).*(1.0./2.0e1)+1.0).*8.240596351581503e-2-in2(:,16).*(in2(:,11).*(1.0./2.0e1)+1.0).*3.435436429814929e-17-in2(:,34).*(in2(:,11).*(1.0./2.0e1)+1.0).*3.323681068388303e-2-in2(:,35).*(in2(:,11).*(1.0./2.0e1)+1.0).*1.827107325922317e-2-in2(:,49).*(in2(:,11).*(1.0./2.0e1)+1.0).*1.72505951906326e-1);-erf(in2(:,12).*(3.0./2.0e1)+in2(:,6).*(in2(:,12).*(1.0./2.0e1)+1.0).*8.216806704978673e-2+in2(:,9).*(in2(:,12).*(1.0./2.0e1)+1.0).*1.682218659854142e-17+in2(:,10).*(in2(:,12).*(1.0./2.0e1)+1.0).*1.973372319359783e-17+in2(:,11).*(in2(:,12).*(1.0./2.0e1)+1.0).*3.903683990982813e-18+in2(:,12).*(in2(:,12).*(1.0./2.0e1)+1.0).*2.813539371412228e-33-in2(:,13).*(in2(:,12).*(1.0./2.0e1)+1.0).*1.747336017045046e-17-in2(:,14).*(in2(:,12).*(1.0./2.0e1)+1.0).*1.269147321280003e-17-in2(:,15).*(in2(:,12).*(1.0./2.0e1)+1.0).*1.623877928057712e-17-in2(:,16).*(in2(:,12).*(1.0./2.0e1)+1.0).*1.843613490785155e-33+in2(:,28).*(in2(:,12).*(1.0./2.0e1)+1.0).*2.512107623619051e-2-in2(:,55).*(in2(:,12).*(1.0./2.0e1)+1.0).*1.432296079988823e-2);-erf(in2(:,13).*(3.0./2.0e1)-in2(:,9).*(in2(:,13).*(1.0./2.0e1)+1.0).*6.73461108777634e-2-in2(:,10).*(in2(:,13).*(1.0./2.0e1)+1.0).*4.734788820410554e-2-in2(:,11).*(in2(:,13).*(1.0./2.0e1)+1.0).*1.347604970571522e-2-in2(:,12).*(in2(:,13).*(1.0./2.0e1)+1.0).*1.519749087590005e-17+in2(:,13).*(in2(:,13).*(1.0./2.0e1)+1.0).*2.682477206651962e-2+in2(:,14).*(in2(:,13).*(1.0./2.0e1)+1.0).*1.288122102415184e-1+in2(:,15).*(in2(:,13).*(1.0./2.0e1)+1.0).*9.200645408218706e-2+in2(:,16).*(in2(:,13).*(1.0./2.0e1)+1.0).*2.309872179996475e-17-in2(:,40).*(in2(:,13).*(1.0./2.0e1)+1.0).*7.80335508170457e-2-in2(:,62).*(in2(:,13).*(1.0./2.0e1)+1.0).*2.399825613240164e-2);-erf(in2(:,14).*(3.0./2.0e1)-in2(:,9).*(in2(:,14).*(1.0./2.0e1)+1.0).*1.114815723579773e-1-in2(:,10).*(in2(:,14).*(1.0./2.0e1)+1.0).*1.641097822423286e-1-in2(:,11).*(in2(:,14).*(1.0./2.0e1)+1.0).*9.881830631748262e-2-in2(:,12).*(in2(:,14).*(1.0./2.0e1)+1.0).*1.950390532377712e-17+in2(:,13).*(in2(:,14).*(1.0./2.0e1)+1.0).*6.551536405350196e-3+in2(:,14).*(in2(:,14).*(1.0./2.0e1)+1.0).*1.098221135558923e-2+in2(:,15).*(in2(:,14).*(1.0./2.0e1)+1.0).*1.351232694994529e-1+in2(:,16).*(in2(:,14).*(1.0./2.0e1)+1.0).*3.294434229211492e-18-in2(:,39).*(in2(:,14).*(1.0./2.0e1)+1.0).*7.423049826156938e-2-in2(:,46).*(in2(:,14).*(1.0./2.0e1)+1.0).*1.069381957651661e-1-in2(:,63).*(in2(:,14).*(1.0./2.0e1)+1.0).*1.447133883067861e-2);-erf(in2(:,15).*(3.0./2.0e1)+in2(:,8).*(in2(:,15).*(1.0./2.0e1)+1.0).*1.28942670881728e-1-in2(:,9).*(in2(:,15).*(1.0./2.0e1)+1.0).*5.533456222808266e-2-in2(:,10).*(in2(:,15).*(1.0./2.0e1)+1.0).*1.34185179504109e-1-in2(:,11).*(in2(:,15).*(1.0./2.0e1)+1.0).*6.893849240065386e-2-in2(:,12).*(in2(:,15).*(1.0./2.0e1)+1.0).*3.414658262112984e-18+in2(:,13).*(in2(:,15).*(1.0./2.0e1)+1.0).*2.747298852612306e-2+in2(:,14).*(in2(:,15).*(1.0./2.0e1)+1.0).*7.452392629843343e-2+in2(:,15).*(in2(:,15).*(1.0./2.0e1)+1.0).*6.003985559181255e-2+in2(:,16).*(in2(:,15).*(1.0./2.0e1)+1.0).*1.438397266444711e-17+in2(:,27).*(in2(:,15).*(1.0./2.0e1)+1.0).*5.887893547322604e-2-in2(:,41).*(in2(:,15).*(1.0./2.0e1)+1.0).*8.687352488839267e-2);-erf(in2(:,16).*(3.0./2.0e1)-in2(:,9).*(in2(:,16).*(1.0./2.0e1)+1.0).*2.346918360226473e-17-in2(:,10).*(in2(:,16).*(1.0./2.0e1)+1.0).*3.436091959418059e-17-in2(:,11).*(in2(:,16).*(1.0./2.0e1)+1.0).*1.179946483769208e-17-in2(:,12).*(in2(:,16).*(1.0./2.0e1)+1.0).*2.240664777687601e-33+in2(:,13).*(in2(:,16).*(1.0./2.0e1)+1.0).*1.398992350037397e-17+in2(:,14).*(in2(:,16).*(1.0./2.0e1)+1.0).*4.143647845299853e-17+in2(:,15).*(in2(:,16).*(1.0./2.0e1)+1.0).*9.032040778515777e-18+in2(:,16).*(in2(:,16).*(1.0./2.0e1)+1.0).*1.202641919627294e-32-in2(:,42).*(in2(:,16).*(1.0./2.0e1)+1.0).*4.600036221229821e-2-in2(:,50).*(in2(:,16).*(1.0./2.0e1)+1.0).*6.728038052866271e-2);-erf(in2(:,17).*(3.0./2.0e1)+in2(:,17).*(in2(:,17).*(1.0./2.0e1)+1.0).*1.783762565943796e-2+in2(:,18).*(in2(:,17).*(1.0./2.0e1)+1.0).*1.385065582385652e-2+in2(:,19).*(in2(:,17).*(1.0./2.0e1)+1.0).*4.739111882253826e-2+in2(:,20).*(in2(:,17).*(1.0./2.0e1)+1.0).*4.891642143552679e-18-in2(:,21).*(in2(:,17).*(1.0./2.0e1)+1.0).*1.001997866192885e-1-in2(:,22).*(in2(:,17).*(1.0./2.0e1)+1.0).*7.919908573630682e-2-in2(:,23).*(in2(:,17).*(1.0./2.0e1)+1.0).*4.188062755500002e-2-in2(:,24).*(in2(:,17).*(1.0./2.0e1)+1.0).*2.398284945360468e-17);-erf(in2(:,18).*(3.0./2.0e1)+in2(:,9).*(in2(:,18).*(1.0./2.0e1)+1.0).*6.611750161764809e-3+in2(:,14).*(in2(:,18).*(1.0./2.0e1)+1.0).*1.921664200426792e-1+in2(:,15).*(in2(:,18).*(1.0./2.0e1)+1.0).*3.696438738684418e-2+in2(:,17).*(in2(:,18).*(1.0./2.0e1)+1.0).*1.367337434597279e-1+in2(:,18).*(in2(:,18).*(1.0./2.0e1)+1.0).*5.404209012637051e-2+in2(:,19).*(in2(:,18).*(1.0./2.0e1)+1.0).*1.019942377678672e-1+in2(:,20).*(in2(:,18).*(1.0./2.0e1)+1.0).*2.452525444925389e-17-in2(:,21).*(in2(:,18).*(1.0./2.0e1)+1.0).*4.678596399559092e-2-in2(:,22).*(in2(:,18).*(1.0./2.0e1)+1.0).*1.959616886865309e-1-in2(:,23).*(in2(:,18).*(1.0./2.0e1)+1.0).*3.538928514851802e-2-in2(:,24).*(in2(:,18).*(1.0./2.0e1)+1.0).*9.197309215958504e-18-in2(:,37).*(in2(:,18).*(1.0./2.0e1)+1.0).*3.037291777057086e-2-in2(:,46).*(in2(:,18).*(1.0./2.0e1)+1.0).*8.635510145821652e-2-in2(:,56).*(in2(:,18).*(1.0./2.0e1)+1.0).*9.369980737661128e-2);-erf(in2(:,19).*(3.0./2.0e1)+in2(:,17).*(in2(:,19).*(1.0./2.0e1)+1.0).*6.462901027074029e-2+in2(:,18).*(in2(:,19).*(1.0./2.0e1)+1.0).*2.915372312377837e-2+in2(:,19).*(in2(:,19).*(1.0./2.0e1)+1.0).*8.205689517351775e-2+in2(:,20).*(in2(:,19).*(1.0./2.0e1)+1.0).*1.118076177097606e-17-in2(:,21).*(in2(:,19).*(1.0./2.0e1)+1.0).*7.582461017218178e-2-in2(:,22).*(in2(:,19).*(1.0./2.0e1)+1.0).*8.667241300130757e-2-in2(:,23).*(in2(:,19).*(1.0./2.0e1)+1.0).*9.261253202761083e-2-in2(:,24).*(in2(:,19).*(1.0./2.0e1)+1.0).*1.782746002279313e-17);-erf(in2(:,20).*(3.0./2.0e1)+in2(:,2).*(in2(:,20).*(1.0./2.0e1)+1.0).*1.213051458070598e-4+in2(:,17).*(in2(:,20).*(1.0./2.0e1)+1.0).*6.750265891388768e-18+in2(:,18).*(in2(:,20).*(1.0./2.0e1)+1.0).*2.357823312662391e-18+in2(:,19).*(in2(:,20).*(1.0./2.0e1)+1.0).*5.713606156453875e-18+in2(:,20).*(in2(:,20).*(1.0./2.0e1)+1.0).*2.135298182243787e-33-in2(:,21).*(in2(:,20).*(1.0./2.0e1)+1.0).*1.46479875080115e-17-in2(:,22).*(in2(:,20).*(1.0./2.0e1)+1.0).*2.384758751628681e-17-in2(:,23).*(in2(:,20).*(1.0./2.0e1)+1.0).*7.950091337210929e-18-in2(:,24).*(in2(:,20).*(1.0./2.0e1)+1.0).*7.598138029656169e-34);-erf(in2(:,21).*(3.0./2.0e1)+in2(:,12).*(in2(:,21).*(1.0./2.0e1)+1.0).*6.156637501118412e-2-in2(:,17).*(in2(:,21).*(1.0./2.0e1)+1.0).*2.647511501864073e-2-in2(:,18).*(in2(:,21).*(1.0./2.0e1)+1.0).*6.193195070583972e-2-in2(:,19).*(in2(:,21).*(1.0./2.0e1)+1.0).*6.276375153405717e-2-in2(:,20).*(in2(:,21).*(1.0./2.0e1)+1.0).*1.58392443711654e-17+in2(:,21).*(in2(:,21).*(1.0./2.0e1)+1.0).*2.96525254930462e-2+in2(:,22).*(in2(:,21).*(1.0./2.0e1)+1.0).*9.981967842577315e-3+in2(:,23).*(in2(:,21).*(1.0./2.0e1)+1.0).*3.164966181419147e-2+in2(:,24).*(in2(:,21).*(1.0./2.0e1)+1.0).*1.417443304190425e-17-in2(:,53).*(in2(:,21).*(1.0./2.0e1)+1.0).*1.513932637398565e-1-in2(:,57).*(in2(:,21).*(1.0./2.0e1)+1.0).*1.170715571094548e-1-in2(:,59).*(in2(:,21).*(1.0./2.0e1)+1.0).*8.557462491790178e-3);-erf(in2(:,22).*(3.0./2.0e1)+in2(:,11).*(in2(:,22).*(1.0./2.0e1)+1.0).*4.946264380676689e-2-in2(:,17).*(in2(:,22).*(1.0./2.0e1)+1.0).*8.220547494822616e-2-in2(:,18).*(in2(:,22).*(1.0./2.0e1)+1.0).*1.21675255268487e-1-in2(:,19).*(in2(:,22).*(1.0./2.0e1)+1.0).*1.011713027056616e-2-in2(:,20).*(in2(:,22).*(1.0./2.0e1)+1.0).*3.471624086032824e-18+in2(:,21).*(in2(:,22).*(1.0./2.0e1)+1.0).*7.163736554997781e-2+in2(:,22).*(in2(:,22).*(1.0./2.0e1)+1.0).*1.508049606340039e-1+in2(:,23).*(in2(:,22).*(1.0./2.0e1)+1.0).*1.142403143334665e-1+in2(:,24).*(in2(:,22).*(1.0./2.0e1)+1.0).*2.835095952755499e-18);-erf(in2(:,23).*(3.0./2.0e1)+in2(:,14).*(in2(:,23).*(1.0./2.0e1)+1.0).*7.57574050868707e-2-in2(:,17).*(in2(:,23).*(1.0./2.0e1)+1.0).*6.598650843886798e-2-in2(:,18).*(in2(:,23).*(1.0./2.0e1)+1.0).*6.165632537104401e-3-in2(:,19).*(in2(:,23).*(1.0./2.0e1)+1.0).*1.856402501911107e-2-in2(:,20).*(in2(:,23).*(1.0./2.0e1)+1.0).*3.035197088397327e-18+in2(:,21).*(in2(:,23).*(1.0./2.0e1)+1.0).*5.820352719025716e-2+in2(:,22).*(in2(:,23).*(1.0./2.0e1)+1.0).*8.348613537562961e-2+in2(:,23).*(in2(:,23).*(1.0./2.0e1)+1.0).*6.123417766597056e-2+in2(:,24).*(in2(:,23).*(1.0./2.0e1)+1.0).*8.190253713099672e-18+in2(:,29).*(in2(:,23).*(1.0./2.0e1)+1.0).*1.323659653569698e-2-in2(:,61).*(in2(:,23).*(1.0./2.0e1)+1.0).*1.177740525961636e-2);-erf(in2(:,24).*(3.0./2.0e1)+in2(:,4).*(in2(:,24).*(1.0./2.0e1)+1.0).*4.297520093347213e-2+in2(:,10).*(in2(:,24).*(1.0./2.0e1)+1.0).*4.172768164530004e-2-in2(:,17).*(in2(:,24).*(1.0./2.0e1)+1.0).*3.431894960466483e-17-in2(:,18).*(in2(:,24).*(1.0./2.0e1)+1.0).*5.806380941697989e-19-in2(:,19).*(in2(:,24).*(1.0./2.0e1)+1.0).*2.230465942260514e-17-in2(:,20).*(in2(:,24).*(1.0./2.0e1)+1.0).*1.953812920297034e-33+in2(:,21).*(in2(:,24).*(1.0./2.0e1)+1.0).*3.472938305726699e-17+in2(:,22).*(in2(:,24).*(1.0./2.0e1)+1.0).*4.626764411020819e-17+in2(:,23).*(in2(:,24).*(1.0./2.0e1)+1.0).*1.046054652674069e-17+in2(:,24).*(in2(:,24).*(1.0./2.0e1)+1.0).*5.670935702949456e-33);-erf(in2(:,25).*(3.0./2.0e1)+in2(:,25).*(in2(:,25).*(1.0./2.0e1)+1.0).*5.802700030511466e-2+in2(:,26).*(in2(:,25).*(1.0./2.0e1)+1.0).*1.32800381117189e-2+in2(:,27).*(in2(:,25).*(1.0./2.0e1)+1.0).*1.326862076283787e-2+in2(:,28).*(in2(:,25).*(1.0./2.0e1)+1.0).*5.832857290119141e-18-in2(:,29).*(in2(:,25).*(1.0./2.0e1)+1.0).*7.716284073572363e-2-in2(:,30).*(in2(:,25).*(1.0./2.0e1)+1.0).*9.096810153061753e-2-in2(:,31).*(in2(:,25).*(1.0./2.0e1)+1.0).*3.342946960359679e-2-in2(:,32).*(in2(:,25).*(1.0./2.0e1)+1.0).*2.7272582474251e-17-in2(:,38).*(in2(:,25).*(1.0./2.0e1)+1.0).*1.64747328317547e-2-in2(:,54).*(in2(:,25).*(1.0./2.0e1)+1.0).*4.577683307318399e-2-in2(:,55).*(in2(:,25).*(1.0./2.0e1)+1.0).*2.202716527665889e-2-in2(:,56).*(in2(:,25).*(1.0./2.0e1)+1.0).*4.04255226272249e-2);-erf(in2(:,26).*(3.0./2.0e1)+in2(:,3).*(in2(:,26).*(1.0./2.0e1)+1.0).*8.898885700813904e-4+in2(:,7).*(in2(:,26).*(1.0./2.0e1)+1.0).*6.929510118261703e-3+in2(:,8).*(in2(:,26).*(1.0./2.0e1)+1.0).*3.094287962363245e-2+in2(:,18).*(in2(:,26).*(1.0./2.0e1)+1.0).*8.533302848329484e-2+in2(:,25).*(in2(:,26).*(1.0./2.0e1)+1.0).*8.755265278141711e-2+in2(:,26).*(in2(:,26).*(1.0./2.0e1)+1.0).*1.761050268924475e-1+in2(:,27).*(in2(:,26).*(1.0./2.0e1)+1.0).*7.408889797706669e-2+in2(:,28).*(in2(:,26).*(1.0./2.0e1)+1.0).*1.900357195573185e-17-in2(:,29).*(in2(:,26).*(1.0./2.0e1)+1.0).*7.176660838582206e-2-in2(:,30).*(in2(:,26).*(1.0./2.0e1)+1.0).*1.142033024880102e-1-in2(:,31).*(in2(:,26).*(1.0./2.0e1)+1.0).*9.712681819543392e-2-in2(:,32).*(in2(:,26).*(1.0./2.0e1)+1.0).*6.151818679764457e-18-in2(:,36).*(in2(:,26).*(1.0./2.0e1)+1.0).*2.619734438674718e-2-in2(:,39).*(in2(:,26).*(1.0./2.0e1)+1.0).*6.076992327036485e-2-in2(:,54).*(in2(:,26).*(1.0./2.0e1)+1.0).*1.894295872650973e-3-in2(:,55).*(in2(:,26).*(1.0./2.0e1)+1.0).*3.284481207759014e-2-in2(:,59).*(in2(:,26).*(1.0./2.0e1)+1.0).*3.210454872182061e-2);-erf(in2(:,27).*(3.0./2.0e1)+in2(:,21).*(in2(:,27).*(1.0./2.0e1)+1.0).*4.73878974599412e-2+in2(:,25).*(in2(:,27).*(1.0./2.0e1)+1.0).*3.437773563873966e-2+in2(:,26).*(in2(:,27).*(1.0./2.0e1)+1.0).*1.240619377388385e-1+in2(:,27).*(in2(:,27).*(1.0./2.0e1)+1.0).*9.315932514020929e-2+in2(:,28).*(in2(:,27).*(1.0./2.0e1)+1.0).*6.621633178641045e-18-in2(:,29).*(in2(:,27).*(1.0./2.0e1)+1.0).*7.47269159000899e-2-in2(:,30).*(in2(:,27).*(1.0./2.0e1)+1.0).*1.423683278989094e-1-in2(:,31).*(in2(:,27).*(1.0./2.0e1)+1.0).*8.873576493424778e-2-in2(:,32).*(in2(:,27).*(1.0./2.0e1)+1.0).*5.792553302413692e-18-in2(:,34).*(in2(:,27).*(1.0./2.0e1)+1.0).*1.769297078472749e-2);-erf(in2(:,28).*(3.0./2.0e1)+in2(:,12).*(in2(:,28).*(1.0./2.0e1)+1.0).*5.479020864175959e-2+in2(:,21).*(in2(:,28).*(1.0./2.0e1)+1.0).*3.226396558192091e-2+in2(:,25).*(in2(:,28).*(1.0./2.0e1)+1.0).*4.965487647711074e-18+in2(:,26).*(in2(:,28).*(1.0./2.0e1)+1.0).*1.437258320109412e-17+in2(:,27).*(in2(:,28).*(1.0./2.0e1)+1.0).*1.517753508164512e-17+in2(:,28).*(in2(:,28).*(1.0./2.0e1)+1.0).*5.468331803678522e-34-in2(:,29).*(in2(:,28).*(1.0./2.0e1)+1.0).*1.220066077524435e-17-in2(:,30).*(in2(:,28).*(1.0./2.0e1)+1.0).*2.327114255494634e-17-in2(:,31).*(in2(:,28).*(1.0./2.0e1)+1.0).*2.405606200665338e-18-in2(:,32).*(in2(:,28).*(1.0./2.0e1)+1.0).*4.249760435411767e-33-in2(:,33).*(in2(:,28).*(1.0./2.0e1)+1.0).*4.395957814050327e-3-in2(:,42).*(in2(:,28).*(1.0./2.0e1)+1.0).*5.531327808355227e-2-in2(:,44).*(in2(:,28).*(1.0./2.0e1)+1.0).*2.672567914524327e-2);-erf(in2(:,29).*(3.0./2.0e1)+in2(:,16).*(in2(:,29).*(1.0./2.0e1)+1.0).*8.434381234764756e-3-in2(:,25).*(in2(:,29).*(1.0./2.0e1)+1.0).*7.167751938106152e-2-in2(:,26).*(in2(:,29).*(1.0./2.0e1)+1.0).*5.422424967384272e-2-in2(:,27).*(in2(:,29).*(1.0./2.0e1)+1.0).*6.000361466599812e-2-in2(:,28).*(in2(:,29).*(1.0./2.0e1)+1.0).*1.177057652709676e-17+in2(:,29).*(in2(:,29).*(1.0./2.0e1)+1.0).*1.137399062007285e-2+in2(:,30).*(in2(:,29).*(1.0./2.0e1)+1.0).*5.856357742170167e-2+in2(:,31).*(in2(:,29).*(1.0./2.0e1)+1.0).*1.410762198848606e-2+in2(:,32).*(in2(:,29).*(1.0./2.0e1)+1.0).*1.400266207545345e-17);-erf(in2(:,30).*(3.0./2.0e1)+in2(:,19).*(in2(:,30).*(1.0./2.0e1)+1.0).*8.325398482915281e-4-in2(:,25).*(in2(:,30).*(1.0./2.0e1)+1.0).*2.693782549688826e-2-in2(:,26).*(in2(:,30).*(1.0./2.0e1)+1.0).*1.148251921481302e-1-in2(:,27).*(in2(:,30).*(1.0./2.0e1)+1.0).*1.408549950359755e-1-in2(:,28).*(in2(:,30).*(1.0./2.0e1)+1.0).*1.690259405850149e-17+in2(:,29).*(in2(:,30).*(1.0./2.0e1)+1.0).*6.16544496801162e-2+in2(:,30).*(in2(:,30).*(1.0./2.0e1)+1.0).*6.202427797811995e-3+in2(:,31).*(in2(:,30).*(1.0./2.0e1)+1.0).*3.064192163439966e-2+in2(:,32).*(in2(:,30).*(1.0./2.0e1)+1.0).*4.423104961327299e-17-in2(:,53).*(in2(:,30).*(1.0./2.0e1)+1.0).*1.374806025867751e-2);-erf(in2(:,31).*(3.0./2.0e1)+in2(:,3).*(in2(:,31).*(1.0./2.0e1)+1.0).*4.350995006986503e-3-in2(:,25).*(in2(:,31).*(1.0./2.0e1)+1.0).*6.087315035747156e-2-in2(:,26).*(in2(:,31).*(1.0./2.0e1)+1.0).*9.894671323570833e-2-in2(:,27).*(in2(:,31).*(1.0./2.0e1)+1.0).*2.172261651082192e-2-in2(:,28).*(in2(:,31).*(1.0./2.0e1)+1.0).*1.054137438234302e-17+in2(:,29).*(in2(:,31).*(1.0./2.0e1)+1.0).*9.291047409005489e-2+in2(:,30).*(in2(:,31).*(1.0./2.0e1)+1.0).*7.376251266410884e-3+in2(:,31).*(in2(:,31).*(1.0./2.0e1)+1.0).*4.951789018796796e-2+in2(:,32).*(in2(:,31).*(1.0./2.0e1)+1.0).*1.61402566788816e-17);-erf(in2(:,32).*(3.0./2.0e1)-in2(:,25).*(in2(:,32).*(1.0./2.0e1)+1.0).*1.455539512054984e-17-in2(:,26).*(in2(:,32).*(1.0./2.0e1)+1.0).*4.619747991108028e-17-in2(:,27).*(in2(:,32).*(1.0./2.0e1)+1.0).*3.305407350605798e-17-in2(:,28).*(in2(:,32).*(1.0./2.0e1)+1.0).*5.859516181454742e-33+in2(:,29).*(in2(:,32).*(1.0./2.0e1)+1.0).*8.726323549448075e-18+in2(:,30).*(in2(:,32).*(1.0./2.0e1)+1.0).*4.395458488396153e-17+in2(:,31).*(in2(:,32).*(1.0./2.0e1)+1.0).*3.435588264990492e-17+in2(:,32).*(in2(:,32).*(1.0./2.0e1)+1.0).*1.889151301369635e-34-in2(:,44).*(in2(:,32).*(1.0./2.0e1)+1.0).*1.228273284993815e-17-in2(:,53).*(in2(:,32).*(1.0./2.0e1)+1.0).*6.320645900421326e-18);-erf(in2(:,33).*(3.0./2.0e1)-in2(:,11).*(in2(:,33).*(1.0./2.0e1)+1.0).*6.477408857562655e-3+in2(:,33).*(in2(:,33).*(1.0./2.0e1)+1.0).*4.408002555500257e-2+in2(:,34).*(in2(:,33).*(1.0./2.0e1)+1.0).*1.165025652418155e-1+in2(:,35).*(in2(:,33).*(1.0./2.0e1)+1.0).*4.983700488288356e-2+in2(:,36).*(in2(:,33).*(1.0./2.0e1)+1.0).*1.338155307749298e-17-in2(:,37).*(in2(:,33).*(1.0./2.0e1)+1.0).*6.295850809281162e-2-in2(:,38).*(in2(:,33).*(1.0./2.0e1)+1.0).*5.256036193045156e-2-in2(:,39).*(in2(:,33).*(1.0./2.0e1)+1.0).*4.473232274010527e-2-in2(:,40).*(in2(:,33).*(1.0./2.0e1)+1.0).*1.836627221336394e-17);-erf(in2(:,34).*(3.0./2.0e1)+in2(:,33).*(in2(:,34).*(1.0./2.0e1)+1.0).*3.383304979320306e-2+in2(:,34).*(in2(:,34).*(1.0./2.0e1)+1.0).*4.308306222091548e-2+in2(:,35).*(in2(:,34).*(1.0./2.0e1)+1.0).*9.281782121238509e-2+in2(:,36).*(in2(:,34).*(1.0./2.0e1)+1.0).*7.113789798071499e-18-in2(:,37).*(in2(:,34).*(1.0./2.0e1)+1.0).*4.057306932978413e-2-in2(:,38).*(in2(:,34).*(1.0./2.0e1)+1.0).*5.593020392859453e-2-in2(:,39).*(in2(:,34).*(1.0./2.0e1)+1.0).*7.842018140233455e-2-in2(:,40).*(in2(:,34).*(1.0./2.0e1)+1.0).*1.329205961461043e-17);-erf(in2(:,35).*(3.0./2.0e1)-in2(:,14).*(in2(:,35).*(1.0./2.0e1)+1.0).*1.55515499252252e-4+in2(:,33).*(in2(:,35).*(1.0./2.0e1)+1.0).*2.901576425055537e-2+in2(:,34).*(in2(:,35).*(1.0./2.0e1)+1.0).*1.241355666265433e-1+in2(:,35).*(in2(:,35).*(1.0./2.0e1)+1.0).*8.984461795626322e-2+in2(:,36).*(in2(:,35).*(1.0./2.0e1)+1.0).*5.628888456595957e-18-in2(:,37).*(in2(:,35).*(1.0./2.0e1)+1.0).*6.647917711328435e-2-in2(:,38).*(in2(:,35).*(1.0./2.0e1)+1.0).*9.716206393913704e-2-in2(:,39).*(in2(:,35).*(1.0./2.0e1)+1.0).*3.931273363542892e-2-in2(:,40).*(in2(:,35).*(1.0./2.0e1)+1.0).*2.559724208428119e-17);-erf(in2(:,36).*(3.0./2.0e1)-in2(:,30).*(in2(:,36).*(1.0./2.0e1)+1.0).*9.240545913939352e-3+in2(:,33).*(in2(:,36).*(1.0./2.0e1)+1.0).*1.616407422422755e-17+in2(:,34).*(in2(:,36).*(1.0./2.0e1)+1.0).*8.10742208017701e-18+in2(:,35).*(in2(:,36).*(1.0./2.0e1)+1.0).*9.383001161652574e-18+in2(:,36).*(in2(:,36).*(1.0./2.0e1)+1.0).*2.602186889907447e-33-in2(:,37).*(in2(:,36).*(1.0./2.0e1)+1.0).*1.305193707305403e-17-in2(:,38).*(in2(:,36).*(1.0./2.0e1)+1.0).*1.021725853861858e-17-in2(:,39).*(in2(:,36).*(1.0./2.0e1)+1.0).*1.625573128550139e-17-in2(:,40).*(in2(:,36).*(1.0./2.0e1)+1.0).*3.645279128672181e-33+in2(:,48).*(in2(:,36).*(1.0./2.0e1)+1.0).*2.600768713713464e-2+in2(:,52).*(in2(:,36).*(1.0./2.0e1)+1.0).*5.052422697904606e-2);-erf(in2(:,37).*(3.0./2.0e1)-in2(:,33).*(in2(:,37).*(1.0./2.0e1)+1.0).*4.327028629987135e-2-in2(:,34).*(in2(:,37).*(1.0./2.0e1)+1.0).*2.739318274052711e-2-in2(:,35).*(in2(:,37).*(1.0./2.0e1)+1.0).*7.918444071329698e-3-in2(:,36).*(in2(:,37).*(1.0./2.0e1)+1.0).*5.763113816588418e-19+in2(:,37).*(in2(:,37).*(1.0./2.0e1)+1.0).*3.583228093791384e-2+in2(:,38).*(in2(:,37).*(1.0./2.0e1)+1.0).*8.758738823196419e-2+in2(:,39).*(in2(:,37).*(1.0./2.0e1)+1.0).*6.15738351717679e-2+in2(:,40).*(in2(:,37).*(1.0./2.0e1)+1.0).*3.344061198015423e-17);-erf(in2(:,38).*(3.0./2.0e1)-in2(:,3).*(in2(:,38).*(1.0./2.0e1)+1.0).*1.05485215633919e-2-in2(:,13).*(in2(:,38).*(1.0./2.0e1)+1.0).*3.321897419824225e-2-in2(:,33).*(in2(:,38).*(1.0./2.0e1)+1.0).*1.232620103622742e-1-in2(:,34).*(in2(:,38).*(1.0./2.0e1)+1.0).*1.969216952752632e-1-in2(:,35).*(in2(:,38).*(1.0./2.0e1)+1.0).*8.553839144048443e-2-in2(:,36).*(in2(:,38).*(1.0./2.0e1)+1.0).*2.46556860107291e-17+in2(:,37).*(in2(:,38).*(1.0./2.0e1)+1.0).*1.029236854295966e-1+in2(:,38).*(in2(:,38).*(1.0./2.0e1)+1.0).*1.324970144196391e-1+in2(:,39).*(in2(:,38).*(1.0./2.0e1)+1.0).*8.377776622533427e-2+in2(:,40).*(in2(:,38).*(1.0./2.0e1)+1.0).*1.80493149713789e-19+in2(:,58).*(in2(:,38).*(1.0./2.0e1)+1.0).*5.103057701674669e-2+in2(:,63).*(in2(:,38).*(1.0./2.0e1)+1.0).*6.607775094370032e-3);-erf(in2(:,39).*(3.0./2.0e1)-in2(:,5).*(in2(:,39).*(1.0./2.0e1)+1.0).*3.426219113109879e-2-in2(:,9).*(in2(:,39).*(1.0./2.0e1)+1.0).*1.15884293305716e-2-in2(:,14).*(in2(:,39).*(1.0./2.0e1)+1.0).*9.662151226856733e-2-in2(:,18).*(in2(:,39).*(1.0./2.0e1)+1.0).*7.708990909438783e-2-in2(:,32).*(in2(:,39).*(1.0./2.0e1)+1.0).*9.732223407491659e-18-in2(:,33).*(in2(:,39).*(1.0./2.0e1)+1.0).*8.568294867673731e-2-in2(:,34).*(in2(:,39).*(1.0./2.0e1)+1.0).*9.300597870078223e-2-in2(:,35).*(in2(:,39).*(1.0./2.0e1)+1.0).*8.731990880814858e-2-in2(:,36).*(in2(:,39).*(1.0./2.0e1)+1.0).*1.373553613210478e-18+in2(:,37).*(in2(:,39).*(1.0./2.0e1)+1.0).*4.462661693164974e-2+in2(:,38).*(in2(:,39).*(1.0./2.0e1)+1.0).*1.668418207416169e-2+in2(:,39).*(in2(:,39).*(1.0./2.0e1)+1.0).*4.737647558449145e-2+in2(:,40).*(in2(:,39).*(1.0./2.0e1)+1.0).*5.109987948310392e-18+in2(:,42).*(in2(:,39).*(1.0./2.0e1)+1.0).*4.018474151058674e-2+in2(:,57).*(in2(:,39).*(1.0./2.0e1)+1.0).*3.766501551733503e-2+in2(:,61).*(in2(:,39).*(1.0./2.0e1)+1.0).*2.55072673728904e-2);-erf(in2(:,40).*(3.0./2.0e1)-in2(:,10).*(in2(:,40).*(1.0./2.0e1)+1.0).*5.154261546049529e-3-in2(:,17).*(in2(:,40).*(1.0./2.0e1)+1.0).*1.283730367743965e-1-in2(:,29).*(in2(:,40).*(1.0./2.0e1)+1.0).*2.849068159993341e-2-in2(:,33).*(in2(:,40).*(1.0./2.0e1)+1.0).*1.581977085947768e-17-in2(:,34).*(in2(:,40).*(1.0./2.0e1)+1.0).*1.140068942305677e-17-in2(:,35).*(in2(:,40).*(1.0./2.0e1)+1.0).*2.599925074662511e-17-in2(:,36).*(in2(:,40).*(1.0./2.0e1)+1.0).*6.020910209557864e-33+in2(:,37).*(in2(:,40).*(1.0./2.0e1)+1.0).*2.15804302695048e-17+in2(:,38).*(in2(:,40).*(1.0./2.0e1)+1.0).*4.678927958628869e-17+in2(:,39).*(in2(:,40).*(1.0./2.0e1)+1.0).*2.954422065709902e-17+in2(:,40).*(in2(:,40).*(1.0./2.0e1)+1.0).*1.481856393803013e-33+in2(:,52).*(in2(:,40).*(1.0./2.0e1)+1.0).*9.011513959420857e-2);-erf(in2(:,41).*(3.0./2.0e1)-in2(:,4).*(in2(:,41).*(1.0./2.0e1)+1.0).*4.736724581859585e-2+in2(:,41).*(in2(:,41).*(1.0./2.0e1)+1.0).*6.32393599697848e-2+in2(:,42).*(in2(:,41).*(1.0./2.0e1)+1.0).*3.887387152402155e-2+in2(:,43).*(in2(:,41).*(1.0./2.0e1)+1.0).*3.508841415454338e-2+in2(:,44).*(in2(:,41).*(1.0./2.0e1)+1.0).*3.977281820326838e-18-in2(:,45).*(in2(:,41).*(1.0./2.0e1)+1.0).*8.476092671136909e-2-in2(:,46).*(in2(:,41).*(1.0./2.0e1)+1.0).*1.478699024727304e-3-in2(:,47).*(in2(:,41).*(1.0./2.0e1)+1.0).*6.918005567231084e-2-in2(:,48).*(in2(:,41).*(1.0./2.0e1)+1.0).*2.561815173741401e-18+in2(:,58).*(in2(:,41).*(1.0./2.0e1)+1.0).*7.439244695311474e-4);-erf(in2(:,42).*(3.0./2.0e1)-in2(:,20).*(in2(:,42).*(1.0./2.0e1)+1.0).*3.230276915399881e-3-in2(:,32).*(in2(:,42).*(1.0./2.0e1)+1.0).*1.479425086901272e-17+in2(:,41).*(in2(:,42).*(1.0./2.0e1)+1.0).*1.256954407538528e-1+in2(:,42).*(in2(:,42).*(1.0./2.0e1)+1.0).*4.574883722109476e-2+in2(:,43).*(in2(:,42).*(1.0./2.0e1)+1.0).*9.600403448055994e-3+in2(:,44).*(in2(:,42).*(1.0./2.0e1)+1.0).*1.49138709201216e-17-in2(:,45).*(in2(:,42).*(1.0./2.0e1)+1.0).*1.324322070428872e-1-in2(:,46).*(in2(:,42).*(1.0./2.0e1)+1.0).*7.563134468364148e-2-in2(:,47).*(in2(:,42).*(1.0./2.0e1)+1.0).*1.250903162787419e-1-in2(:,48).*(in2(:,42).*(1.0./2.0e1)+1.0).*1.174718703862174e-17+in2(:,53).*(in2(:,42).*(1.0./2.0e1)+1.0).*7.428229439131843e-2+in2(:,58).*(in2(:,42).*(1.0./2.0e1)+1.0).*6.862920708236914e-3);-erf(in2(:,43).*(3.0./2.0e1)-in2(:,10).*(in2(:,43).*(1.0./2.0e1)+1.0).*1.14832595399792e-1-in2(:,28).*(in2(:,43).*(1.0./2.0e1)+1.0).*1.048977743216548e-2+in2(:,41).*(in2(:,43).*(1.0./2.0e1)+1.0).*1.200559030541199e-2+in2(:,42).*(in2(:,43).*(1.0./2.0e1)+1.0).*1.241324178771654e-1+in2(:,43).*(in2(:,43).*(1.0./2.0e1)+1.0).*8.206870887089983e-2+in2(:,44).*(in2(:,43).*(1.0./2.0e1)+1.0).*1.048370045709639e-17-in2(:,45).*(in2(:,43).*(1.0./2.0e1)+1.0).*3.099912456904102e-2-in2(:,46).*(in2(:,43).*(1.0./2.0e1)+1.0).*4.086664217448681e-3-in2(:,47).*(in2(:,43).*(1.0./2.0e1)+1.0).*6.513172946871201e-2-in2(:,48).*(in2(:,43).*(1.0./2.0e1)+1.0).*6.91913896759581e-18);-erf(in2(:,44).*(3.0./2.0e1)+in2(:,41).*(in2(:,44).*(1.0./2.0e1)+1.0).*9.875419713421195e-18+in2(:,42).*(in2(:,44).*(1.0./2.0e1)+1.0).*1.904753693584705e-17+in2(:,43).*(in2(:,44).*(1.0./2.0e1)+1.0).*1.489724869261735e-17+in2(:,44).*(in2(:,44).*(1.0./2.0e1)+1.0).*1.118312837765146e-33-in2(:,45).*(in2(:,44).*(1.0./2.0e1)+1.0).*1.916701648170537e-19-in2(:,46).*(in2(:,44).*(1.0./2.0e1)+1.0).*6.038616430509341e-18-in2(:,47).*(in2(:,44).*(1.0./2.0e1)+1.0).*3.479744801542694e-18-in2(:,48).*(in2(:,44).*(1.0./2.0e1)+1.0).*4.325107168799558e-33);-erf(in2(:,45).*(3.0./2.0e1)-in2(:,22).*(in2(:,45).*(1.0./2.0e1)+1.0).*1.508546306977886e-5-in2(:,26).*(in2(:,45).*(1.0./2.0e1)+1.0).*8.401397348795517e-2+in2(:,34).*(in2(:,45).*(1.0./2.0e1)+1.0).*3.59903444190064e-3-in2(:,41).*(in2(:,45).*(1.0./2.0e1)+1.0).*8.048445818343347e-2-in2(:,42).*(in2(:,45).*(1.0./2.0e1)+1.0).*6.842781627840457e-2-in2(:,43).*(in2(:,45).*(1.0./2.0e1)+1.0).*9.148783853020702e-3-in2(:,44).*(in2(:,45).*(1.0./2.0e1)+1.0).*6.035991183173718e-18+in2(:,45).*(in2(:,45).*(1.0./2.0e1)+1.0).*1.492525690903476e-2+in2(:,46).*(in2(:,45).*(1.0./2.0e1)+1.0).*4.999859205014167e-2+in2(:,47).*(in2(:,45).*(1.0./2.0e1)+1.0).*8.333170257084657e-2+in2(:,48).*(in2(:,45).*(1.0./2.0e1)+1.0).*3.291191163330941e-17+in2(:,57).*(in2(:,45).*(1.0./2.0e1)+1.0).*7.740482991402331e-2);-erf(in2(:,46).*(3.0./2.0e1)-in2(:,13).*(in2(:,46).*(1.0./2.0e1)+1.0).*5.845385729369587e-2-in2(:,41).*(in2(:,46).*(1.0./2.0e1)+1.0).*2.067705882414607e-3-in2(:,42).*(in2(:,46).*(1.0./2.0e1)+1.0).*1.984027975455437e-1-in2(:,43).*(in2(:,46).*(1.0./2.0e1)+1.0).*1.094805199427041e-1-in2(:,44).*(in2(:,46).*(1.0./2.0e1)+1.0).*1.200563173507961e-17+in2(:,45).*(in2(:,46).*(1.0./2.0e1)+1.0).*5.924777256359193e-2+in2(:,46).*(in2(:,46).*(1.0./2.0e1)+1.0).*7.179177907980673e-2+in2(:,47).*(in2(:,46).*(1.0./2.0e1)+1.0).*5.119295430898408e-2+in2(:,48).*(in2(:,46).*(1.0./2.0e1)+1.0).*3.226226629891917e-17+in2(:,62).*(in2(:,46).*(1.0./2.0e1)+1.0).*9.851588164508956e-3);-erf(in2(:,47).*(3.0./2.0e1)-in2(:,23).*(in2(:,47).*(1.0./2.0e1)+1.0).*3.89553201491433e-2+in2(:,38).*(in2(:,47).*(1.0./2.0e1)+1.0).*3.278269121867893e-3-in2(:,41).*(in2(:,47).*(1.0./2.0e1)+1.0).*7.609735371371004e-2-in2(:,42).*(in2(:,47).*(1.0./2.0e1)+1.0).*9.029900647950076e-2-in2(:,43).*(in2(:,47).*(1.0./2.0e1)+1.0).*8.444232497949762e-2-in2(:,44).*(in2(:,47).*(1.0./2.0e1)+1.0).*7.998564105512524e-18+in2(:,45).*(in2(:,47).*(1.0./2.0e1)+1.0).*7.6610395743341e-2+in2(:,46).*(in2(:,47).*(1.0./2.0e1)+1.0).*1.373421605815184e-1+in2(:,47).*(in2(:,47).*(1.0./2.0e1)+1.0).*9.805367964540126e-2+in2(:,48).*(in2(:,47).*(1.0./2.0e1)+1.0).*2.736079102385573e-17+in2(:,51).*(in2(:,47).*(1.0./2.0e1)+1.0).*1.393166398515541e-1+in2(:,59).*(in2(:,47).*(1.0./2.0e1)+1.0).*5.880394343021392e-2+in2(:,62).*(in2(:,47).*(1.0./2.0e1)+1.0).*5.15075301537647e-3);-erf(in2(:,48).*(3.0./2.0e1)-in2(:,20).*(in2(:,48).*(1.0./2.0e1)+1.0).*3.721548782970641e-2-in2(:,41).*(in2(:,48).*(1.0./2.0e1)+1.0).*2.300870552576028e-17-in2(:,42).*(in2(:,48).*(1.0./2.0e1)+1.0).*3.700557830489578e-17-in2(:,43).*(in2(:,48).*(1.0./2.0e1)+1.0).*2.359865339314553e-17-in2(:,44).*(in2(:,48).*(1.0./2.0e1)+1.0).*4.931521723112676e-33+in2(:,45).*(in2(:,48).*(1.0./2.0e1)+1.0).*5.03990658323618e-18+in2(:,46).*(in2(:,48).*(1.0./2.0e1)+1.0).*4.331685002855733e-17+in2(:,47).*(in2(:,48).*(1.0./2.0e1)+1.0).*9.359755501529636e-18+in2(:,48).*(in2(:,48).*(1.0./2.0e1)+1.0).*6.525165061861911e-33+in2(:,58).*(in2(:,48).*(1.0./2.0e1)+1.0).*9.137424448763541e-2);-erf(in2(:,49).*(3.0./2.0e1)-in2(:,21).*(in2(:,49).*(1.0./2.0e1)+1.0).*9.857908176928684e-2+in2(:,49).*(in2(:,49).*(1.0./2.0e1)+1.0).*2.474320716763176e-2+in2(:,50).*(in2(:,49).*(1.0./2.0e1)+1.0).*7.176448535004341e-2+in2(:,51).*(in2(:,49).*(1.0./2.0e1)+1.0).*7.840362440233302e-2+in2(:,52).*(in2(:,49).*(1.0./2.0e1)+1.0).*9.708607316832438e-18-in2(:,53).*(in2(:,49).*(1.0./2.0e1)+1.0).*2.740790608493292e-2-in2(:,54).*(in2(:,49).*(1.0./2.0e1)+1.0).*8.071037279207702e-2-in2(:,55).*(in2(:,49).*(1.0./2.0e1)+1.0).*4.400643764021951e-2-in2(:,56).*(in2(:,49).*(1.0./2.0e1)+1.0).*1.701492775957736e-17+in2(:,64).*(in2(:,49).*(1.0./2.0e1)+1.0).*2.682242948835593e-17);-erf(in2(:,50).*(3.0./2.0e1)-in2(:,5).*(in2(:,50).*(1.0./2.0e1)+1.0).*6.894177068903391e-2-in2(:,16).*(in2(:,50).*(1.0./2.0e1)+1.0).*9.655729269920701e-2+in2(:,49).*(in2(:,50).*(1.0./2.0e1)+1.0).*1.240333657887553e-1+in2(:,50).*(in2(:,50).*(1.0./2.0e1)+1.0).*1.40751209192934e-1+in2(:,51).*(in2(:,50).*(1.0./2.0e1)+1.0).*6.722788406372257e-2+in2(:,52).*(in2(:,50).*(1.0./2.0e1)+1.0).*2.432596035326128e-17-in2(:,53).*(in2(:,50).*(1.0./2.0e1)+1.0).*1.207234500365136e-2-in2(:,54).*(in2(:,50).*(1.0./2.0e1)+1.0).*7.079620595304885e-2-in2(:,55).*(in2(:,50).*(1.0./2.0e1)+1.0).*1.699450962775171e-2-in2(:,56).*(in2(:,50).*(1.0./2.0e1)+1.0).*2.084075194454823e-17+in2(:,64).*(in2(:,50).*(1.0./2.0e1)+1.0).*6.512744097192034e-18);-erf(in2(:,51).*(3.0./2.0e1)-in2(:,10).*(in2(:,51).*(1.0./2.0e1)+1.0).*2.82435175833341e-2-in2(:,16).*(in2(:,51).*(1.0./2.0e1)+1.0).*2.983201911797166e-2+in2(:,49).*(in2(:,51).*(1.0./2.0e1)+1.0).*7.5564769016303e-2+in2(:,50).*(in2(:,51).*(1.0./2.0e1)+1.0).*8.450550562253886e-2+in2(:,51).*(in2(:,51).*(1.0./2.0e1)+1.0).*8.082120102969016e-2+in2(:,52).*(in2(:,51).*(1.0./2.0e1)+1.0).*5.686456808468156e-18-in2(:,53).*(in2(:,51).*(1.0./2.0e1)+1.0).*5.901044201182928e-2-in2(:,54).*(in2(:,51).*(1.0./2.0e1)+1.0).*7.231609634703651e-2-in2(:,55).*(in2(:,51).*(1.0./2.0e1)+1.0).*5.359494099714684e-2-in2(:,56).*(in2(:,51).*(1.0./2.0e1)+1.0).*2.034038302273041e-17);-erf(in2(:,52).*(3.0./2.0e1)-in2(:,15).*(in2(:,52).*(1.0./2.0e1)+1.0).*5.120902145517747e-2-in2(:,23).*(in2(:,52).*(1.0./2.0e1)+1.0).*5.11350616589288e-2-in2(:,26).*(in2(:,52).*(1.0./2.0e1)+1.0).*8.972551913020023e-2-in2(:,27).*(in2(:,52).*(1.0./2.0e1)+1.0).*3.179887245590516e-2+in2(:,46).*(in2(:,52).*(1.0./2.0e1)+1.0).*2.94220255305881e-2+in2(:,49).*(in2(:,52).*(1.0./2.0e1)+1.0).*6.492498697889435e-18+in2(:,50).*(in2(:,52).*(1.0./2.0e1)+1.0).*1.473061307072669e-17+in2(:,51).*(in2(:,52).*(1.0./2.0e1)+1.0).*7.558665057256497e-18+in2(:,52).*(in2(:,52).*(1.0./2.0e1)+1.0).*1.106700373479543e-33-in2(:,53).*(in2(:,52).*(1.0./2.0e1)+1.0).*3.676226223975795e-18-in2(:,54).*(in2(:,52).*(1.0./2.0e1)+1.0).*2.464050699156708e-17-in2(:,55).*(in2(:,52).*(1.0./2.0e1)+1.0).*3.499503617659552e-18-in2(:,56).*(in2(:,52).*(1.0./2.0e1)+1.0).*3.013036189826127e-33+in2(:,62).*(in2(:,52).*(1.0./2.0e1)+1.0).*2.263334307145609e-2);-erf(in2(:,53).*(3.0./2.0e1)-in2(:,49).*(in2(:,53).*(1.0./2.0e1)+1.0).*6.574207844711497e-2-in2(:,50).*(in2(:,53).*(1.0./2.0e1)+1.0).*1.351821608569858e-1-in2(:,51).*(in2(:,53).*(1.0./2.0e1)+1.0).*4.667503714255333e-2-in2(:,52).*(in2(:,53).*(1.0./2.0e1)+1.0).*4.614634407974253e-18+in2(:,53).*(in2(:,53).*(1.0./2.0e1)+1.0).*1.587358892880781e-2+in2(:,54).*(in2(:,53).*(1.0./2.0e1)+1.0).*3.190976663469747e-2+in2(:,55).*(in2(:,53).*(1.0./2.0e1)+1.0).*7.836033969199643e-2+in2(:,56).*(in2(:,53).*(1.0./2.0e1)+1.0).*1.328815694723028e-17);-erf(in2(:,54).*(3.0./2.0e1)-in2(:,17).*(in2(:,54).*(1.0./2.0e1)+1.0).*2.822047994137444e-3-in2(:,49).*(in2(:,54).*(1.0./2.0e1)+1.0).*1.200979997243381e-1-in2(:,50).*(in2(:,54).*(1.0./2.0e1)+1.0).*4.118986310356282e-2-in2(:,51).*(in2(:,54).*(1.0./2.0e1)+1.0).*4.366315958105379e-2-in2(:,52).*(in2(:,54).*(1.0./2.0e1)+1.0).*1.027839725995038e-17+in2(:,53).*(in2(:,54).*(1.0./2.0e1)+1.0).*1.04856327276143e-1+in2(:,54).*(in2(:,54).*(1.0./2.0e1)+1.0).*5.948713840397096e-2+in2(:,55).*(in2(:,54).*(1.0./2.0e1)+1.0).*1.157346826256587e-1+in2(:,56).*(in2(:,54).*(1.0./2.0e1)+1.0).*2.848246108862831e-17+in2(:,63).*(in2(:,54).*(1.0./2.0e1)+1.0).*1.267821129429252e-2);-erf(in2(:,55).*(3.0./2.0e1)-in2(:,10).*(in2(:,55).*(1.0./2.0e1)+1.0).*4.908391084411631e-2-in2(:,26).*(in2(:,55).*(1.0./2.0e1)+1.0).*4.036054488878428e-2+in2(:,40).*(in2(:,55).*(1.0./2.0e1)+1.0).*2.483145586383008e-2+in2(:,47).*(in2(:,55).*(1.0./2.0e1)+1.0).*8.674241580784237e-2-in2(:,49).*(in2(:,55).*(1.0./2.0e1)+1.0).*8.743785113741707e-2-in2(:,50).*(in2(:,55).*(1.0./2.0e1)+1.0).*1.265517588410084e-1-in2(:,51).*(in2(:,55).*(1.0./2.0e1)+1.0).*7.151151727964981e-2-in2(:,52).*(in2(:,55).*(1.0./2.0e1)+1.0).*2.316878230497685e-18+in2(:,53).*(in2(:,55).*(1.0./2.0e1)+1.0).*1.850275844432574e-2+in2(:,54).*(in2(:,55).*(1.0./2.0e1)+1.0).*1.310801084549697e-2+in2(:,55).*(in2(:,55).*(1.0./2.0e1)+1.0).*7.978518057024689e-3+in2(:,56).*(in2(:,55).*(1.0./2.0e1)+1.0).*1.499559736614852e-17);-erf(in2(:,56).*(3.0./2.0e1)-in2(:,1).*(in2(:,56).*(1.0./2.0e1)+1.0).*1.119927463685712e-2-in2(:,22).*(in2(:,56).*(1.0./2.0e1)+1.0).*7.838964882830754e-2+in2(:,42).*(in2(:,56).*(1.0./2.0e1)+1.0).*4.115872181970816e-2-in2(:,49).*(in2(:,56).*(1.0./2.0e1)+1.0).*3.02652120010873e-17-in2(:,50).*(in2(:,56).*(1.0./2.0e1)+1.0).*4.025622747450175e-17-in2(:,51).*(in2(:,56).*(1.0./2.0e1)+1.0).*1.240029476641861e-17-in2(:,52).*(in2(:,56).*(1.0./2.0e1)+1.0).*4.551656614790261e-33+in2(:,53).*(in2(:,56).*(1.0./2.0e1)+1.0).*5.619634183728522e-18+in2(:,54).*(in2(:,56).*(1.0./2.0e1)+1.0).*3.443000973625418e-17+in2(:,55).*(in2(:,56).*(1.0./2.0e1)+1.0).*1.40060539345841e-17+in2(:,56).*(in2(:,56).*(1.0./2.0e1)+1.0).*1.512825254661356e-33);-erf(in2(:,57).*(3.0./2.0e1)-in2(:,11).*(in2(:,57).*(1.0./2.0e1)+1.0).*7.480796732804948e-2+in2(:,34).*(in2(:,57).*(1.0./2.0e1)+1.0).*1.825312539812768e-2+in2(:,44).*(in2(:,57).*(1.0./2.0e1)+1.0).*6.501515728295869e-2+in2(:,57).*(in2(:,57).*(1.0./2.0e1)+1.0).*1.295090709595031e-2+in2(:,58).*(in2(:,57).*(1.0./2.0e1)+1.0).*1.592452335141693e-2+in2(:,59).*(in2(:,57).*(1.0./2.0e1)+1.0).*9.529129822350146e-2+in2(:,60).*(in2(:,57).*(1.0./2.0e1)+1.0).*1.00353789261408e-17-in2(:,61).*(in2(:,57).*(1.0./2.0e1)+1.0).*2.783888204306757e-2-in2(:,62).*(in2(:,57).*(1.0./2.0e1)+1.0).*2.703593297723952e-2-in2(:,63).*(in2(:,57).*(1.0./2.0e1)+1.0).*4.101221800476773e-2-in2(:,64).*(in2(:,57).*(1.0./2.0e1)+1.0).*1.569640112971119e-17);-erf(in2(:,58).*(3.0./2.0e1)-in2(:,11).*(in2(:,58).*(1.0./2.0e1)+1.0).*1.031376131926105e-1-in2(:,31).*(in2(:,58).*(1.0./2.0e1)+1.0).*5.193686421285055e-3+in2(:,38).*(in2(:,58).*(1.0./2.0e1)+1.0).*1.765498087376976e-2+in2(:,57).*(in2(:,58).*(1.0./2.0e1)+1.0).*3.943088545481338e-2+in2(:,58).*(in2(:,58).*(1.0./2.0e1)+1.0).*1.388759041253312e-1+in2(:,59).*(in2(:,58).*(1.0./2.0e1)+1.0).*9.778193073518422e-2+in2(:,60).*(in2(:,58).*(1.0./2.0e1)+1.0).*1.966751563309325e-17-in2(:,61).*(in2(:,58).*(1.0./2.0e1)+1.0).*1.242733998673855e-1-in2(:,62).*(in2(:,58).*(1.0./2.0e1)+1.0).*1.702677527359229e-1-in2(:,63).*(in2(:,58).*(1.0./2.0e1)+1.0).*1.061432417260699e-1-in2(:,64).*(in2(:,58).*(1.0./2.0e1)+1.0).*3.84598986562603e-18);-erf(in2(:,59).*(3.0./2.0e1)-in2(:,3).*(in2(:,59).*(1.0./2.0e1)+1.0).*2.489155436474396e-2-in2(:,5).*(in2(:,59).*(1.0./2.0e1)+1.0).*3.72233822660077e-2-in2(:,9).*(in2(:,59).*(1.0./2.0e1)+1.0).*3.329603480113025e-2-in2(:,18).*(in2(:,59).*(1.0./2.0e1)+1.0).*6.061570163718779e-3-in2(:,24).*(in2(:,59).*(1.0./2.0e1)+1.0).*2.012549477423741e-3-in2(:,25).*(in2(:,59).*(1.0./2.0e1)+1.0).*1.613863336832132e-3+in2(:,45).*(in2(:,59).*(1.0./2.0e1)+1.0).*9.299671414362611e-2+in2(:,54).*(in2(:,59).*(1.0./2.0e1)+1.0).*3.099360405880699e-3+in2(:,57).*(in2(:,59).*(1.0./2.0e1)+1.0).*4.233631888608366e-2+in2(:,58).*(in2(:,59).*(1.0./2.0e1)+1.0).*1.408758829349345e-1+in2(:,59).*(in2(:,59).*(1.0./2.0e1)+1.0).*8.460646499591547e-2+in2(:,60).*(in2(:,59).*(1.0./2.0e1)+1.0).*4.899019278472999e-18-in2(:,61).*(in2(:,59).*(1.0./2.0e1)+1.0).*3.000310493726665e-2-in2(:,62).*(in2(:,59).*(1.0./2.0e1)+1.0).*1.08481684777812e-1-in2(:,63).*(in2(:,59).*(1.0./2.0e1)+1.0).*3.355046432960932e-2-in2(:,64).*(in2(:,59).*(1.0./2.0e1)+1.0).*2.703991288438674e-18);-erf(in2(:,60).*(3.0./2.0e1)+in2(:,38).*(in2(:,60).*(1.0./2.0e1)+1.0).*1.796485504983708e-2+in2(:,57).*(in2(:,60).*(1.0./2.0e1)+1.0).*1.455641521909379e-17+in2(:,58).*(in2(:,60).*(1.0./2.0e1)+1.0).*4.259849100331935e-18+in2(:,59).*(in2(:,60).*(1.0./2.0e1)+1.0).*1.401742409502955e-17+in2(:,60).*(in2(:,60).*(1.0./2.0e1)+1.0).*2.829061481755543e-33-in2(:,61).*(in2(:,60).*(1.0./2.0e1)+1.0).*7.930137966830369e-18-in2(:,62).*(in2(:,60).*(1.0./2.0e1)+1.0).*8.71555130215685e-18-in2(:,63).*(in2(:,60).*(1.0./2.0e1)+1.0).*6.949053962336861e-18-in2(:,64).*(in2(:,60).*(1.0./2.0e1)+1.0).*3.674262933383157e-33);-erf(in2(:,61).*(3.0./2.0e1)-in2(:,16).*(in2(:,61).*(1.0./2.0e1)+1.0).*4.066636233209656e-2+in2(:,49).*(in2(:,61).*(1.0./2.0e1)+1.0).*3.492924098609727e-2+in2(:,54).*(in2(:,61).*(1.0./2.0e1)+1.0).*1.381969337248854e-2-in2(:,57).*(in2(:,61).*(1.0./2.0e1)+1.0).*4.534221575590793e-2-in2(:,58).*(in2(:,61).*(1.0./2.0e1)+1.0).*1.209103421023022e-1-in2(:,59).*(in2(:,61).*(1.0./2.0e1)+1.0).*6.864835913792819e-2-in2(:,60).*(in2(:,61).*(1.0./2.0e1)+1.0).*1.126758477806807e-17+in2(:,61).*(in2(:,61).*(1.0./2.0e1)+1.0).*1.66750879961244e-2+in2(:,62).*(in2(:,61).*(1.0./2.0e1)+1.0).*1.004549669841348e-1+in2(:,63).*(in2(:,61).*(1.0./2.0e1)+1.0).*2.051734931790909e-2+in2(:,64).*(in2(:,61).*(1.0./2.0e1)+1.0).*2.710679050765018e-18);-erf(in2(:,62).*(3.0./2.0e1)-in2(:,18).*(in2(:,62).*(1.0./2.0e1)+1.0).*1.78619330592731e-2-in2(:,25).*(in2(:,62).*(1.0./2.0e1)+1.0).*1.943190076796052e-2-in2(:,57).*(in2(:,62).*(1.0./2.0e1)+1.0).*1.152545097597477e-1-in2(:,58).*(in2(:,62).*(1.0./2.0e1)+1.0).*7.870461430853969e-2-in2(:,59).*(in2(:,62).*(1.0./2.0e1)+1.0).*5.535412166204887e-2-in2(:,60).*(in2(:,62).*(1.0./2.0e1)+1.0).*1.825113928338218e-18+in2(:,61).*(in2(:,62).*(1.0./2.0e1)+1.0).*7.650001022154299e-2+in2(:,62).*(in2(:,62).*(1.0./2.0e1)+1.0).*6.291043604403004e-2+in2(:,63).*(in2(:,62).*(1.0./2.0e1)+1.0).*2.420489176686986e-2+in2(:,64).*(in2(:,62).*(1.0./2.0e1)+1.0).*7.170186978859028e-18);-erf(in2(:,63).*(3.0./2.0e1)-in2(:,11).*(in2(:,63).*(1.0./2.0e1)+1.0).*2.355904578670434e-3-in2(:,26).*(in2(:,63).*(1.0./2.0e1)+1.0).*1.003127522218053e-2-in2(:,57).*(in2(:,63).*(1.0./2.0e1)+1.0).*2.551002620247973e-2-in2(:,58).*(in2(:,63).*(1.0./2.0e1)+1.0).*9.740784677395445e-2-in2(:,59).*(in2(:,63).*(1.0./2.0e1)+1.0).*7.052283304073764e-2-in2(:,60).*(in2(:,63).*(1.0./2.0e1)+1.0).*1.465651950645551e-17+in2(:,61).*(in2(:,63).*(1.0./2.0e1)+1.0).*3.802147402723738e-2+in2(:,62).*(in2(:,63).*(1.0./2.0e1)+1.0).*7.320238176740972e-2+in2(:,63).*(in2(:,63).*(1.0./2.0e1)+1.0).*4.396268079917013e-2+in2(:,64).*(in2(:,63).*(1.0./2.0e1)+1.0).*3.444352548681069e-17);-erf(in2(:,64).*(3.0./2.0e1)-in2(:,15).*(in2(:,64).*(1.0./2.0e1)+1.0).*2.856001192225517e-17+in2(:,44).*(in2(:,64).*(1.0./2.0e1)+1.0).*2.140220335099263e-17+in2(:,52).*(in2(:,64).*(1.0./2.0e1)+1.0).*1.701109246254243e-17-in2(:,57).*(in2(:,64).*(1.0./2.0e1)+1.0).*2.530591542930475e-17-in2(:,58).*(in2(:,64).*(1.0./2.0e1)+1.0).*3.102390525648194e-17-in2(:,59).*(in2(:,64).*(1.0./2.0e1)+1.0).*2.772392850656886e-17-in2(:,60).*(in2(:,64).*(1.0./2.0e1)+1.0).*4.507862019525818e-33+in2(:,61).*(in2(:,64).*(1.0./2.0e1)+1.0).*2.911662616791729e-17+in2(:,62).*(in2(:,64).*(1.0./2.0e1)+1.0).*2.892478572723488e-17+in2(:,63).*(in2(:,64).*(1.0./2.0e1)+1.0).*3.471311068691882e-18+in2(:,64).*(in2(:,64).*(1.0./2.0e1)+1.0).*1.042571923704379e-32)]

Error in odearguments (line 87)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

Error in ode45 (line 113)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

[Walter Roberson]

odefun = matlabFunction(rhs(:), 'vars', {sym('t') X.'});

However, I find that in practice the outputs are linear and dependent upon init, so with the all-ones init you are passing, the outputs are all the same.

It is not clear to me at the moment why the outputs should all turn out the same.

[Walter Roberson]

It might have been the random Q that I choose for testing purposes. All of the results were above 250 before the erf(), and erf() of that is 1 to within working position...

Yes, a different (much smaller) Q made a notable difference.

[holistic]

Thank you once again Walter!

This works now. With regard to the initial conditions in my case, I just choose some arbitrary ones. But normally, the initial conditions won't be the same for all "units".