Category: Matlab
Category Archives: Matlab
I get this error: Incorrect number or types of inputs or outputs for function ‘phaseShift’ when I try to compile my code.
Hi everyone. Could someone please explain what could be the cause of this error when I use the phaseShift function. Many thanks in advance.
%%Create Variables
clc
clear all
close all
freq = 5.33e9;
c = physconst("lightspeed");
lambda = c/freq;
%%Create pcbStack object
pcbobj = pcbStack;
%%Create board shape
%Creating BoardShape metal layer.
%Creating Rectangle shape.
Rectangle = antenna.Rectangle;
Rectangle.Name = "Rectangle";
Rectangle.Center = [0 0];
Rectangle.Length = 0.02725;
Rectangle.Width = 0.02725;
Rectangle = rotate(Rectangle,0,[Rectangle.Center,-1],[Rectangle.Center,1]);
BoardShape = Rectangle;
pcbobj.BoardShape = BoardShape;
%%Create Stackup
%Creating PatchLayer metal layer.
%Creating Patch shape.
Patch = antenna.Polygon;
Patch.Name = "Patch";
Patch.Vertices = [
-0.00775 0.0064 0;
0.00775 0.0064 0;
0.00775 -0.0064 0;
0.000535 -0.0064 0;
0.000535 -0.0015 0;
0.000375 -0.0015 0;
0.000375 -0.013625 0;
-0.000375 -0.013625 0;
-0.000375 -0.0015 0;
-0.000535 -0.0015 0;
-0.000535 -0.0064 0;
-0.00775 -0.0064 0];
Patch = rotate(Patch,0,[0 -0.00430416666666667 -1],[0 -0.00430416666666667 1]);
PatchLayer = Patch;
%Creating DielectricLayer dielectric layer.
DielectricLayer = dielectric("Name",’Custom’,"EpsilonR",4.66,"LossTangent",0.026,"Thickness",0.0006);
%Creating GroundPlane metal layer.
%Creating Ground shape.
Ground = antenna.Rectangle;
Ground.Name = "Ground";
Ground.Center = [0 0];
Ground.Length = 0.02725;
Ground.Width = 0.02725;
Ground = rotate(Ground,0,[Ground.Center,-1],[Ground.Center,1]);
GroundPlane = Ground;
%%Create Feed
% feedloc = [[0 -0.006237 3 1];…
% ];
feedloc = [[0 -0.0064 3 1];…
];
%%Create Metal
metalobj = metal;
metalobj.Name = ‘Copper’;
metalobj.Conductivity = 5.96e7;
metalobj.Thickness = 35e-6;
pcbobj.Conductor = metalobj;
%%Assign propertie
pcbobj.BoardThickness = 0.00053;
pcbobj.Layers = {PatchLayer,DielectricLayer,GroundPlane,};
pcbobj.FeedLocations = feedloc;
pcbobj.FeedDiameter = 0.00030672;
pcbobj.ViaDiameter = 0.00075;
pcbobj.FeedViaModel = ‘strip’;
pcbobj.FeedVoltage = 1;
pcbobj.FeedPhase = 0;
figure(‘Position’, [100 750 600 500]);
show(pcbobj)
figure(‘Position’, [700 750 600 500]);
pattern(pcbobj,freq)
figure(‘Position’, [1300 750 600 500]);
current(pcbobj,freq)
% antennaArray = rectangularArray;
antennaArray = phased.URA;
antennaArray.Element = pcbobj;
antennaArray.Size = [10 10];
antennaArray.ElementSpacing = lambda/2;
antennaArray.Lattice = "Rectangular";
antennaArray.ArrayNormal = "x";
figure(‘Position’, [100 150 600 500]);
pattern(antennaArray, freq)
pos_rec = getElementPosition(antennaArray);
pos_yx_rec = pos_rec([2,3],:);
figure(‘Position’, [700 150 600 500]);
plot(pos_yx_rec(1,:), pos_yx_rec(2,:), ‘.’, ‘MarkerSize’, 10)
% axis([-0.2 0.2 -0.2 0.2])
xlabel(‘x’); ylabel(‘y’)
title(‘Rectangular Lattice’)
figure(‘Position’, [1300 150 600 500]);
az_angle = -180:0.25:180;
pattern(antennaArray,freq, az_angle, 0, CoordinateSystem="rectangular")
axis([-90 90 -35 25])
ps = phaseShift(antennaArray, freq, [30 45]);Hi everyone. Could someone please explain what could be the cause of this error when I use the phaseShift function. Many thanks in advance.
%%Create Variables
clc
clear all
close all
freq = 5.33e9;
c = physconst("lightspeed");
lambda = c/freq;
%%Create pcbStack object
pcbobj = pcbStack;
%%Create board shape
%Creating BoardShape metal layer.
%Creating Rectangle shape.
Rectangle = antenna.Rectangle;
Rectangle.Name = "Rectangle";
Rectangle.Center = [0 0];
Rectangle.Length = 0.02725;
Rectangle.Width = 0.02725;
Rectangle = rotate(Rectangle,0,[Rectangle.Center,-1],[Rectangle.Center,1]);
BoardShape = Rectangle;
pcbobj.BoardShape = BoardShape;
%%Create Stackup
%Creating PatchLayer metal layer.
%Creating Patch shape.
Patch = antenna.Polygon;
Patch.Name = "Patch";
Patch.Vertices = [
-0.00775 0.0064 0;
0.00775 0.0064 0;
0.00775 -0.0064 0;
0.000535 -0.0064 0;
0.000535 -0.0015 0;
0.000375 -0.0015 0;
0.000375 -0.013625 0;
-0.000375 -0.013625 0;
-0.000375 -0.0015 0;
-0.000535 -0.0015 0;
-0.000535 -0.0064 0;
-0.00775 -0.0064 0];
Patch = rotate(Patch,0,[0 -0.00430416666666667 -1],[0 -0.00430416666666667 1]);
PatchLayer = Patch;
%Creating DielectricLayer dielectric layer.
DielectricLayer = dielectric("Name",’Custom’,"EpsilonR",4.66,"LossTangent",0.026,"Thickness",0.0006);
%Creating GroundPlane metal layer.
%Creating Ground shape.
Ground = antenna.Rectangle;
Ground.Name = "Ground";
Ground.Center = [0 0];
Ground.Length = 0.02725;
Ground.Width = 0.02725;
Ground = rotate(Ground,0,[Ground.Center,-1],[Ground.Center,1]);
GroundPlane = Ground;
%%Create Feed
% feedloc = [[0 -0.006237 3 1];…
% ];
feedloc = [[0 -0.0064 3 1];…
];
%%Create Metal
metalobj = metal;
metalobj.Name = ‘Copper’;
metalobj.Conductivity = 5.96e7;
metalobj.Thickness = 35e-6;
pcbobj.Conductor = metalobj;
%%Assign propertie
pcbobj.BoardThickness = 0.00053;
pcbobj.Layers = {PatchLayer,DielectricLayer,GroundPlane,};
pcbobj.FeedLocations = feedloc;
pcbobj.FeedDiameter = 0.00030672;
pcbobj.ViaDiameter = 0.00075;
pcbobj.FeedViaModel = ‘strip’;
pcbobj.FeedVoltage = 1;
pcbobj.FeedPhase = 0;
figure(‘Position’, [100 750 600 500]);
show(pcbobj)
figure(‘Position’, [700 750 600 500]);
pattern(pcbobj,freq)
figure(‘Position’, [1300 750 600 500]);
current(pcbobj,freq)
% antennaArray = rectangularArray;
antennaArray = phased.URA;
antennaArray.Element = pcbobj;
antennaArray.Size = [10 10];
antennaArray.ElementSpacing = lambda/2;
antennaArray.Lattice = "Rectangular";
antennaArray.ArrayNormal = "x";
figure(‘Position’, [100 150 600 500]);
pattern(antennaArray, freq)
pos_rec = getElementPosition(antennaArray);
pos_yx_rec = pos_rec([2,3],:);
figure(‘Position’, [700 150 600 500]);
plot(pos_yx_rec(1,:), pos_yx_rec(2,:), ‘.’, ‘MarkerSize’, 10)
% axis([-0.2 0.2 -0.2 0.2])
xlabel(‘x’); ylabel(‘y’)
title(‘Rectangular Lattice’)
figure(‘Position’, [1300 150 600 500]);
az_angle = -180:0.25:180;
pattern(antennaArray,freq, az_angle, 0, CoordinateSystem="rectangular")
axis([-90 90 -35 25])
ps = phaseShift(antennaArray, freq, [30 45]); Hi everyone. Could someone please explain what could be the cause of this error when I use the phaseShift function. Many thanks in advance.
%%Create Variables
clc
clear all
close all
freq = 5.33e9;
c = physconst("lightspeed");
lambda = c/freq;
%%Create pcbStack object
pcbobj = pcbStack;
%%Create board shape
%Creating BoardShape metal layer.
%Creating Rectangle shape.
Rectangle = antenna.Rectangle;
Rectangle.Name = "Rectangle";
Rectangle.Center = [0 0];
Rectangle.Length = 0.02725;
Rectangle.Width = 0.02725;
Rectangle = rotate(Rectangle,0,[Rectangle.Center,-1],[Rectangle.Center,1]);
BoardShape = Rectangle;
pcbobj.BoardShape = BoardShape;
%%Create Stackup
%Creating PatchLayer metal layer.
%Creating Patch shape.
Patch = antenna.Polygon;
Patch.Name = "Patch";
Patch.Vertices = [
-0.00775 0.0064 0;
0.00775 0.0064 0;
0.00775 -0.0064 0;
0.000535 -0.0064 0;
0.000535 -0.0015 0;
0.000375 -0.0015 0;
0.000375 -0.013625 0;
-0.000375 -0.013625 0;
-0.000375 -0.0015 0;
-0.000535 -0.0015 0;
-0.000535 -0.0064 0;
-0.00775 -0.0064 0];
Patch = rotate(Patch,0,[0 -0.00430416666666667 -1],[0 -0.00430416666666667 1]);
PatchLayer = Patch;
%Creating DielectricLayer dielectric layer.
DielectricLayer = dielectric("Name",’Custom’,"EpsilonR",4.66,"LossTangent",0.026,"Thickness",0.0006);
%Creating GroundPlane metal layer.
%Creating Ground shape.
Ground = antenna.Rectangle;
Ground.Name = "Ground";
Ground.Center = [0 0];
Ground.Length = 0.02725;
Ground.Width = 0.02725;
Ground = rotate(Ground,0,[Ground.Center,-1],[Ground.Center,1]);
GroundPlane = Ground;
%%Create Feed
% feedloc = [[0 -0.006237 3 1];…
% ];
feedloc = [[0 -0.0064 3 1];…
];
%%Create Metal
metalobj = metal;
metalobj.Name = ‘Copper’;
metalobj.Conductivity = 5.96e7;
metalobj.Thickness = 35e-6;
pcbobj.Conductor = metalobj;
%%Assign propertie
pcbobj.BoardThickness = 0.00053;
pcbobj.Layers = {PatchLayer,DielectricLayer,GroundPlane,};
pcbobj.FeedLocations = feedloc;
pcbobj.FeedDiameter = 0.00030672;
pcbobj.ViaDiameter = 0.00075;
pcbobj.FeedViaModel = ‘strip’;
pcbobj.FeedVoltage = 1;
pcbobj.FeedPhase = 0;
figure(‘Position’, [100 750 600 500]);
show(pcbobj)
figure(‘Position’, [700 750 600 500]);
pattern(pcbobj,freq)
figure(‘Position’, [1300 750 600 500]);
current(pcbobj,freq)
% antennaArray = rectangularArray;
antennaArray = phased.URA;
antennaArray.Element = pcbobj;
antennaArray.Size = [10 10];
antennaArray.ElementSpacing = lambda/2;
antennaArray.Lattice = "Rectangular";
antennaArray.ArrayNormal = "x";
figure(‘Position’, [100 150 600 500]);
pattern(antennaArray, freq)
pos_rec = getElementPosition(antennaArray);
pos_yx_rec = pos_rec([2,3],:);
figure(‘Position’, [700 150 600 500]);
plot(pos_yx_rec(1,:), pos_yx_rec(2,:), ‘.’, ‘MarkerSize’, 10)
% axis([-0.2 0.2 -0.2 0.2])
xlabel(‘x’); ylabel(‘y’)
title(‘Rectangular Lattice’)
figure(‘Position’, [1300 150 600 500]);
az_angle = -180:0.25:180;
pattern(antennaArray,freq, az_angle, 0, CoordinateSystem="rectangular")
axis([-90 90 -35 25])
ps = phaseShift(antennaArray, freq, [30 45]); phaseshift MATLAB Answers — New Questions
Mex error while running a code matlab 2021a
Hi, i’m new to matlab so i can’t figure out where the error is. Could somone tell me how I can rectify it
I’m using this code which sombody had already made.
Thanks in advance
function Create_Data()
% This function creates the model data depicted in Figure 4-6 of
%
% Characterization of K-Complexes and Slow Wave Activity in a Neural Mass Model
% A Weigenand, M Schellenberger Costa, H-VV Ngo, JC Claussen, T Martinetz
% PLoS Computational Biology. 2014;10:e1003923
% Move to source folder(assuming it contains the Figures folder
cd ..;
% Check if the executable exists and compile if needed
if(exist(‘Cortex_mex.mesa64’, ‘file’)==0)
mex CXXFLAGS="$CXXFLAGS -std=c++11 -O3" Cortex_mex.cpp Cortical_Column.cpp;
end
% Add the path to the simulation routine
addpath(pwd);
% Go back into figures folder
cd Figures;
% Time series
Data_Time_Series(1);
Data_Time_Series(2);
% Time series
Data_Stimulation();
end
First I got this error
>> Create_Data
Building with ‘MinGW64 Compiler (C++)’.
Error using mex
g++: error: -fexceptions: No such file or directory
So, I put D:MatlabFilescortex_mex.cpp
>> D:MatlabFilescortex_mex.cpp
D:MatlabFilescortex_mex.cpp
↑
Invalid use of operator.
But it still said thisHi, i’m new to matlab so i can’t figure out where the error is. Could somone tell me how I can rectify it
I’m using this code which sombody had already made.
Thanks in advance
function Create_Data()
% This function creates the model data depicted in Figure 4-6 of
%
% Characterization of K-Complexes and Slow Wave Activity in a Neural Mass Model
% A Weigenand, M Schellenberger Costa, H-VV Ngo, JC Claussen, T Martinetz
% PLoS Computational Biology. 2014;10:e1003923
% Move to source folder(assuming it contains the Figures folder
cd ..;
% Check if the executable exists and compile if needed
if(exist(‘Cortex_mex.mesa64’, ‘file’)==0)
mex CXXFLAGS="$CXXFLAGS -std=c++11 -O3" Cortex_mex.cpp Cortical_Column.cpp;
end
% Add the path to the simulation routine
addpath(pwd);
% Go back into figures folder
cd Figures;
% Time series
Data_Time_Series(1);
Data_Time_Series(2);
% Time series
Data_Stimulation();
end
First I got this error
>> Create_Data
Building with ‘MinGW64 Compiler (C++)’.
Error using mex
g++: error: -fexceptions: No such file or directory
So, I put D:MatlabFilescortex_mex.cpp
>> D:MatlabFilescortex_mex.cpp
D:MatlabFilescortex_mex.cpp
↑
Invalid use of operator.
But it still said this Hi, i’m new to matlab so i can’t figure out where the error is. Could somone tell me how I can rectify it
I’m using this code which sombody had already made.
Thanks in advance
function Create_Data()
% This function creates the model data depicted in Figure 4-6 of
%
% Characterization of K-Complexes and Slow Wave Activity in a Neural Mass Model
% A Weigenand, M Schellenberger Costa, H-VV Ngo, JC Claussen, T Martinetz
% PLoS Computational Biology. 2014;10:e1003923
% Move to source folder(assuming it contains the Figures folder
cd ..;
% Check if the executable exists and compile if needed
if(exist(‘Cortex_mex.mesa64’, ‘file’)==0)
mex CXXFLAGS="$CXXFLAGS -std=c++11 -O3" Cortex_mex.cpp Cortical_Column.cpp;
end
% Add the path to the simulation routine
addpath(pwd);
% Go back into figures folder
cd Figures;
% Time series
Data_Time_Series(1);
Data_Time_Series(2);
% Time series
Data_Stimulation();
end
First I got this error
>> Create_Data
Building with ‘MinGW64 Compiler (C++)’.
Error using mex
g++: error: -fexceptions: No such file or directory
So, I put D:MatlabFilescortex_mex.cpp
>> D:MatlabFilescortex_mex.cpp
D:MatlabFilescortex_mex.cpp
↑
Invalid use of operator.
But it still said this mex compiler, mex error MATLAB Answers — New Questions
Protected model fails to build
Hello,
I have created a simulink model which implements an inlined s-function which uses a .tlc file and external .c and.h files. I have configured everything correctly as the simulation is working, and when I perform a compile this also goes without issues.
However, when a try to save te model as a protected model the compilation fails as it cannot fild the external .c and .h files. Event hough i have configured it correctly under "custom code" in the model settings.
I would expect that saving it as a protected model would also use these settings, but it does not seem to be the case. What am i dowing wrong?Hello,
I have created a simulink model which implements an inlined s-function which uses a .tlc file and external .c and.h files. I have configured everything correctly as the simulation is working, and when I perform a compile this also goes without issues.
However, when a try to save te model as a protected model the compilation fails as it cannot fild the external .c and .h files. Event hough i have configured it correctly under "custom code" in the model settings.
I would expect that saving it as a protected model would also use these settings, but it does not seem to be the case. What am i dowing wrong? Hello,
I have created a simulink model which implements an inlined s-function which uses a .tlc file and external .c and.h files. I have configured everything correctly as the simulation is working, and when I perform a compile this also goes without issues.
However, when a try to save te model as a protected model the compilation fails as it cannot fild the external .c and .h files. Event hough i have configured it correctly under "custom code" in the model settings.
I would expect that saving it as a protected model would also use these settings, but it does not seem to be the case. What am i dowing wrong? protected model, simulink MATLAB Answers — New Questions
Can we use custom markers to plot trajectory
Is there a way I can use a figure as custom marker to represent the trajectory of an unamanned aerial system. I also need it to represent the orientation of the UAS at any point on the trajectory.
I have attached examples that I found in some papers.Is there a way I can use a figure as custom marker to represent the trajectory of an unamanned aerial system. I also need it to represent the orientation of the UAS at any point on the trajectory.
I have attached examples that I found in some papers. Is there a way I can use a figure as custom marker to represent the trajectory of an unamanned aerial system. I also need it to represent the orientation of the UAS at any point on the trajectory.
I have attached examples that I found in some papers. figure, plot, plotting MATLAB Answers — New Questions
How to connect simpowersystem generic battery with simscape blocks and solve algebraic loop er?
Hello,
I’m new to Matlab-Simulink. I cannot directly connect the generic battery block to the simelectronics blocks in the model. To achieve this, I first connected the voltage measurement block to the generic battery terminals, and then its output signal to the controlled voltage sensor block. However, the voltage I receive this way is the initial voltage value of the generic battery and does not decrease over time depending on the load. To ensure this voltage decrease, I connected a controlled current source as a load to the generic battery. I connected the signal I received from the current sensor on the Simscape side to the input of the current source. But this gave the following algebraic loop error.
Is there a different way to connect the generic battery block to the simscape blocks without getting an algebraic loop error?
Thanks for any help.
Best Regards
Block diagram ‘DronePropellerMotorBattery_v5’ contains 1 algebraic loop(s). To see more details about the loops use the command Simulink.BlockDiagram.getAlgebraicLoops(‘DronePropellerMotorBattery_v5’) or the command line Simulink debugger by typing sldebug(‘DronePropellerMotorBattery_v5’) in the MATLAB command window. To eliminate this message, set Algebraic loop to "none".
Component:Simulink | Category:Block diagram warning
Found algebraic loop containing: DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/INPUT_1_1_1DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/STATE_1 (discontinuity) DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/OUTPUT_1_0DronePropellerMotorBattery_v5/powergui/EquivalentModel1/State-SpaceDronePropellerMotorBattery_v5/Voltage Measurement/do not delete this gain (algebraic variable)
Component:Simulink | Category:Model
Discontinuities detected within algebraic loop(s), may have trouble solving
Component:Simulink | Category:Model warning
Convergence problem when solving algebraic loop containing ‘DronePropellerMotorBattery_v5/Solver Configuration’ at time 0.0. Simulink will try to solve this loop using Simulink 3 (R11) strategy. Use feature(‘ModeIterationsInAlgLoops’,0) to disable the strategy introduced in Simulink 4 (R12)
Component:Simulink | Category:Block warning
Solver is encountering difficulty in simulating model ‘DronePropellerMotorBattery_v5’ at time 1.1125369292536007e-308. Simulink will continue to simulate with warnings. Please check the model for errors.
Caused by:
Solver was unable to reduce the step size without violating the minimum step size value of 2.22507E-308 for 1 consecutive times at time 0. This usually may be caused by violating algebraic constraints in the differential-algebraic system or by the high stiffness of the system. Try tightening the error tolerances, and/or the tolerances for computing consistent conditions. If the problem persists, please check the system or increase the solver Number of consecutive min steps violation parameter.
Suggested Actions
Open Solver Profiler to explore such issuesOpen
Component:Simulink | Category:Block diagram warning
An error occurred during simulation and the simulation was terminated
Caused by:
Algebraic state in algebraic loop containing ‘DronePropellerMotorBattery_v5/Voltage Measurement/do not delete this gain’ computed at time 2.2250738585072014E-308 is Inf or NaN. There may be a singularity in the solution. If the model is correct, try reducing the step size (either by reducing the fixed step size or by tightening the error tolerances) or tweak the initial guess of algebraic loop variable values.Hello,
I’m new to Matlab-Simulink. I cannot directly connect the generic battery block to the simelectronics blocks in the model. To achieve this, I first connected the voltage measurement block to the generic battery terminals, and then its output signal to the controlled voltage sensor block. However, the voltage I receive this way is the initial voltage value of the generic battery and does not decrease over time depending on the load. To ensure this voltage decrease, I connected a controlled current source as a load to the generic battery. I connected the signal I received from the current sensor on the Simscape side to the input of the current source. But this gave the following algebraic loop error.
Is there a different way to connect the generic battery block to the simscape blocks without getting an algebraic loop error?
Thanks for any help.
Best Regards
Block diagram ‘DronePropellerMotorBattery_v5’ contains 1 algebraic loop(s). To see more details about the loops use the command Simulink.BlockDiagram.getAlgebraicLoops(‘DronePropellerMotorBattery_v5’) or the command line Simulink debugger by typing sldebug(‘DronePropellerMotorBattery_v5’) in the MATLAB command window. To eliminate this message, set Algebraic loop to "none".
Component:Simulink | Category:Block diagram warning
Found algebraic loop containing: DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/INPUT_1_1_1DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/STATE_1 (discontinuity) DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/OUTPUT_1_0DronePropellerMotorBattery_v5/powergui/EquivalentModel1/State-SpaceDronePropellerMotorBattery_v5/Voltage Measurement/do not delete this gain (algebraic variable)
Component:Simulink | Category:Model
Discontinuities detected within algebraic loop(s), may have trouble solving
Component:Simulink | Category:Model warning
Convergence problem when solving algebraic loop containing ‘DronePropellerMotorBattery_v5/Solver Configuration’ at time 0.0. Simulink will try to solve this loop using Simulink 3 (R11) strategy. Use feature(‘ModeIterationsInAlgLoops’,0) to disable the strategy introduced in Simulink 4 (R12)
Component:Simulink | Category:Block warning
Solver is encountering difficulty in simulating model ‘DronePropellerMotorBattery_v5’ at time 1.1125369292536007e-308. Simulink will continue to simulate with warnings. Please check the model for errors.
Caused by:
Solver was unable to reduce the step size without violating the minimum step size value of 2.22507E-308 for 1 consecutive times at time 0. This usually may be caused by violating algebraic constraints in the differential-algebraic system or by the high stiffness of the system. Try tightening the error tolerances, and/or the tolerances for computing consistent conditions. If the problem persists, please check the system or increase the solver Number of consecutive min steps violation parameter.
Suggested Actions
Open Solver Profiler to explore such issuesOpen
Component:Simulink | Category:Block diagram warning
An error occurred during simulation and the simulation was terminated
Caused by:
Algebraic state in algebraic loop containing ‘DronePropellerMotorBattery_v5/Voltage Measurement/do not delete this gain’ computed at time 2.2250738585072014E-308 is Inf or NaN. There may be a singularity in the solution. If the model is correct, try reducing the step size (either by reducing the fixed step size or by tightening the error tolerances) or tweak the initial guess of algebraic loop variable values. Hello,
I’m new to Matlab-Simulink. I cannot directly connect the generic battery block to the simelectronics blocks in the model. To achieve this, I first connected the voltage measurement block to the generic battery terminals, and then its output signal to the controlled voltage sensor block. However, the voltage I receive this way is the initial voltage value of the generic battery and does not decrease over time depending on the load. To ensure this voltage decrease, I connected a controlled current source as a load to the generic battery. I connected the signal I received from the current sensor on the Simscape side to the input of the current source. But this gave the following algebraic loop error.
Is there a different way to connect the generic battery block to the simscape blocks without getting an algebraic loop error?
Thanks for any help.
Best Regards
Block diagram ‘DronePropellerMotorBattery_v5’ contains 1 algebraic loop(s). To see more details about the loops use the command Simulink.BlockDiagram.getAlgebraicLoops(‘DronePropellerMotorBattery_v5’) or the command line Simulink debugger by typing sldebug(‘DronePropellerMotorBattery_v5’) in the MATLAB command window. To eliminate this message, set Algebraic loop to "none".
Component:Simulink | Category:Block diagram warning
Found algebraic loop containing: DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/INPUT_1_1_1DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/STATE_1 (discontinuity) DronePropellerMotorBattery_v5/Solver Configuration/EVAL_KEY/OUTPUT_1_0DronePropellerMotorBattery_v5/powergui/EquivalentModel1/State-SpaceDronePropellerMotorBattery_v5/Voltage Measurement/do not delete this gain (algebraic variable)
Component:Simulink | Category:Model
Discontinuities detected within algebraic loop(s), may have trouble solving
Component:Simulink | Category:Model warning
Convergence problem when solving algebraic loop containing ‘DronePropellerMotorBattery_v5/Solver Configuration’ at time 0.0. Simulink will try to solve this loop using Simulink 3 (R11) strategy. Use feature(‘ModeIterationsInAlgLoops’,0) to disable the strategy introduced in Simulink 4 (R12)
Component:Simulink | Category:Block warning
Solver is encountering difficulty in simulating model ‘DronePropellerMotorBattery_v5’ at time 1.1125369292536007e-308. Simulink will continue to simulate with warnings. Please check the model for errors.
Caused by:
Solver was unable to reduce the step size without violating the minimum step size value of 2.22507E-308 for 1 consecutive times at time 0. This usually may be caused by violating algebraic constraints in the differential-algebraic system or by the high stiffness of the system. Try tightening the error tolerances, and/or the tolerances for computing consistent conditions. If the problem persists, please check the system or increase the solver Number of consecutive min steps violation parameter.
Suggested Actions
Open Solver Profiler to explore such issuesOpen
Component:Simulink | Category:Block diagram warning
An error occurred during simulation and the simulation was terminated
Caused by:
Algebraic state in algebraic loop containing ‘DronePropellerMotorBattery_v5/Voltage Measurement/do not delete this gain’ computed at time 2.2250738585072014E-308 is Inf or NaN. There may be a singularity in the solution. If the model is correct, try reducing the step size (either by reducing the fixed step size or by tightening the error tolerances) or tweak the initial guess of algebraic loop variable values. simulink, model, simscape, simpowersystems MATLAB Answers — New Questions
how to model a varying load in Simulink using specialized power system library in single phase
Hello, I’m modelling a Full-Active Hybrid Energy Storage System in Simulink using specialized power system library, For getting proper results, i need a varying load. The varying load prsent in this library is in 3-phase.. I need a single phase varying load. Can anyone help?Hello, I’m modelling a Full-Active Hybrid Energy Storage System in Simulink using specialized power system library, For getting proper results, i need a varying load. The varying load prsent in this library is in 3-phase.. I need a single phase varying load. Can anyone help? Hello, I’m modelling a Full-Active Hybrid Energy Storage System in Simulink using specialized power system library, For getting proper results, i need a varying load. The varying load prsent in this library is in 3-phase.. I need a single phase varying load. Can anyone help? simulink, simpowersystems, load MATLAB Answers — New Questions
How can I solve this non-linear system?
EquationOne = @(x1,x2) ((x1.*6.6912e+5+x2.*1.879456e+6+x1.^2.*x2.^2.*1.204072e+6+x1.^2.*x2.^3.*9.70736e+5-x1.^3.*x2.^2.*2.95528e+5+x1.^2.*x2.^4.*3.58764e+5+x1.^3.*x2.^3.*7.1352e+4-x1.^4.*x2.^2.*2.8628e+5+x1.^2.*x2.^5.*6.6096e+4+x1.^3.*x2.^4.*5.3724e+4-x1.^4.*x2.^3.*4.5354e+4-x1.^5.*x2.^2.*6.447e+4+x1.^2.*x2.^6.*5.643e+3+x1.^3.*x2.^5.*9.252e+3-x1.^4.*x2.^4.*9.0-x1.^5.*x2.^3.*9.558e+3-x1.^6.*x2.^2.*5.571e+3+x1.^2.*x2.^7.*1.62e+2+x1.^3.*x2.^6.*4.86e+2+x1.^4.*x2.^5.*3.24e+2-x1.^5.*x2.^4.*3.24e+2-x1.^6.*x2.^3.*4.86e+2-x1.^7.*x2.^2.*1.62e+2+x1.*x2.*2.355808e+6+x1.*x2.^2.*3.214896e+6+x1.^2.*x2.*4.26544e+5+x1.*x2.^3.*2.192712e+6-x1.^3.*x2.*9.14456e+5+x1.*x2.^4.*8.23048e+5-x1.^4.*x2.*5.97748e+5+x1.*x2.^5.*1.722e+5-x1.^5.*x2.*1.39216e+5+x1.*x2.^6.*1.8666e+4-x1.^6.*x2.*1.1316e+4+x1.*x2.^7.*8.1e+2+x1.^7.*x2.*7.2e+1-x1.^2.*1.29728e+5-x1.^3.*6.41072e+5+x2.^2.*2.563296e+6-x1.^4.*3.81992e+5+x2.^3.*1.860112e+6-x1.^5.*7.7072e+4+x2.^4.*7.89904e+5+x1.^6.*5.3e+3+x2.^5.*2.01492e+5+x1.^7.*4.224e+3+x2.^6.*3.0168e+4+x1.^8.*4.32e+2+x2.^7.*2.43e+3+x2.^8.*8.1e+1+5.64928e+5).*8.1e+1)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*8.8672e+4+x2.*1.52096e+5+x1.^2.*x2.^2.*1.31664e+5+x1.^2.*x2.^3.*6.6072e+4+x1.^3.*x2.^2.*1.576e+4+x1.^2.*x2.^4.*1.5606e+4+x1.^3.*x2.^3.*1.1916e+4-x1.^4.*x2.^2.*6.816e+3+x1.^2.*x2.^5.*1.629e+3+x1.^3.*x2.^4.*2.52e+3+x1.^4.*x2.^3.*1.17e+2-x1.^5.*x2.^2.*2.106e+3+x1.^2.*x2.^6.*5.4e+1+x1.^3.*x2.^5.*1.62e+2+x1.^4.*x2.^4.*1.08e+2-x1.^5.*x2.^3.*1.08e+2-x1.^6.*x2.^2.*1.62e+2+x1.*x2.*2.32736e+5+x1.*x2.^2.*2.50536e+5+x1.^2.*x2.*1.15128e+5+x1.*x2.^3.*1.32696e+5-x1.^3.*x2.*1.1464e+4+x1.*x2.^4.*3.6756e+4-x1.^4.*x2.*2.716e+4+x1.*x2.^5.*5.058e+3-x1.^5.*x2.*8.958e+3+x1.*x2.^6.*2.7e+2-x1.^6.*x2.*1.197e+3-x1.^7.*x2.*5.4e+1+x1.^2.*4.256e+4-x1.^3.*1.5032e+4+x2.^2.*1.78144e+5-x1.^4.*2.1812e+4+x2.^3.*1.07928e+5-x1.^5.*8.596e+3+x2.^4.*3.666e+4-x1.^6.*1.554e+3+x2.^5.*6.96e+3-x1.^7.*1.08e+2+x2.^6.*6.84e+2+x2.^7.*2.7e+1+5.3152e+4).*8.1e+2)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*2.5536e+4+x2.*2.1872e+4+x1.^2.*x2.^2.*1.0552e+4+x1.^2.*x2.^3.*1.83e+3+x1.^3.*x2.^2.*2.034e+3+x1.^2.*x2.^4.*9.0e+1+x1.^3.*x2.^3.*2.16e+2+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*3.7368e+4+x1.*x2.^2.*2.0396e+4+x1.^2.*x2.*2.3436e+4+x1.*x2.^3.*4.872e+3+x1.^3.*x2.*5.12e+3+x1.*x2.^4.*4.23e+2+x1.^4.*x2.*4.5e+1-x1.^5.*x2.*7.2e+1+x1.^2.*1.7564e+4+x1.^3.*2.996e+3+x2.^2.*1.51e+4-x1.^4.*1.234e+3+x2.^3.*5.028e+3-x1.^5.*5.37e+2+x2.^4.*7.86e+2-x1.^6.*5.4e+1+x2.^5.*4.5e+1+1.2688e+4).*4.32e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*4.9136e+4+x2.*5.1168e+4+x1.^2.*x2.^2.*4.76e+4+x1.^2.*x2.^3.*1.4742e+4+x1.^3.*x2.^2.*1.2622e+4+x1.^2.*x2.^4.*1.929e+3+x1.^3.*x2.^3.*3.102e+3+x1.^4.*x2.^2.*1.05e+3+x1.^2.*x2.^5.*7.2e+1+x1.^3.*x2.^4.*2.34e+2+x1.^4.*x2.^3.*2.16e+2-x1.^5.*x2.^2.*3.6e+1+x1.*x2.*9.4896e+4+x1.*x2.^2.*7.5476e+4+x1.^2.*x2.*6.922e+4+x1.*x2.^3.*2.8916e+4+x1.^3.*x2.*1.9388e+4+x1.*x2.^4.*5.313e+3-x1.^4.*x2.*3.79e+2+x1.*x2.^5.*3.69e+2-x1.^5.*x2.*1.119e+3-x1.^6.*x2.*1.44e+2+x1.^2.*4.156e+4+x1.^3.*1.2268e+4+x2.^2.*4.8376e+4-x1.^4.*2.69e+3+x2.^3.*2.3164e+4-x1.^5.*2.719e+3+x2.^4.*5.89e+3-x1.^6.*6.63e+2+x2.^5.*7.41e+2-x1.^7.*5.4e+1+x2.^6.*3.6e+1+2.24e+4).*2.16e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*4.876e+3+x2.*9.092e+3+x1.^2.*x2.^2.*1.8e+2+x1.*x2.*5.204e+3+x1.*x2.^2.*1.257e+3+x1.^2.*x2.*1.221e+3+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*6.38e+2-x1.^3.*2.85e+2+x2.^2.*3.226e+3-x1.^4.*5.4e+1+x2.^3.*4.23e+2+x2.^4.*1.8e+1+7.9e+3).*3.84e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*1.3388e+4+x2.*1.298e+4+x1.^2.*x2.^2.*1.869e+3+x1.^2.*x2.^3.*1.44e+2+x1.^3.*x2.^2.*1.8e+2+x1.*x2.*1.6168e+4+x1.*x2.^2.*6.067e+3+x1.^2.*x2.*6.353e+3+x1.*x2.^3.*8.19e+2+x1.^3.*x2.*8.25e+2+x1.*x2.^4.*1.8e+1+x1.^2.*5.082e+3-x1.^3.*1.15e+2+x2.^2.*6.43e+3-x1.^4.*4.11e+2+x2.^3.*1.359e+3-x1.^5.*5.4e+1+x2.^4.*9.0e+1+9.448e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).*(x1.*3.1424e+4+x2.*2.016e+4+x1.^2.*x2.^2.*2.1836e+4+x1.^2.*x2.^3.*5.256e+3+x1.^3.*x2.^2.*7.756e+3+x1.^2.*x2.^4.*5.7e+2+x1.^3.*x2.^3.*1.266e+3+x1.^4.*x2.^2.*1.257e+3+x1.^2.*x2.^5.*1.8e+1+x1.^3.*x2.^4.*7.2e+1+x1.^4.*x2.^3.*1.08e+2+x1.^5.*x2.^2.*7.2e+1+x1.*x2.*4.608e+4+x1.*x2.^2.*2.9192e+4+x1.^2.*x2.*4.4224e+4+x1.*x2.^3.*9.452e+3+x1.^3.*x2.*2.1188e+4+x1.*x2.^4.*1.542e+3+x1.^4.*x2.*5.096e+3+x1.*x2.^5.*9.9e+1+x1.^5.*x2.*5.55e+2+x1.^6.*x2.*1.8e+1+x1.^2.*3.748e+4+x1.^3.*2.2984e+4+x2.^2.*1.6168e+4+x1.^4.*7.528e+3+x2.^3.*6.88e+3+x1.^5.*1.262e+3+x2.^4.*1.612e+3+x1.^6.*8.4e+1+x2.^5.*1.92e+2+x2.^6.*9.0+1.1072e+4).*3.24e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).^2.*(x1.*5.0e+3+x2.*3.792e+3+x1.^2.*x2.^2.*5.34e+2+x1.^2.*x2.^3.*3.6e+1+x1.^3.*x2.^2.*5.4e+1+x1.*x2.*5.192e+3+x1.*x2.^2.*1.766e+3+x1.^2.*x2.*2.246e+3+x1.*x2.^3.*2.34e+2+x1.^3.*x2.*4.38e+2+x1.*x2.^4.*9.0+x1.^4.*x2.*3.6e+1+x1.^2.*2.628e+3+x1.^3.*7.14e+2+x2.^2.*1.652e+3+x1.^4.*1.2e+2+x2.^3.*3.06e+2+x1.^5.*9.0+x2.^4.*1.8e+1+3.152e+3).*1.296e+4)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*1.8224e+4+x2.*9.824e+3+x1.^2.*x2.^2.*4.916e+3+x1.^2.*x2.^3.*5.82e+2+x1.^3.*x2.^2.*1.218e+3+x1.^2.*x2.^4.*1.8e+1+x1.^3.*x2.^3.*7.2e+1+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*1.9984e+4+x1.*x2.^2.*8.332e+3+x1.^2.*x2.*1.5908e+4+x1.*x2.^3.*1.512e+3+x1.^3.*x2.*5.88e+3+x1.*x2.^4.*9.9e+1+x1.^4.*x2.*1.041e+3+x1.^5.*x2.*7.2e+1+x1.^2.*1.7896e+4+x1.^3.*8.66e+3+x2.^2.*5.432e+3+x1.^4.*2.274e+3+x2.^3.*1.452e+3+x1.^5.*3.15e+2+x2.^4.*1.86e+2+x1.^6.*1.8e+1+x2.^5.*9.0+7.296e+3).*6.48e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*4.264e+3+x2.*6.752e+3+x1.^2.*x2.^2.*1.08e+2+x1.*x2.*3.944e+3+x1.*x2.^2.*9.33e+2+x1.^2.*x2.*8.61e+2+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*1.07e+3+x1.^3.*1.83e+2+x2.^2.*2.362e+3+x1.^4.*1.8e+1+x2.^3.*3.51e+2+x2.^4.*1.8e+1+6.352e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).^2.*(x1.*1.688e+3+x2.*2.488e+3+x1.^2.*x2.^2.*5.4e+1+x1.*x2.*1.6e+3+x1.*x2.^2.*4.26e+2+x1.^2.*x2.*3.66e+2+x1.*x2.^3.*3.6e+1+x1.^3.*x2.*3.6e+1+x1.^2.*4.96e+2+x1.^3.*1.02e+2+x2.^2.*9.68e+2+x1.^4.*9.0+x2.^3.*1.62e+2+x2.^4.*9.0+2.216e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*1.1096e+4+x2.*8.888e+3+x1.^2.*x2.^2.*1.122e+3+x1.^2.*x2.^3.*7.2e+1+x1.^3.*x2.^2.*1.08e+2+x1.*x2.*1.1488e+4+x1.*x2.^2.*3.904e+3+x1.^2.*x2.*4.772e+3+x1.*x2.^3.*5.04e+2+x1.^3.*x2.*9.12e+2+x1.*x2.^4.*1.8e+1+x1.^4.*x2.*7.2e+1+x1.^2.*5.7e+3+x1.^3.*1.544e+3+x2.^2.*3.988e+3+x1.^4.*2.49e+2+x2.^3.*7.56e+2+x1.^5.*1.8e+1+x2.^4.*4.5e+1+7.24e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8);
EquationTwo = @(x1,x2) ((x1.*6.6912e+5+x2.*1.879456e+6+x1.^2.*x2.^2.*1.204072e+6+x1.^2.*x2.^3.*9.70736e+5-x1.^3.*x2.^2.*2.95528e+5+x1.^2.*x2.^4.*3.58764e+5+x1.^3.*x2.^3.*7.1352e+4-x1.^4.*x2.^2.*2.8628e+5+x1.^2.*x2.^5.*6.6096e+4+x1.^3.*x2.^4.*5.3724e+4-x1.^4.*x2.^3.*4.5354e+4-x1.^5.*x2.^2.*6.447e+4+x1.^2.*x2.^6.*5.643e+3+x1.^3.*x2.^5.*9.252e+3-x1.^4.*x2.^4.*9.0-x1.^5.*x2.^3.*9.558e+3-x1.^6.*x2.^2.*5.571e+3+x1.^2.*x2.^7.*1.62e+2+x1.^3.*x2.^6.*4.86e+2+x1.^4.*x2.^5.*3.24e+2-x1.^5.*x2.^4.*3.24e+2-x1.^6.*x2.^3.*4.86e+2-x1.^7.*x2.^2.*1.62e+2+x1.*x2.*2.355808e+6+x1.*x2.^2.*3.214896e+6+x1.^2.*x2.*4.26544e+5+x1.*x2.^3.*2.192712e+6-x1.^3.*x2.*9.14456e+5+x1.*x2.^4.*8.23048e+5-x1.^4.*x2.*5.97748e+5+x1.*x2.^5.*1.722e+5-x1.^5.*x2.*1.39216e+5+x1.*x2.^6.*1.8666e+4-x1.^6.*x2.*1.1316e+4+x1.*x2.^7.*8.1e+2+x1.^7.*x2.*7.2e+1-x1.^2.*1.29728e+5-x1.^3.*6.41072e+5+x2.^2.*2.563296e+6-x1.^4.*3.81992e+5+x2.^3.*1.860112e+6-x1.^5.*7.7072e+4+x2.^4.*7.89904e+5+x1.^6.*5.3e+3+x2.^5.*2.01492e+5+x1.^7.*4.224e+3+x2.^6.*3.0168e+4+x1.^8.*4.32e+2+x2.^7.*2.43e+3+x2.^8.*8.1e+1+5.64928e+5).*8.1e+1)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*2.5536e+4+x2.*2.1872e+4+x1.^2.*x2.^2.*1.0552e+4+x1.^2.*x2.^3.*1.83e+3+x1.^3.*x2.^2.*2.034e+3+x1.^2.*x2.^4.*9.0e+1+x1.^3.*x2.^3.*2.16e+2+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*3.7368e+4+x1.*x2.^2.*2.0396e+4+x1.^2.*x2.*2.3436e+4+x1.*x2.^3.*4.872e+3+x1.^3.*x2.*5.12e+3+x1.*x2.^4.*4.23e+2+x1.^4.*x2.*4.5e+1-x1.^5.*x2.*7.2e+1+x1.^2.*1.7564e+4+x1.^3.*2.996e+3+x2.^2.*1.51e+4-x1.^4.*1.234e+3+x2.^3.*5.028e+3-x1.^5.*5.37e+2+x2.^4.*7.86e+2-x1.^6.*5.4e+1+x2.^5.*4.5e+1+1.2688e+4).*4.32e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*4.9136e+4+x2.*5.1168e+4+x1.^2.*x2.^2.*4.76e+4+x1.^2.*x2.^3.*1.4742e+4+x1.^3.*x2.^2.*1.2622e+4+x1.^2.*x2.^4.*1.929e+3+x1.^3.*x2.^3.*3.102e+3+x1.^4.*x2.^2.*1.05e+3+x1.^2.*x2.^5.*7.2e+1+x1.^3.*x2.^4.*2.34e+2+x1.^4.*x2.^3.*2.16e+2-x1.^5.*x2.^2.*3.6e+1+x1.*x2.*9.4896e+4+x1.*x2.^2.*7.5476e+4+x1.^2.*x2.*6.922e+4+x1.*x2.^3.*2.8916e+4+x1.^3.*x2.*1.9388e+4+x1.*x2.^4.*5.313e+3-x1.^4.*x2.*3.79e+2+x1.*x2.^5.*3.69e+2-x1.^5.*x2.*1.119e+3-x1.^6.*x2.*1.44e+2+x1.^2.*4.156e+4+x1.^3.*1.2268e+4+x2.^2.*4.8376e+4-x1.^4.*2.69e+3+x2.^3.*2.3164e+4-x1.^5.*2.719e+3+x2.^4.*5.89e+3-x1.^6.*6.63e+2+x2.^5.*7.41e+2-x1.^7.*5.4e+1+x2.^6.*3.6e+1+2.24e+4).*2.16e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*4.876e+3+x2.*9.092e+3+x1.^2.*x2.^2.*1.8e+2+x1.*x2.*5.204e+3+x1.*x2.^2.*1.257e+3+x1.^2.*x2.*1.221e+3+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*6.38e+2-x1.^3.*2.85e+2+x2.^2.*3.226e+3-x1.^4.*5.4e+1+x2.^3.*4.23e+2+x2.^4.*1.8e+1+7.9e+3).*3.84e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*1.3388e+4+x2.*1.298e+4+x1.^2.*x2.^2.*1.869e+3+x1.^2.*x2.^3.*1.44e+2+x1.^3.*x2.^2.*1.8e+2+x1.*x2.*1.6168e+4+x1.*x2.^2.*6.067e+3+x1.^2.*x2.*6.353e+3+x1.*x2.^3.*8.19e+2+x1.^3.*x2.*8.25e+2+x1.*x2.^4.*1.8e+1+x1.^2.*5.082e+3-x1.^3.*1.15e+2+x2.^2.*6.43e+3-x1.^4.*4.11e+2+x2.^3.*1.359e+3-x1.^5.*5.4e+1+x2.^4.*9.0e+1+9.448e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).*(x1.*3.1424e+4+x2.*2.016e+4+x1.^2.*x2.^2.*2.1836e+4+x1.^2.*x2.^3.*5.256e+3+x1.^3.*x2.^2.*7.756e+3+x1.^2.*x2.^4.*5.7e+2+x1.^3.*x2.^3.*1.266e+3+x1.^4.*x2.^2.*1.257e+3+x1.^2.*x2.^5.*1.8e+1+x1.^3.*x2.^4.*7.2e+1+x1.^4.*x2.^3.*1.08e+2+x1.^5.*x2.^2.*7.2e+1+x1.*x2.*4.608e+4+x1.*x2.^2.*2.9192e+4+x1.^2.*x2.*4.4224e+4+x1.*x2.^3.*9.452e+3+x1.^3.*x2.*2.1188e+4+x1.*x2.^4.*1.542e+3+x1.^4.*x2.*5.096e+3+x1.*x2.^5.*9.9e+1+x1.^5.*x2.*5.55e+2+x1.^6.*x2.*1.8e+1+x1.^2.*3.748e+4+x1.^3.*2.2984e+4+x2.^2.*1.6168e+4+x1.^4.*7.528e+3+x2.^3.*6.88e+3+x1.^5.*1.262e+3+x2.^4.*1.612e+3+x1.^6.*8.4e+1+x2.^5.*1.92e+2+x2.^6.*9.0+1.1072e+4).*3.24e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).^2.*(x1.*5.0e+3+x2.*3.792e+3+x1.^2.*x2.^2.*5.34e+2+x1.^2.*x2.^3.*3.6e+1+x1.^3.*x2.^2.*5.4e+1+x1.*x2.*5.192e+3+x1.*x2.^2.*1.766e+3+x1.^2.*x2.*2.246e+3+x1.*x2.^3.*2.34e+2+x1.^3.*x2.*4.38e+2+x1.*x2.^4.*9.0+x1.^4.*x2.*3.6e+1+x1.^2.*2.628e+3+x1.^3.*7.14e+2+x2.^2.*1.652e+3+x1.^4.*1.2e+2+x2.^3.*3.06e+2+x1.^5.*9.0+x2.^4.*1.8e+1+3.152e+3).*1.296e+4)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x2+1.0).*(x1.*4.264e+3+x2.*6.752e+3+x1.^2.*x2.^2.*1.08e+2+x1.*x2.*3.944e+3+x1.*x2.^2.*9.33e+2+x1.^2.*x2.*8.61e+2+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*1.07e+3+x1.^3.*1.83e+2+x2.^2.*2.362e+3+x1.^4.*1.8e+1+x2.^3.*3.51e+2+x2.^4.*1.8e+1+6.352e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x2+1.0).*(x1.*1.1096e+4+x2.*8.888e+3+x1.^2.*x2.^2.*1.122e+3+x1.^2.*x2.^3.*7.2e+1+x1.^3.*x2.^2.*1.08e+2+x1.*x2.*1.1488e+4+x1.*x2.^2.*3.904e+3+x1.^2.*x2.*4.772e+3+x1.*x2.^3.*5.04e+2+x1.^3.*x2.*9.12e+2+x1.*x2.^4.*1.8e+1+x1.^4.*x2.*7.2e+1+x1.^2.*5.7e+3+x1.^3.*1.544e+3+x2.^2.*3.988e+3+x1.^4.*2.49e+2+x2.^3.*7.56e+2+x1.^5.*1.8e+1+x2.^4.*4.5e+1+7.24e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8);EquationOne = @(x1,x2) ((x1.*6.6912e+5+x2.*1.879456e+6+x1.^2.*x2.^2.*1.204072e+6+x1.^2.*x2.^3.*9.70736e+5-x1.^3.*x2.^2.*2.95528e+5+x1.^2.*x2.^4.*3.58764e+5+x1.^3.*x2.^3.*7.1352e+4-x1.^4.*x2.^2.*2.8628e+5+x1.^2.*x2.^5.*6.6096e+4+x1.^3.*x2.^4.*5.3724e+4-x1.^4.*x2.^3.*4.5354e+4-x1.^5.*x2.^2.*6.447e+4+x1.^2.*x2.^6.*5.643e+3+x1.^3.*x2.^5.*9.252e+3-x1.^4.*x2.^4.*9.0-x1.^5.*x2.^3.*9.558e+3-x1.^6.*x2.^2.*5.571e+3+x1.^2.*x2.^7.*1.62e+2+x1.^3.*x2.^6.*4.86e+2+x1.^4.*x2.^5.*3.24e+2-x1.^5.*x2.^4.*3.24e+2-x1.^6.*x2.^3.*4.86e+2-x1.^7.*x2.^2.*1.62e+2+x1.*x2.*2.355808e+6+x1.*x2.^2.*3.214896e+6+x1.^2.*x2.*4.26544e+5+x1.*x2.^3.*2.192712e+6-x1.^3.*x2.*9.14456e+5+x1.*x2.^4.*8.23048e+5-x1.^4.*x2.*5.97748e+5+x1.*x2.^5.*1.722e+5-x1.^5.*x2.*1.39216e+5+x1.*x2.^6.*1.8666e+4-x1.^6.*x2.*1.1316e+4+x1.*x2.^7.*8.1e+2+x1.^7.*x2.*7.2e+1-x1.^2.*1.29728e+5-x1.^3.*6.41072e+5+x2.^2.*2.563296e+6-x1.^4.*3.81992e+5+x2.^3.*1.860112e+6-x1.^5.*7.7072e+4+x2.^4.*7.89904e+5+x1.^6.*5.3e+3+x2.^5.*2.01492e+5+x1.^7.*4.224e+3+x2.^6.*3.0168e+4+x1.^8.*4.32e+2+x2.^7.*2.43e+3+x2.^8.*8.1e+1+5.64928e+5).*8.1e+1)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*8.8672e+4+x2.*1.52096e+5+x1.^2.*x2.^2.*1.31664e+5+x1.^2.*x2.^3.*6.6072e+4+x1.^3.*x2.^2.*1.576e+4+x1.^2.*x2.^4.*1.5606e+4+x1.^3.*x2.^3.*1.1916e+4-x1.^4.*x2.^2.*6.816e+3+x1.^2.*x2.^5.*1.629e+3+x1.^3.*x2.^4.*2.52e+3+x1.^4.*x2.^3.*1.17e+2-x1.^5.*x2.^2.*2.106e+3+x1.^2.*x2.^6.*5.4e+1+x1.^3.*x2.^5.*1.62e+2+x1.^4.*x2.^4.*1.08e+2-x1.^5.*x2.^3.*1.08e+2-x1.^6.*x2.^2.*1.62e+2+x1.*x2.*2.32736e+5+x1.*x2.^2.*2.50536e+5+x1.^2.*x2.*1.15128e+5+x1.*x2.^3.*1.32696e+5-x1.^3.*x2.*1.1464e+4+x1.*x2.^4.*3.6756e+4-x1.^4.*x2.*2.716e+4+x1.*x2.^5.*5.058e+3-x1.^5.*x2.*8.958e+3+x1.*x2.^6.*2.7e+2-x1.^6.*x2.*1.197e+3-x1.^7.*x2.*5.4e+1+x1.^2.*4.256e+4-x1.^3.*1.5032e+4+x2.^2.*1.78144e+5-x1.^4.*2.1812e+4+x2.^3.*1.07928e+5-x1.^5.*8.596e+3+x2.^4.*3.666e+4-x1.^6.*1.554e+3+x2.^5.*6.96e+3-x1.^7.*1.08e+2+x2.^6.*6.84e+2+x2.^7.*2.7e+1+5.3152e+4).*8.1e+2)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*2.5536e+4+x2.*2.1872e+4+x1.^2.*x2.^2.*1.0552e+4+x1.^2.*x2.^3.*1.83e+3+x1.^3.*x2.^2.*2.034e+3+x1.^2.*x2.^4.*9.0e+1+x1.^3.*x2.^3.*2.16e+2+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*3.7368e+4+x1.*x2.^2.*2.0396e+4+x1.^2.*x2.*2.3436e+4+x1.*x2.^3.*4.872e+3+x1.^3.*x2.*5.12e+3+x1.*x2.^4.*4.23e+2+x1.^4.*x2.*4.5e+1-x1.^5.*x2.*7.2e+1+x1.^2.*1.7564e+4+x1.^3.*2.996e+3+x2.^2.*1.51e+4-x1.^4.*1.234e+3+x2.^3.*5.028e+3-x1.^5.*5.37e+2+x2.^4.*7.86e+2-x1.^6.*5.4e+1+x2.^5.*4.5e+1+1.2688e+4).*4.32e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*4.9136e+4+x2.*5.1168e+4+x1.^2.*x2.^2.*4.76e+4+x1.^2.*x2.^3.*1.4742e+4+x1.^3.*x2.^2.*1.2622e+4+x1.^2.*x2.^4.*1.929e+3+x1.^3.*x2.^3.*3.102e+3+x1.^4.*x2.^2.*1.05e+3+x1.^2.*x2.^5.*7.2e+1+x1.^3.*x2.^4.*2.34e+2+x1.^4.*x2.^3.*2.16e+2-x1.^5.*x2.^2.*3.6e+1+x1.*x2.*9.4896e+4+x1.*x2.^2.*7.5476e+4+x1.^2.*x2.*6.922e+4+x1.*x2.^3.*2.8916e+4+x1.^3.*x2.*1.9388e+4+x1.*x2.^4.*5.313e+3-x1.^4.*x2.*3.79e+2+x1.*x2.^5.*3.69e+2-x1.^5.*x2.*1.119e+3-x1.^6.*x2.*1.44e+2+x1.^2.*4.156e+4+x1.^3.*1.2268e+4+x2.^2.*4.8376e+4-x1.^4.*2.69e+3+x2.^3.*2.3164e+4-x1.^5.*2.719e+3+x2.^4.*5.89e+3-x1.^6.*6.63e+2+x2.^5.*7.41e+2-x1.^7.*5.4e+1+x2.^6.*3.6e+1+2.24e+4).*2.16e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*4.876e+3+x2.*9.092e+3+x1.^2.*x2.^2.*1.8e+2+x1.*x2.*5.204e+3+x1.*x2.^2.*1.257e+3+x1.^2.*x2.*1.221e+3+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*6.38e+2-x1.^3.*2.85e+2+x2.^2.*3.226e+3-x1.^4.*5.4e+1+x2.^3.*4.23e+2+x2.^4.*1.8e+1+7.9e+3).*3.84e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*1.3388e+4+x2.*1.298e+4+x1.^2.*x2.^2.*1.869e+3+x1.^2.*x2.^3.*1.44e+2+x1.^3.*x2.^2.*1.8e+2+x1.*x2.*1.6168e+4+x1.*x2.^2.*6.067e+3+x1.^2.*x2.*6.353e+3+x1.*x2.^3.*8.19e+2+x1.^3.*x2.*8.25e+2+x1.*x2.^4.*1.8e+1+x1.^2.*5.082e+3-x1.^3.*1.15e+2+x2.^2.*6.43e+3-x1.^4.*4.11e+2+x2.^3.*1.359e+3-x1.^5.*5.4e+1+x2.^4.*9.0e+1+9.448e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).*(x1.*3.1424e+4+x2.*2.016e+4+x1.^2.*x2.^2.*2.1836e+4+x1.^2.*x2.^3.*5.256e+3+x1.^3.*x2.^2.*7.756e+3+x1.^2.*x2.^4.*5.7e+2+x1.^3.*x2.^3.*1.266e+3+x1.^4.*x2.^2.*1.257e+3+x1.^2.*x2.^5.*1.8e+1+x1.^3.*x2.^4.*7.2e+1+x1.^4.*x2.^3.*1.08e+2+x1.^5.*x2.^2.*7.2e+1+x1.*x2.*4.608e+4+x1.*x2.^2.*2.9192e+4+x1.^2.*x2.*4.4224e+4+x1.*x2.^3.*9.452e+3+x1.^3.*x2.*2.1188e+4+x1.*x2.^4.*1.542e+3+x1.^4.*x2.*5.096e+3+x1.*x2.^5.*9.9e+1+x1.^5.*x2.*5.55e+2+x1.^6.*x2.*1.8e+1+x1.^2.*3.748e+4+x1.^3.*2.2984e+4+x2.^2.*1.6168e+4+x1.^4.*7.528e+3+x2.^3.*6.88e+3+x1.^5.*1.262e+3+x2.^4.*1.612e+3+x1.^6.*8.4e+1+x2.^5.*1.92e+2+x2.^6.*9.0+1.1072e+4).*3.24e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).^2.*(x1.*5.0e+3+x2.*3.792e+3+x1.^2.*x2.^2.*5.34e+2+x1.^2.*x2.^3.*3.6e+1+x1.^3.*x2.^2.*5.4e+1+x1.*x2.*5.192e+3+x1.*x2.^2.*1.766e+3+x1.^2.*x2.*2.246e+3+x1.*x2.^3.*2.34e+2+x1.^3.*x2.*4.38e+2+x1.*x2.^4.*9.0+x1.^4.*x2.*3.6e+1+x1.^2.*2.628e+3+x1.^3.*7.14e+2+x2.^2.*1.652e+3+x1.^4.*1.2e+2+x2.^3.*3.06e+2+x1.^5.*9.0+x2.^4.*1.8e+1+3.152e+3).*1.296e+4)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*1.8224e+4+x2.*9.824e+3+x1.^2.*x2.^2.*4.916e+3+x1.^2.*x2.^3.*5.82e+2+x1.^3.*x2.^2.*1.218e+3+x1.^2.*x2.^4.*1.8e+1+x1.^3.*x2.^3.*7.2e+1+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*1.9984e+4+x1.*x2.^2.*8.332e+3+x1.^2.*x2.*1.5908e+4+x1.*x2.^3.*1.512e+3+x1.^3.*x2.*5.88e+3+x1.*x2.^4.*9.9e+1+x1.^4.*x2.*1.041e+3+x1.^5.*x2.*7.2e+1+x1.^2.*1.7896e+4+x1.^3.*8.66e+3+x2.^2.*5.432e+3+x1.^4.*2.274e+3+x2.^3.*1.452e+3+x1.^5.*3.15e+2+x2.^4.*1.86e+2+x1.^6.*1.8e+1+x2.^5.*9.0+7.296e+3).*6.48e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*4.264e+3+x2.*6.752e+3+x1.^2.*x2.^2.*1.08e+2+x1.*x2.*3.944e+3+x1.*x2.^2.*9.33e+2+x1.^2.*x2.*8.61e+2+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*1.07e+3+x1.^3.*1.83e+2+x2.^2.*2.362e+3+x1.^4.*1.8e+1+x2.^3.*3.51e+2+x2.^4.*1.8e+1+6.352e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).^2.*(x1.*1.688e+3+x2.*2.488e+3+x1.^2.*x2.^2.*5.4e+1+x1.*x2.*1.6e+3+x1.*x2.^2.*4.26e+2+x1.^2.*x2.*3.66e+2+x1.*x2.^3.*3.6e+1+x1.^3.*x2.*3.6e+1+x1.^2.*4.96e+2+x1.^3.*1.02e+2+x2.^2.*9.68e+2+x1.^4.*9.0+x2.^3.*1.62e+2+x2.^4.*9.0+2.216e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*1.1096e+4+x2.*8.888e+3+x1.^2.*x2.^2.*1.122e+3+x1.^2.*x2.^3.*7.2e+1+x1.^3.*x2.^2.*1.08e+2+x1.*x2.*1.1488e+4+x1.*x2.^2.*3.904e+3+x1.^2.*x2.*4.772e+3+x1.*x2.^3.*5.04e+2+x1.^3.*x2.*9.12e+2+x1.*x2.^4.*1.8e+1+x1.^4.*x2.*7.2e+1+x1.^2.*5.7e+3+x1.^3.*1.544e+3+x2.^2.*3.988e+3+x1.^4.*2.49e+2+x2.^3.*7.56e+2+x1.^5.*1.8e+1+x2.^4.*4.5e+1+7.24e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8);
EquationTwo = @(x1,x2) ((x1.*6.6912e+5+x2.*1.879456e+6+x1.^2.*x2.^2.*1.204072e+6+x1.^2.*x2.^3.*9.70736e+5-x1.^3.*x2.^2.*2.95528e+5+x1.^2.*x2.^4.*3.58764e+5+x1.^3.*x2.^3.*7.1352e+4-x1.^4.*x2.^2.*2.8628e+5+x1.^2.*x2.^5.*6.6096e+4+x1.^3.*x2.^4.*5.3724e+4-x1.^4.*x2.^3.*4.5354e+4-x1.^5.*x2.^2.*6.447e+4+x1.^2.*x2.^6.*5.643e+3+x1.^3.*x2.^5.*9.252e+3-x1.^4.*x2.^4.*9.0-x1.^5.*x2.^3.*9.558e+3-x1.^6.*x2.^2.*5.571e+3+x1.^2.*x2.^7.*1.62e+2+x1.^3.*x2.^6.*4.86e+2+x1.^4.*x2.^5.*3.24e+2-x1.^5.*x2.^4.*3.24e+2-x1.^6.*x2.^3.*4.86e+2-x1.^7.*x2.^2.*1.62e+2+x1.*x2.*2.355808e+6+x1.*x2.^2.*3.214896e+6+x1.^2.*x2.*4.26544e+5+x1.*x2.^3.*2.192712e+6-x1.^3.*x2.*9.14456e+5+x1.*x2.^4.*8.23048e+5-x1.^4.*x2.*5.97748e+5+x1.*x2.^5.*1.722e+5-x1.^5.*x2.*1.39216e+5+x1.*x2.^6.*1.8666e+4-x1.^6.*x2.*1.1316e+4+x1.*x2.^7.*8.1e+2+x1.^7.*x2.*7.2e+1-x1.^2.*1.29728e+5-x1.^3.*6.41072e+5+x2.^2.*2.563296e+6-x1.^4.*3.81992e+5+x2.^3.*1.860112e+6-x1.^5.*7.7072e+4+x2.^4.*7.89904e+5+x1.^6.*5.3e+3+x2.^5.*2.01492e+5+x1.^7.*4.224e+3+x2.^6.*3.0168e+4+x1.^8.*4.32e+2+x2.^7.*2.43e+3+x2.^8.*8.1e+1+5.64928e+5).*8.1e+1)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*2.5536e+4+x2.*2.1872e+4+x1.^2.*x2.^2.*1.0552e+4+x1.^2.*x2.^3.*1.83e+3+x1.^3.*x2.^2.*2.034e+3+x1.^2.*x2.^4.*9.0e+1+x1.^3.*x2.^3.*2.16e+2+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*3.7368e+4+x1.*x2.^2.*2.0396e+4+x1.^2.*x2.*2.3436e+4+x1.*x2.^3.*4.872e+3+x1.^3.*x2.*5.12e+3+x1.*x2.^4.*4.23e+2+x1.^4.*x2.*4.5e+1-x1.^5.*x2.*7.2e+1+x1.^2.*1.7564e+4+x1.^3.*2.996e+3+x2.^2.*1.51e+4-x1.^4.*1.234e+3+x2.^3.*5.028e+3-x1.^5.*5.37e+2+x2.^4.*7.86e+2-x1.^6.*5.4e+1+x2.^5.*4.5e+1+1.2688e+4).*4.32e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*4.9136e+4+x2.*5.1168e+4+x1.^2.*x2.^2.*4.76e+4+x1.^2.*x2.^3.*1.4742e+4+x1.^3.*x2.^2.*1.2622e+4+x1.^2.*x2.^4.*1.929e+3+x1.^3.*x2.^3.*3.102e+3+x1.^4.*x2.^2.*1.05e+3+x1.^2.*x2.^5.*7.2e+1+x1.^3.*x2.^4.*2.34e+2+x1.^4.*x2.^3.*2.16e+2-x1.^5.*x2.^2.*3.6e+1+x1.*x2.*9.4896e+4+x1.*x2.^2.*7.5476e+4+x1.^2.*x2.*6.922e+4+x1.*x2.^3.*2.8916e+4+x1.^3.*x2.*1.9388e+4+x1.*x2.^4.*5.313e+3-x1.^4.*x2.*3.79e+2+x1.*x2.^5.*3.69e+2-x1.^5.*x2.*1.119e+3-x1.^6.*x2.*1.44e+2+x1.^2.*4.156e+4+x1.^3.*1.2268e+4+x2.^2.*4.8376e+4-x1.^4.*2.69e+3+x2.^3.*2.3164e+4-x1.^5.*2.719e+3+x2.^4.*5.89e+3-x1.^6.*6.63e+2+x2.^5.*7.41e+2-x1.^7.*5.4e+1+x2.^6.*3.6e+1+2.24e+4).*2.16e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*4.876e+3+x2.*9.092e+3+x1.^2.*x2.^2.*1.8e+2+x1.*x2.*5.204e+3+x1.*x2.^2.*1.257e+3+x1.^2.*x2.*1.221e+3+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*6.38e+2-x1.^3.*2.85e+2+x2.^2.*3.226e+3-x1.^4.*5.4e+1+x2.^3.*4.23e+2+x2.^4.*1.8e+1+7.9e+3).*3.84e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*1.3388e+4+x2.*1.298e+4+x1.^2.*x2.^2.*1.869e+3+x1.^2.*x2.^3.*1.44e+2+x1.^3.*x2.^2.*1.8e+2+x1.*x2.*1.6168e+4+x1.*x2.^2.*6.067e+3+x1.^2.*x2.*6.353e+3+x1.*x2.^3.*8.19e+2+x1.^3.*x2.*8.25e+2+x1.*x2.^4.*1.8e+1+x1.^2.*5.082e+3-x1.^3.*1.15e+2+x2.^2.*6.43e+3-x1.^4.*4.11e+2+x2.^3.*1.359e+3-x1.^5.*5.4e+1+x2.^4.*9.0e+1+9.448e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).*(x1.*3.1424e+4+x2.*2.016e+4+x1.^2.*x2.^2.*2.1836e+4+x1.^2.*x2.^3.*5.256e+3+x1.^3.*x2.^2.*7.756e+3+x1.^2.*x2.^4.*5.7e+2+x1.^3.*x2.^3.*1.266e+3+x1.^4.*x2.^2.*1.257e+3+x1.^2.*x2.^5.*1.8e+1+x1.^3.*x2.^4.*7.2e+1+x1.^4.*x2.^3.*1.08e+2+x1.^5.*x2.^2.*7.2e+1+x1.*x2.*4.608e+4+x1.*x2.^2.*2.9192e+4+x1.^2.*x2.*4.4224e+4+x1.*x2.^3.*9.452e+3+x1.^3.*x2.*2.1188e+4+x1.*x2.^4.*1.542e+3+x1.^4.*x2.*5.096e+3+x1.*x2.^5.*9.9e+1+x1.^5.*x2.*5.55e+2+x1.^6.*x2.*1.8e+1+x1.^2.*3.748e+4+x1.^3.*2.2984e+4+x2.^2.*1.6168e+4+x1.^4.*7.528e+3+x2.^3.*6.88e+3+x1.^5.*1.262e+3+x2.^4.*1.612e+3+x1.^6.*8.4e+1+x2.^5.*1.92e+2+x2.^6.*9.0+1.1072e+4).*3.24e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).^2.*(x1.*5.0e+3+x2.*3.792e+3+x1.^2.*x2.^2.*5.34e+2+x1.^2.*x2.^3.*3.6e+1+x1.^3.*x2.^2.*5.4e+1+x1.*x2.*5.192e+3+x1.*x2.^2.*1.766e+3+x1.^2.*x2.*2.246e+3+x1.*x2.^3.*2.34e+2+x1.^3.*x2.*4.38e+2+x1.*x2.^4.*9.0+x1.^4.*x2.*3.6e+1+x1.^2.*2.628e+3+x1.^3.*7.14e+2+x2.^2.*1.652e+3+x1.^4.*1.2e+2+x2.^3.*3.06e+2+x1.^5.*9.0+x2.^4.*1.8e+1+3.152e+3).*1.296e+4)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x2+1.0).*(x1.*4.264e+3+x2.*6.752e+3+x1.^2.*x2.^2.*1.08e+2+x1.*x2.*3.944e+3+x1.*x2.^2.*9.33e+2+x1.^2.*x2.*8.61e+2+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*1.07e+3+x1.^3.*1.83e+2+x2.^2.*2.362e+3+x1.^4.*1.8e+1+x2.^3.*3.51e+2+x2.^4.*1.8e+1+6.352e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x2+1.0).*(x1.*1.1096e+4+x2.*8.888e+3+x1.^2.*x2.^2.*1.122e+3+x1.^2.*x2.^3.*7.2e+1+x1.^3.*x2.^2.*1.08e+2+x1.*x2.*1.1488e+4+x1.*x2.^2.*3.904e+3+x1.^2.*x2.*4.772e+3+x1.*x2.^3.*5.04e+2+x1.^3.*x2.*9.12e+2+x1.*x2.^4.*1.8e+1+x1.^4.*x2.*7.2e+1+x1.^2.*5.7e+3+x1.^3.*1.544e+3+x2.^2.*3.988e+3+x1.^4.*2.49e+2+x2.^3.*7.56e+2+x1.^5.*1.8e+1+x2.^4.*4.5e+1+7.24e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8); EquationOne = @(x1,x2) ((x1.*6.6912e+5+x2.*1.879456e+6+x1.^2.*x2.^2.*1.204072e+6+x1.^2.*x2.^3.*9.70736e+5-x1.^3.*x2.^2.*2.95528e+5+x1.^2.*x2.^4.*3.58764e+5+x1.^3.*x2.^3.*7.1352e+4-x1.^4.*x2.^2.*2.8628e+5+x1.^2.*x2.^5.*6.6096e+4+x1.^3.*x2.^4.*5.3724e+4-x1.^4.*x2.^3.*4.5354e+4-x1.^5.*x2.^2.*6.447e+4+x1.^2.*x2.^6.*5.643e+3+x1.^3.*x2.^5.*9.252e+3-x1.^4.*x2.^4.*9.0-x1.^5.*x2.^3.*9.558e+3-x1.^6.*x2.^2.*5.571e+3+x1.^2.*x2.^7.*1.62e+2+x1.^3.*x2.^6.*4.86e+2+x1.^4.*x2.^5.*3.24e+2-x1.^5.*x2.^4.*3.24e+2-x1.^6.*x2.^3.*4.86e+2-x1.^7.*x2.^2.*1.62e+2+x1.*x2.*2.355808e+6+x1.*x2.^2.*3.214896e+6+x1.^2.*x2.*4.26544e+5+x1.*x2.^3.*2.192712e+6-x1.^3.*x2.*9.14456e+5+x1.*x2.^4.*8.23048e+5-x1.^4.*x2.*5.97748e+5+x1.*x2.^5.*1.722e+5-x1.^5.*x2.*1.39216e+5+x1.*x2.^6.*1.8666e+4-x1.^6.*x2.*1.1316e+4+x1.*x2.^7.*8.1e+2+x1.^7.*x2.*7.2e+1-x1.^2.*1.29728e+5-x1.^3.*6.41072e+5+x2.^2.*2.563296e+6-x1.^4.*3.81992e+5+x2.^3.*1.860112e+6-x1.^5.*7.7072e+4+x2.^4.*7.89904e+5+x1.^6.*5.3e+3+x2.^5.*2.01492e+5+x1.^7.*4.224e+3+x2.^6.*3.0168e+4+x1.^8.*4.32e+2+x2.^7.*2.43e+3+x2.^8.*8.1e+1+5.64928e+5).*8.1e+1)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*8.8672e+4+x2.*1.52096e+5+x1.^2.*x2.^2.*1.31664e+5+x1.^2.*x2.^3.*6.6072e+4+x1.^3.*x2.^2.*1.576e+4+x1.^2.*x2.^4.*1.5606e+4+x1.^3.*x2.^3.*1.1916e+4-x1.^4.*x2.^2.*6.816e+3+x1.^2.*x2.^5.*1.629e+3+x1.^3.*x2.^4.*2.52e+3+x1.^4.*x2.^3.*1.17e+2-x1.^5.*x2.^2.*2.106e+3+x1.^2.*x2.^6.*5.4e+1+x1.^3.*x2.^5.*1.62e+2+x1.^4.*x2.^4.*1.08e+2-x1.^5.*x2.^3.*1.08e+2-x1.^6.*x2.^2.*1.62e+2+x1.*x2.*2.32736e+5+x1.*x2.^2.*2.50536e+5+x1.^2.*x2.*1.15128e+5+x1.*x2.^3.*1.32696e+5-x1.^3.*x2.*1.1464e+4+x1.*x2.^4.*3.6756e+4-x1.^4.*x2.*2.716e+4+x1.*x2.^5.*5.058e+3-x1.^5.*x2.*8.958e+3+x1.*x2.^6.*2.7e+2-x1.^6.*x2.*1.197e+3-x1.^7.*x2.*5.4e+1+x1.^2.*4.256e+4-x1.^3.*1.5032e+4+x2.^2.*1.78144e+5-x1.^4.*2.1812e+4+x2.^3.*1.07928e+5-x1.^5.*8.596e+3+x2.^4.*3.666e+4-x1.^6.*1.554e+3+x2.^5.*6.96e+3-x1.^7.*1.08e+2+x2.^6.*6.84e+2+x2.^7.*2.7e+1+5.3152e+4).*8.1e+2)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*2.5536e+4+x2.*2.1872e+4+x1.^2.*x2.^2.*1.0552e+4+x1.^2.*x2.^3.*1.83e+3+x1.^3.*x2.^2.*2.034e+3+x1.^2.*x2.^4.*9.0e+1+x1.^3.*x2.^3.*2.16e+2+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*3.7368e+4+x1.*x2.^2.*2.0396e+4+x1.^2.*x2.*2.3436e+4+x1.*x2.^3.*4.872e+3+x1.^3.*x2.*5.12e+3+x1.*x2.^4.*4.23e+2+x1.^4.*x2.*4.5e+1-x1.^5.*x2.*7.2e+1+x1.^2.*1.7564e+4+x1.^3.*2.996e+3+x2.^2.*1.51e+4-x1.^4.*1.234e+3+x2.^3.*5.028e+3-x1.^5.*5.37e+2+x2.^4.*7.86e+2-x1.^6.*5.4e+1+x2.^5.*4.5e+1+1.2688e+4).*4.32e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*4.9136e+4+x2.*5.1168e+4+x1.^2.*x2.^2.*4.76e+4+x1.^2.*x2.^3.*1.4742e+4+x1.^3.*x2.^2.*1.2622e+4+x1.^2.*x2.^4.*1.929e+3+x1.^3.*x2.^3.*3.102e+3+x1.^4.*x2.^2.*1.05e+3+x1.^2.*x2.^5.*7.2e+1+x1.^3.*x2.^4.*2.34e+2+x1.^4.*x2.^3.*2.16e+2-x1.^5.*x2.^2.*3.6e+1+x1.*x2.*9.4896e+4+x1.*x2.^2.*7.5476e+4+x1.^2.*x2.*6.922e+4+x1.*x2.^3.*2.8916e+4+x1.^3.*x2.*1.9388e+4+x1.*x2.^4.*5.313e+3-x1.^4.*x2.*3.79e+2+x1.*x2.^5.*3.69e+2-x1.^5.*x2.*1.119e+3-x1.^6.*x2.*1.44e+2+x1.^2.*4.156e+4+x1.^3.*1.2268e+4+x2.^2.*4.8376e+4-x1.^4.*2.69e+3+x2.^3.*2.3164e+4-x1.^5.*2.719e+3+x2.^4.*5.89e+3-x1.^6.*6.63e+2+x2.^5.*7.41e+2-x1.^7.*5.4e+1+x2.^6.*3.6e+1+2.24e+4).*2.16e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*4.876e+3+x2.*9.092e+3+x1.^2.*x2.^2.*1.8e+2+x1.*x2.*5.204e+3+x1.*x2.^2.*1.257e+3+x1.^2.*x2.*1.221e+3+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*6.38e+2-x1.^3.*2.85e+2+x2.^2.*3.226e+3-x1.^4.*5.4e+1+x2.^3.*4.23e+2+x2.^4.*1.8e+1+7.9e+3).*3.84e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x1.*1.3388e+4+x2.*1.298e+4+x1.^2.*x2.^2.*1.869e+3+x1.^2.*x2.^3.*1.44e+2+x1.^3.*x2.^2.*1.8e+2+x1.*x2.*1.6168e+4+x1.*x2.^2.*6.067e+3+x1.^2.*x2.*6.353e+3+x1.*x2.^3.*8.19e+2+x1.^3.*x2.*8.25e+2+x1.*x2.^4.*1.8e+1+x1.^2.*5.082e+3-x1.^3.*1.15e+2+x2.^2.*6.43e+3-x1.^4.*4.11e+2+x2.^3.*1.359e+3-x1.^5.*5.4e+1+x2.^4.*9.0e+1+9.448e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).*(x1.*3.1424e+4+x2.*2.016e+4+x1.^2.*x2.^2.*2.1836e+4+x1.^2.*x2.^3.*5.256e+3+x1.^3.*x2.^2.*7.756e+3+x1.^2.*x2.^4.*5.7e+2+x1.^3.*x2.^3.*1.266e+3+x1.^4.*x2.^2.*1.257e+3+x1.^2.*x2.^5.*1.8e+1+x1.^3.*x2.^4.*7.2e+1+x1.^4.*x2.^3.*1.08e+2+x1.^5.*x2.^2.*7.2e+1+x1.*x2.*4.608e+4+x1.*x2.^2.*2.9192e+4+x1.^2.*x2.*4.4224e+4+x1.*x2.^3.*9.452e+3+x1.^3.*x2.*2.1188e+4+x1.*x2.^4.*1.542e+3+x1.^4.*x2.*5.096e+3+x1.*x2.^5.*9.9e+1+x1.^5.*x2.*5.55e+2+x1.^6.*x2.*1.8e+1+x1.^2.*3.748e+4+x1.^3.*2.2984e+4+x2.^2.*1.6168e+4+x1.^4.*7.528e+3+x2.^3.*6.88e+3+x1.^5.*1.262e+3+x2.^4.*1.612e+3+x1.^6.*8.4e+1+x2.^5.*1.92e+2+x2.^6.*9.0+1.1072e+4).*3.24e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).^2.*(x1.*5.0e+3+x2.*3.792e+3+x1.^2.*x2.^2.*5.34e+2+x1.^2.*x2.^3.*3.6e+1+x1.^3.*x2.^2.*5.4e+1+x1.*x2.*5.192e+3+x1.*x2.^2.*1.766e+3+x1.^2.*x2.*2.246e+3+x1.*x2.^3.*2.34e+2+x1.^3.*x2.*4.38e+2+x1.*x2.^4.*9.0+x1.^4.*x2.*3.6e+1+x1.^2.*2.628e+3+x1.^3.*7.14e+2+x2.^2.*1.652e+3+x1.^4.*1.2e+2+x2.^3.*3.06e+2+x1.^5.*9.0+x2.^4.*1.8e+1+3.152e+3).*1.296e+4)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*1.8224e+4+x2.*9.824e+3+x1.^2.*x2.^2.*4.916e+3+x1.^2.*x2.^3.*5.82e+2+x1.^3.*x2.^2.*1.218e+3+x1.^2.*x2.^4.*1.8e+1+x1.^3.*x2.^3.*7.2e+1+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*1.9984e+4+x1.*x2.^2.*8.332e+3+x1.^2.*x2.*1.5908e+4+x1.*x2.^3.*1.512e+3+x1.^3.*x2.*5.88e+3+x1.*x2.^4.*9.9e+1+x1.^4.*x2.*1.041e+3+x1.^5.*x2.*7.2e+1+x1.^2.*1.7896e+4+x1.^3.*8.66e+3+x2.^2.*5.432e+3+x1.^4.*2.274e+3+x2.^3.*1.452e+3+x1.^5.*3.15e+2+x2.^4.*1.86e+2+x1.^6.*1.8e+1+x2.^5.*9.0+7.296e+3).*6.48e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*4.264e+3+x2.*6.752e+3+x1.^2.*x2.^2.*1.08e+2+x1.*x2.*3.944e+3+x1.*x2.^2.*9.33e+2+x1.^2.*x2.*8.61e+2+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*1.07e+3+x1.^3.*1.83e+2+x2.^2.*2.362e+3+x1.^4.*1.8e+1+x2.^3.*3.51e+2+x2.^4.*1.8e+1+6.352e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).^2.*(x1.*1.688e+3+x2.*2.488e+3+x1.^2.*x2.^2.*5.4e+1+x1.*x2.*1.6e+3+x1.*x2.^2.*4.26e+2+x1.^2.*x2.*3.66e+2+x1.*x2.^3.*3.6e+1+x1.^3.*x2.*3.6e+1+x1.^2.*4.96e+2+x1.^3.*1.02e+2+x2.^2.*9.68e+2+x1.^4.*9.0+x2.^3.*1.62e+2+x2.^4.*9.0+2.216e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x1.*(x2+1.0).*(x1.*1.1096e+4+x2.*8.888e+3+x1.^2.*x2.^2.*1.122e+3+x1.^2.*x2.^3.*7.2e+1+x1.^3.*x2.^2.*1.08e+2+x1.*x2.*1.1488e+4+x1.*x2.^2.*3.904e+3+x1.^2.*x2.*4.772e+3+x1.*x2.^3.*5.04e+2+x1.^3.*x2.*9.12e+2+x1.*x2.^4.*1.8e+1+x1.^4.*x2.*7.2e+1+x1.^2.*5.7e+3+x1.^3.*1.544e+3+x2.^2.*3.988e+3+x1.^4.*2.49e+2+x2.^3.*7.56e+2+x1.^5.*1.8e+1+x2.^4.*4.5e+1+7.24e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8);
EquationTwo = @(x1,x2) ((x1.*6.6912e+5+x2.*1.879456e+6+x1.^2.*x2.^2.*1.204072e+6+x1.^2.*x2.^3.*9.70736e+5-x1.^3.*x2.^2.*2.95528e+5+x1.^2.*x2.^4.*3.58764e+5+x1.^3.*x2.^3.*7.1352e+4-x1.^4.*x2.^2.*2.8628e+5+x1.^2.*x2.^5.*6.6096e+4+x1.^3.*x2.^4.*5.3724e+4-x1.^4.*x2.^3.*4.5354e+4-x1.^5.*x2.^2.*6.447e+4+x1.^2.*x2.^6.*5.643e+3+x1.^3.*x2.^5.*9.252e+3-x1.^4.*x2.^4.*9.0-x1.^5.*x2.^3.*9.558e+3-x1.^6.*x2.^2.*5.571e+3+x1.^2.*x2.^7.*1.62e+2+x1.^3.*x2.^6.*4.86e+2+x1.^4.*x2.^5.*3.24e+2-x1.^5.*x2.^4.*3.24e+2-x1.^6.*x2.^3.*4.86e+2-x1.^7.*x2.^2.*1.62e+2+x1.*x2.*2.355808e+6+x1.*x2.^2.*3.214896e+6+x1.^2.*x2.*4.26544e+5+x1.*x2.^3.*2.192712e+6-x1.^3.*x2.*9.14456e+5+x1.*x2.^4.*8.23048e+5-x1.^4.*x2.*5.97748e+5+x1.*x2.^5.*1.722e+5-x1.^5.*x2.*1.39216e+5+x1.*x2.^6.*1.8666e+4-x1.^6.*x2.*1.1316e+4+x1.*x2.^7.*8.1e+2+x1.^7.*x2.*7.2e+1-x1.^2.*1.29728e+5-x1.^3.*6.41072e+5+x2.^2.*2.563296e+6-x1.^4.*3.81992e+5+x2.^3.*1.860112e+6-x1.^5.*7.7072e+4+x2.^4.*7.89904e+5+x1.^6.*5.3e+3+x2.^5.*2.01492e+5+x1.^7.*4.224e+3+x2.^6.*3.0168e+4+x1.^8.*4.32e+2+x2.^7.*2.43e+3+x2.^8.*8.1e+1+5.64928e+5).*8.1e+1)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*2.5536e+4+x2.*2.1872e+4+x1.^2.*x2.^2.*1.0552e+4+x1.^2.*x2.^3.*1.83e+3+x1.^3.*x2.^2.*2.034e+3+x1.^2.*x2.^4.*9.0e+1+x1.^3.*x2.^3.*2.16e+2+x1.^4.*x2.^2.*1.08e+2+x1.*x2.*3.7368e+4+x1.*x2.^2.*2.0396e+4+x1.^2.*x2.*2.3436e+4+x1.*x2.^3.*4.872e+3+x1.^3.*x2.*5.12e+3+x1.*x2.^4.*4.23e+2+x1.^4.*x2.*4.5e+1-x1.^5.*x2.*7.2e+1+x1.^2.*1.7564e+4+x1.^3.*2.996e+3+x2.^2.*1.51e+4-x1.^4.*1.234e+3+x2.^3.*5.028e+3-x1.^5.*5.37e+2+x2.^4.*7.86e+2-x1.^6.*5.4e+1+x2.^5.*4.5e+1+1.2688e+4).*4.32e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*4.9136e+4+x2.*5.1168e+4+x1.^2.*x2.^2.*4.76e+4+x1.^2.*x2.^3.*1.4742e+4+x1.^3.*x2.^2.*1.2622e+4+x1.^2.*x2.^4.*1.929e+3+x1.^3.*x2.^3.*3.102e+3+x1.^4.*x2.^2.*1.05e+3+x1.^2.*x2.^5.*7.2e+1+x1.^3.*x2.^4.*2.34e+2+x1.^4.*x2.^3.*2.16e+2-x1.^5.*x2.^2.*3.6e+1+x1.*x2.*9.4896e+4+x1.*x2.^2.*7.5476e+4+x1.^2.*x2.*6.922e+4+x1.*x2.^3.*2.8916e+4+x1.^3.*x2.*1.9388e+4+x1.*x2.^4.*5.313e+3-x1.^4.*x2.*3.79e+2+x1.*x2.^5.*3.69e+2-x1.^5.*x2.*1.119e+3-x1.^6.*x2.*1.44e+2+x1.^2.*4.156e+4+x1.^3.*1.2268e+4+x2.^2.*4.8376e+4-x1.^4.*2.69e+3+x2.^3.*2.3164e+4-x1.^5.*2.719e+3+x2.^4.*5.89e+3-x1.^6.*6.63e+2+x2.^5.*7.41e+2-x1.^7.*5.4e+1+x2.^6.*3.6e+1+2.24e+4).*2.16e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*4.876e+3+x2.*9.092e+3+x1.^2.*x2.^2.*1.8e+2+x1.*x2.*5.204e+3+x1.*x2.^2.*1.257e+3+x1.^2.*x2.*1.221e+3+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*6.38e+2-x1.^3.*2.85e+2+x2.^2.*3.226e+3-x1.^4.*5.4e+1+x2.^3.*4.23e+2+x2.^4.*1.8e+1+7.9e+3).*3.84e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x1.*1.3388e+4+x2.*1.298e+4+x1.^2.*x2.^2.*1.869e+3+x1.^2.*x2.^3.*1.44e+2+x1.^3.*x2.^2.*1.8e+2+x1.*x2.*1.6168e+4+x1.*x2.^2.*6.067e+3+x1.^2.*x2.*6.353e+3+x1.*x2.^3.*8.19e+2+x1.^3.*x2.*8.25e+2+x1.*x2.^4.*1.8e+1+x1.^2.*5.082e+3-x1.^3.*1.15e+2+x2.^2.*6.43e+3-x1.^4.*4.11e+2+x2.^3.*1.359e+3-x1.^5.*5.4e+1+x2.^4.*9.0e+1+9.448e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).*(x1.*3.1424e+4+x2.*2.016e+4+x1.^2.*x2.^2.*2.1836e+4+x1.^2.*x2.^3.*5.256e+3+x1.^3.*x2.^2.*7.756e+3+x1.^2.*x2.^4.*5.7e+2+x1.^3.*x2.^3.*1.266e+3+x1.^4.*x2.^2.*1.257e+3+x1.^2.*x2.^5.*1.8e+1+x1.^3.*x2.^4.*7.2e+1+x1.^4.*x2.^3.*1.08e+2+x1.^5.*x2.^2.*7.2e+1+x1.*x2.*4.608e+4+x1.*x2.^2.*2.9192e+4+x1.^2.*x2.*4.4224e+4+x1.*x2.^3.*9.452e+3+x1.^3.*x2.*2.1188e+4+x1.*x2.^4.*1.542e+3+x1.^4.*x2.*5.096e+3+x1.*x2.^5.*9.9e+1+x1.^5.*x2.*5.55e+2+x1.^6.*x2.*1.8e+1+x1.^2.*3.748e+4+x1.^3.*2.2984e+4+x2.^2.*1.6168e+4+x1.^4.*7.528e+3+x2.^3.*6.88e+3+x1.^5.*1.262e+3+x2.^4.*1.612e+3+x1.^6.*8.4e+1+x2.^5.*1.92e+2+x2.^6.*9.0+1.1072e+4).*3.24e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)+((x2+1.0).^2.*(x1.*5.0e+3+x2.*3.792e+3+x1.^2.*x2.^2.*5.34e+2+x1.^2.*x2.^3.*3.6e+1+x1.^3.*x2.^2.*5.4e+1+x1.*x2.*5.192e+3+x1.*x2.^2.*1.766e+3+x1.^2.*x2.*2.246e+3+x1.*x2.^3.*2.34e+2+x1.^3.*x2.*4.38e+2+x1.*x2.^4.*9.0+x1.^4.*x2.*3.6e+1+x1.^2.*2.628e+3+x1.^3.*7.14e+2+x2.^2.*1.652e+3+x1.^4.*1.2e+2+x2.^3.*3.06e+2+x1.^5.*9.0+x2.^4.*1.8e+1+3.152e+3).*1.296e+4)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x2+1.0).*(x1.*4.264e+3+x2.*6.752e+3+x1.^2.*x2.^2.*1.08e+2+x1.*x2.*3.944e+3+x1.*x2.^2.*9.33e+2+x1.^2.*x2.*8.61e+2+x1.*x2.^3.*7.2e+1+x1.^3.*x2.*7.2e+1+x1.^2.*1.07e+3+x1.^3.*1.83e+2+x2.^2.*2.362e+3+x1.^4.*1.8e+1+x2.^3.*3.51e+2+x2.^4.*1.8e+1+6.352e+3).*5.76e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8)-(x2.*(x2+1.0).*(x1.*1.1096e+4+x2.*8.888e+3+x1.^2.*x2.^2.*1.122e+3+x1.^2.*x2.^3.*7.2e+1+x1.^3.*x2.^2.*1.08e+2+x1.*x2.*1.1488e+4+x1.*x2.^2.*3.904e+3+x1.^2.*x2.*4.772e+3+x1.*x2.^3.*5.04e+2+x1.^3.*x2.*9.12e+2+x1.*x2.^4.*1.8e+1+x1.^4.*x2.*7.2e+1+x1.^2.*5.7e+3+x1.^3.*1.544e+3+x2.^2.*3.988e+3+x1.^4.*2.49e+2+x2.^3.*7.56e+2+x1.^5.*1.8e+1+x2.^4.*4.5e+1+7.24e+3).*8.64e+3)./(x1.*2.8601024e+8+x2.*4.02569632e+8+x1.^2.*x2.^2.*2.93262264e+8+x1.^2.*x2.^3.*1.18193472e+8+x1.^3.*x2.^2.*6.5908584e+7+x1.^2.*x2.^4.*2.5987608e+7+x1.^3.*x2.^3.*2.2991904e+7+x1.^4.*x2.^2.*1.74744e+6+x1.^2.*x2.^5.*3.105702e+6+x1.^3.*x2.^4.*4.133268e+6+x1.^4.*x2.^3.*8.25552e+5-x1.^5.*x2.^2.*1.46367e+6+x1.^2.*x2.^6.*1.86381e+5+x1.^3.*x2.^5.*3.71304e+5+x1.^4.*x2.^4.*1.45557e+5-x1.^5.*x2.^3.*2.09466e+5-x1.^6.*x2.^2.*1.74717e+5+x1.^2.*x2.^7.*4.374e+3+x1.^3.*x2.^6.*1.3122e+4+x1.^4.*x2.^5.*8.748e+3-x1.^5.*x2.^4.*8.748e+3-x1.^6.*x2.^3.*1.3122e+4-x1.^7.*x2.^2.*4.374e+3+x1.*x2.*6.13279456e+8+x1.*x2.^2.*5.36958912e+8+x1.^2.*x2.*3.71391648e+8+x1.*x2.^3.*2.48018184e+8+x1.^3.*x2.*9.0650088e+7+x1.*x2.^4.*6.4734336e+7-x1.^4.*x2.*6.82956e+5+x1.*x2.^5.*9.51102e+6-x1.^5.*x2.*4.158972e+6+x1.*x2.^6.*7.22682e+5-x1.^6.*x2.*5.83362e+5+x1.*x2.^7.*2.187e+4-x1.^7.*x2.*1.2636e+4+x1.^2.*1.84163584e+8+x1.^3.*4.2898656e+7+x2.^2.*3.94143392e+8-x1.^4.*6.436944e+6+x2.^3.*2.10249984e+8-x1.^5.*5.102424e+6+x2.^4.*6.6124008e+7-x1.^6.*7.02e+5+x2.^5.*1.2447324e+7+x1.^7.*4.6008e+4+x2.^6.*1.365336e+6+x1.^8.*1.1664e+4+x2.^7.*8.262e+4+x2.^8.*2.187e+3+1.73017856e+8); matlab function MATLAB Answers — New Questions
How to delay start or reset time in Signal Editor?
In my simulation a capacitor is getting discharged to the resistor.
When the capacitor voltage reaches a predetermined threshold, I need to start variable Duty Cycle PWM.
Duty Cycle value is coming from the Signal Editor.
How can I do that?
In other words – how can I start this circuitry not at time 0, but when the condition is met?In my simulation a capacitor is getting discharged to the resistor.
When the capacitor voltage reaches a predetermined threshold, I need to start variable Duty Cycle PWM.
Duty Cycle value is coming from the Signal Editor.
How can I do that?
In other words – how can I start this circuitry not at time 0, but when the condition is met? In my simulation a capacitor is getting discharged to the resistor.
When the capacitor voltage reaches a predetermined threshold, I need to start variable Duty Cycle PWM.
Duty Cycle value is coming from the Signal Editor.
How can I do that?
In other words – how can I start this circuitry not at time 0, but when the condition is met? signal editor, time reset MATLAB Answers — New Questions
Combining Scripts with AppDesigner
I am using Matlab to control some linear stages and a scope. the application needs to run in a couple of modes, 1) manual mode (for debugging setup) and 2) automated mode (for functional tests). I would normally use a "main loop" (ex. main.m) to manage hardware activities and a GUI to set parameters and display results. How do I do this AppDesign? I’ve not seen a way to do my normal approach, is there a better approach for AppDesigner.I am using Matlab to control some linear stages and a scope. the application needs to run in a couple of modes, 1) manual mode (for debugging setup) and 2) automated mode (for functional tests). I would normally use a "main loop" (ex. main.m) to manage hardware activities and a GUI to set parameters and display results. How do I do this AppDesign? I’ve not seen a way to do my normal approach, is there a better approach for AppDesigner. I am using Matlab to control some linear stages and a scope. the application needs to run in a couple of modes, 1) manual mode (for debugging setup) and 2) automated mode (for functional tests). I would normally use a "main loop" (ex. main.m) to manage hardware activities and a GUI to set parameters and display results. How do I do this AppDesign? I’ve not seen a way to do my normal approach, is there a better approach for AppDesigner. m scripts, appdesigner MATLAB Answers — New Questions
MATLAB does not detect NVIDIA T400 GPU despite updated drivers
Hi everyone,
I am experiencing an issue with MATLAB (version R2024a) on Windows 11 Pro where it fails to detect my NVIDIA T400 4GB GPU. I have already followed several troubleshooting steps without success. Here are the details of my system and the steps I’ve taken so far:
System Details:
Operating System: Windows 11 Pro Version 10.0 (Build 22631)
MATLAB Version: R2024a (Update 3)
GPU: NVIDIA T400 4GB
NVIDIA Driver Version: 516.94
CUDA Version: 11.7
Steps Taken:
Updated the NVIDIA drivers to the latest version (516.94).
Verified that the GPU is correctly recognized by the operating system using nvidia-smi:bash
C:UsersUtente>nvidia-smi Wed May 29 12:02:48 2024
+—————————————————————————–+
| NVIDIA-SMI 516.94 Driver Version: 516.94 CUDA Version: 11.7 |
|——————————-+———————-+———————-+
| GPU Name TCC/WDDM | Bus-Id Disp.A | Volatile Uncorr. ECC |
| Fan Temp Perf Pwr:Usage/Cap| Memory-Usage | GPU-Util Compute M. |
| | | MIG M. |
|===============================+======================+======================|
| 0 NVIDIA T400 4GB WDDM | 00000000:01:00.0 Off | N/A |
| 38% 37C P8 N/A / 31W | 1MiB / 4096MiB | 0% Default |
| | | N/A |
+——————————-+———————-+———————-+
+—————————————————————————–+
| Processes: |
| GPU GI CI PID Type Process name GPU Memory |
| ID ID Usage |
|=============================================================================|
| 0 N/A N/A 10668 C+G …024abinwin64MATLAB.exe N/A |
+—————————————————————————–+
Ran MATLAB as an administrator to ensure it has the necessary permissions.
Checked that the CUDA environment variables are correctly set.
Set the Windows power plan to "High performance".
Ran gpuDeviceCount in MATLAB, but it does not detect any available GPU = gpuDeviceCount; disp([‘Number of GPUs available: ‘, num2str(n)]); % Output: Number of GPUs available: 0
Question:What can I do to resolve this issue and have MATLAB correctly detect my NVIDIA GPU? Are there any additional configuration steps I should follow or any other settings to check? Any suggestions or similar experiences would be greatly appreciated.
Thanks in advance for your help!Hi everyone,
I am experiencing an issue with MATLAB (version R2024a) on Windows 11 Pro where it fails to detect my NVIDIA T400 4GB GPU. I have already followed several troubleshooting steps without success. Here are the details of my system and the steps I’ve taken so far:
System Details:
Operating System: Windows 11 Pro Version 10.0 (Build 22631)
MATLAB Version: R2024a (Update 3)
GPU: NVIDIA T400 4GB
NVIDIA Driver Version: 516.94
CUDA Version: 11.7
Steps Taken:
Updated the NVIDIA drivers to the latest version (516.94).
Verified that the GPU is correctly recognized by the operating system using nvidia-smi:bash
C:UsersUtente>nvidia-smi Wed May 29 12:02:48 2024
+—————————————————————————–+
| NVIDIA-SMI 516.94 Driver Version: 516.94 CUDA Version: 11.7 |
|——————————-+———————-+———————-+
| GPU Name TCC/WDDM | Bus-Id Disp.A | Volatile Uncorr. ECC |
| Fan Temp Perf Pwr:Usage/Cap| Memory-Usage | GPU-Util Compute M. |
| | | MIG M. |
|===============================+======================+======================|
| 0 NVIDIA T400 4GB WDDM | 00000000:01:00.0 Off | N/A |
| 38% 37C P8 N/A / 31W | 1MiB / 4096MiB | 0% Default |
| | | N/A |
+——————————-+———————-+———————-+
+—————————————————————————–+
| Processes: |
| GPU GI CI PID Type Process name GPU Memory |
| ID ID Usage |
|=============================================================================|
| 0 N/A N/A 10668 C+G …024abinwin64MATLAB.exe N/A |
+—————————————————————————–+
Ran MATLAB as an administrator to ensure it has the necessary permissions.
Checked that the CUDA environment variables are correctly set.
Set the Windows power plan to "High performance".
Ran gpuDeviceCount in MATLAB, but it does not detect any available GPU = gpuDeviceCount; disp([‘Number of GPUs available: ‘, num2str(n)]); % Output: Number of GPUs available: 0
Question:What can I do to resolve this issue and have MATLAB correctly detect my NVIDIA GPU? Are there any additional configuration steps I should follow or any other settings to check? Any suggestions or similar experiences would be greatly appreciated.
Thanks in advance for your help! Hi everyone,
I am experiencing an issue with MATLAB (version R2024a) on Windows 11 Pro where it fails to detect my NVIDIA T400 4GB GPU. I have already followed several troubleshooting steps without success. Here are the details of my system and the steps I’ve taken so far:
System Details:
Operating System: Windows 11 Pro Version 10.0 (Build 22631)
MATLAB Version: R2024a (Update 3)
GPU: NVIDIA T400 4GB
NVIDIA Driver Version: 516.94
CUDA Version: 11.7
Steps Taken:
Updated the NVIDIA drivers to the latest version (516.94).
Verified that the GPU is correctly recognized by the operating system using nvidia-smi:bash
C:UsersUtente>nvidia-smi Wed May 29 12:02:48 2024
+—————————————————————————–+
| NVIDIA-SMI 516.94 Driver Version: 516.94 CUDA Version: 11.7 |
|——————————-+———————-+———————-+
| GPU Name TCC/WDDM | Bus-Id Disp.A | Volatile Uncorr. ECC |
| Fan Temp Perf Pwr:Usage/Cap| Memory-Usage | GPU-Util Compute M. |
| | | MIG M. |
|===============================+======================+======================|
| 0 NVIDIA T400 4GB WDDM | 00000000:01:00.0 Off | N/A |
| 38% 37C P8 N/A / 31W | 1MiB / 4096MiB | 0% Default |
| | | N/A |
+——————————-+———————-+———————-+
+—————————————————————————–+
| Processes: |
| GPU GI CI PID Type Process name GPU Memory |
| ID ID Usage |
|=============================================================================|
| 0 N/A N/A 10668 C+G …024abinwin64MATLAB.exe N/A |
+—————————————————————————–+
Ran MATLAB as an administrator to ensure it has the necessary permissions.
Checked that the CUDA environment variables are correctly set.
Set the Windows power plan to "High performance".
Ran gpuDeviceCount in MATLAB, but it does not detect any available GPU = gpuDeviceCount; disp([‘Number of GPUs available: ‘, num2str(n)]); % Output: Number of GPUs available: 0
Question:What can I do to resolve this issue and have MATLAB correctly detect my NVIDIA GPU? Are there any additional configuration steps I should follow or any other settings to check? Any suggestions or similar experiences would be greatly appreciated.
Thanks in advance for your help! gpu MATLAB Answers — New Questions
what meaning “We are unable to offer you a trial. If you believe you are receiving this message in error, please”
what meaning
"We are unable to offer you a trial. If you believe you are receiving this message in error, please"what meaning
"We are unable to offer you a trial. If you believe you are receiving this message in error, please" what meaning
"We are unable to offer you a trial. If you believe you are receiving this message in error, please" MATLAB Answers — New Questions
Voronoi boundaries
Hi guys
the program voronoiDiagram gives the vertices of a voronoi decomposition of a set of scattered datapoints. These vertices do, however, extend beyond the rectangular domain containing the datapoints. Does anyone know of a routine/method/program which allows the vertices to extend no further than the boundary of the domain?
I have been trying, in vain, to set up my own routine based on the voronoiDiagram output, but it has become a bit more unmanageable than I had foreseen.
Thanks in advance!Hi guys
the program voronoiDiagram gives the vertices of a voronoi decomposition of a set of scattered datapoints. These vertices do, however, extend beyond the rectangular domain containing the datapoints. Does anyone know of a routine/method/program which allows the vertices to extend no further than the boundary of the domain?
I have been trying, in vain, to set up my own routine based on the voronoiDiagram output, but it has become a bit more unmanageable than I had foreseen.
Thanks in advance! Hi guys
the program voronoiDiagram gives the vertices of a voronoi decomposition of a set of scattered datapoints. These vertices do, however, extend beyond the rectangular domain containing the datapoints. Does anyone know of a routine/method/program which allows the vertices to extend no further than the boundary of the domain?
I have been trying, in vain, to set up my own routine based on the voronoiDiagram output, but it has become a bit more unmanageable than I had foreseen.
Thanks in advance! voronoi, voronoidiagram, vertices, domain MATLAB Answers — New Questions
Regarding Matlab Simulink – R2019 Version
I have a Simulink model from R2019b and now I want to upgrade my MATLAB to R2023b, so I want to know how to use the previous version Simulink model in latest MATLAB, are there any steps to follow, I tried to check but couldn’t get any clarity on this. Any help is greatly appreciated.I have a Simulink model from R2019b and now I want to upgrade my MATLAB to R2023b, so I want to know how to use the previous version Simulink model in latest MATLAB, are there any steps to follow, I tried to check but couldn’t get any clarity on this. Any help is greatly appreciated. I have a Simulink model from R2019b and now I want to upgrade my MATLAB to R2023b, so I want to know how to use the previous version Simulink model in latest MATLAB, are there any steps to follow, I tried to check but couldn’t get any clarity on this. Any help is greatly appreciated. matlab, simulink MATLAB Answers — New Questions
changes the the simulated mass with changing the absolute and relative tolerance
Hi,
The mass of the tissue under the effect of a drug goes to very small values, i.e. in the order 10^-12, which is unrealistic. and then the tissue regrow, which is also unrealistic because the smallest mass of a cell is in the order of 10^-9. What should i do about that?
Also, if I increase the absulote and relative tolerance, the simulation results change and this is unreliable. How can I have more reliable results?
Any thoughts or suggestions are appreciated.
ThanksHi,
The mass of the tissue under the effect of a drug goes to very small values, i.e. in the order 10^-12, which is unrealistic. and then the tissue regrow, which is also unrealistic because the smallest mass of a cell is in the order of 10^-9. What should i do about that?
Also, if I increase the absulote and relative tolerance, the simulation results change and this is unreliable. How can I have more reliable results?
Any thoughts or suggestions are appreciated.
Thanks Hi,
The mass of the tissue under the effect of a drug goes to very small values, i.e. in the order 10^-12, which is unrealistic. and then the tissue regrow, which is also unrealistic because the smallest mass of a cell is in the order of 10^-9. What should i do about that?
Also, if I increase the absulote and relative tolerance, the simulation results change and this is unreliable. How can I have more reliable results?
Any thoughts or suggestions are appreciated.
Thanks simulation results change with changing settings MATLAB Answers — New Questions
I run the code below and expect to get a uniform histogram. It’s not. I don’t understand why not.
% I have a two column array, the first column is time and the second column is state. I ran interp1 to observe my data at uniform time points. The results were unexpected. In trying to identify the problem, I have run this simpler code expecting the deterministic uniformly spaced points to be uniformly distributed. They are not. Can anyone explain why not? Thanks.
% Parameters
thismesh = 100; % Number of bins
deltat = 1 / thismesh;
lasttime = 1000; % Example last time
% Generate uniform time steps
xx = 0:deltat:lasttime;
% Calculate fractional times
fractime = (xx – floor(xx)).’;
% Plot histogram to check uniformity
figure;
histogram(fractime, thismesh);
title(‘Histogram of fractional times’);
xlabel(‘Fractional time’);
ylabel(‘Count’);
% Display counts in each bin
[counts, edges] = histcounts(fractime, thismesh);
disp(‘Counts per bin:’);
disp(counts);
% Expected count per bin
expected_count = length(fractime) / thismesh;
disp(‘Expected count per bin:’);
disp(expected_count);
% Plot actual counts vs expected uniform count
figure;
bar(counts);
hold on;
yline(expected_count, ‘r’, ‘LineWidth’, 2);
hold off;
title(‘Actual counts vs Expected uniform count’);
xlabel(‘Bin’);
ylabel(‘Count’);
legend(‘Actual counts’, ‘Expected count’);
[SL: formatted code as code and ran that code]% I have a two column array, the first column is time and the second column is state. I ran interp1 to observe my data at uniform time points. The results were unexpected. In trying to identify the problem, I have run this simpler code expecting the deterministic uniformly spaced points to be uniformly distributed. They are not. Can anyone explain why not? Thanks.
% Parameters
thismesh = 100; % Number of bins
deltat = 1 / thismesh;
lasttime = 1000; % Example last time
% Generate uniform time steps
xx = 0:deltat:lasttime;
% Calculate fractional times
fractime = (xx – floor(xx)).’;
% Plot histogram to check uniformity
figure;
histogram(fractime, thismesh);
title(‘Histogram of fractional times’);
xlabel(‘Fractional time’);
ylabel(‘Count’);
% Display counts in each bin
[counts, edges] = histcounts(fractime, thismesh);
disp(‘Counts per bin:’);
disp(counts);
% Expected count per bin
expected_count = length(fractime) / thismesh;
disp(‘Expected count per bin:’);
disp(expected_count);
% Plot actual counts vs expected uniform count
figure;
bar(counts);
hold on;
yline(expected_count, ‘r’, ‘LineWidth’, 2);
hold off;
title(‘Actual counts vs Expected uniform count’);
xlabel(‘Bin’);
ylabel(‘Count’);
legend(‘Actual counts’, ‘Expected count’);
[SL: formatted code as code and ran that code] % I have a two column array, the first column is time and the second column is state. I ran interp1 to observe my data at uniform time points. The results were unexpected. In trying to identify the problem, I have run this simpler code expecting the deterministic uniformly spaced points to be uniformly distributed. They are not. Can anyone explain why not? Thanks.
% Parameters
thismesh = 100; % Number of bins
deltat = 1 / thismesh;
lasttime = 1000; % Example last time
% Generate uniform time steps
xx = 0:deltat:lasttime;
% Calculate fractional times
fractime = (xx – floor(xx)).’;
% Plot histogram to check uniformity
figure;
histogram(fractime, thismesh);
title(‘Histogram of fractional times’);
xlabel(‘Fractional time’);
ylabel(‘Count’);
% Display counts in each bin
[counts, edges] = histcounts(fractime, thismesh);
disp(‘Counts per bin:’);
disp(counts);
% Expected count per bin
expected_count = length(fractime) / thismesh;
disp(‘Expected count per bin:’);
disp(expected_count);
% Plot actual counts vs expected uniform count
figure;
bar(counts);
hold on;
yline(expected_count, ‘r’, ‘LineWidth’, 2);
hold off;
title(‘Actual counts vs Expected uniform count’);
xlabel(‘Bin’);
ylabel(‘Count’);
legend(‘Actual counts’, ‘Expected count’);
[SL: formatted code as code and ran that code] histogram, uniform distribution, deterministic MATLAB Answers — New Questions
what is i in the answer?
Hello,
I don’t know what is i in the answer of my script. Complex number? how can I simplify it? I need to get rid of i.
syms x(t) s
Dx = diff(x(t),t);
laplace(sin(t)*Dx,t,s)
ans =
– (s*laplace(x(t), t, s – 1i)*1i)/2 + (s*laplace(x(t), t, s + 1i)*1i)/2 – laplace(x(t), t, s – 1i)/2 – laplace(x(t), t, s + 1i)/2
Thank you for helping meHello,
I don’t know what is i in the answer of my script. Complex number? how can I simplify it? I need to get rid of i.
syms x(t) s
Dx = diff(x(t),t);
laplace(sin(t)*Dx,t,s)
ans =
– (s*laplace(x(t), t, s – 1i)*1i)/2 + (s*laplace(x(t), t, s + 1i)*1i)/2 – laplace(x(t), t, s – 1i)/2 – laplace(x(t), t, s + 1i)/2
Thank you for helping me Hello,
I don’t know what is i in the answer of my script. Complex number? how can I simplify it? I need to get rid of i.
syms x(t) s
Dx = diff(x(t),t);
laplace(sin(t)*Dx,t,s)
ans =
– (s*laplace(x(t), t, s – 1i)*1i)/2 + (s*laplace(x(t), t, s + 1i)*1i)/2 – laplace(x(t), t, s – 1i)/2 – laplace(x(t), t, s + 1i)/2
Thank you for helping me laplace, complex number, differential equations, differential, function, laplace transform MATLAB Answers — New Questions
SoC Blockset model not running on ZCU216
I am using SoC Blockset to create a model to run on the Xilinx Zynq Ultrascale+ ZCU216 RFSoC board. I have simulated and run this example on the board, so I know that the DACs and ADCs works: Transmit and Receive Tone Using Xilinx RFSoC Device – Part 2 Deployment – MATLAB & Simulink (mathworks.com)
I created a new Simulink model using the RFSoC template under SoC Blockset. I added my logic (generating a tone on the DAC and reading it on the ADC) and was able to simulate successfully. I am able to also build successfully in external mode, but when I load the bitstream, I do not see any output on my scope. The counter does not go up and no output shows up. Does anyone know what is wrong?I am using SoC Blockset to create a model to run on the Xilinx Zynq Ultrascale+ ZCU216 RFSoC board. I have simulated and run this example on the board, so I know that the DACs and ADCs works: Transmit and Receive Tone Using Xilinx RFSoC Device – Part 2 Deployment – MATLAB & Simulink (mathworks.com)
I created a new Simulink model using the RFSoC template under SoC Blockset. I added my logic (generating a tone on the DAC and reading it on the ADC) and was able to simulate successfully. I am able to also build successfully in external mode, but when I load the bitstream, I do not see any output on my scope. The counter does not go up and no output shows up. Does anyone know what is wrong? I am using SoC Blockset to create a model to run on the Xilinx Zynq Ultrascale+ ZCU216 RFSoC board. I have simulated and run this example on the board, so I know that the DACs and ADCs works: Transmit and Receive Tone Using Xilinx RFSoC Device – Part 2 Deployment – MATLAB & Simulink (mathworks.com)
I created a new Simulink model using the RFSoC template under SoC Blockset. I added my logic (generating a tone on the DAC and reading it on the ADC) and was able to simulate successfully. I am able to also build successfully in external mode, but when I load the bitstream, I do not see any output on my scope. The counter does not go up and no output shows up. Does anyone know what is wrong? xilinx, zynq, simulink, soc blockset MATLAB Answers — New Questions
“Error using uislider (line 53) ‘range’ is not a valid STYLE for uislider. STYLE must be ‘slider’ or ‘range’.”
This error seems like it contradicts itself… it’s coming from my attempt to create a uislider in the following line:
sld = uislider(g,"range","Limits",[0 1], "Value",4/9, "ValueChangedFcn",@(src,event)updateFractal(src,event,iterations, pattern,resolution,ax,ax2));
(of course, ignore the function call, unless it could somehow be relevant)
I don’t see why I would be unable to create a range slider in this instance. At the very least, the error message is very unclear. It works as a normal slider, but not as a range slider.This error seems like it contradicts itself… it’s coming from my attempt to create a uislider in the following line:
sld = uislider(g,"range","Limits",[0 1], "Value",4/9, "ValueChangedFcn",@(src,event)updateFractal(src,event,iterations, pattern,resolution,ax,ax2));
(of course, ignore the function call, unless it could somehow be relevant)
I don’t see why I would be unable to create a range slider in this instance. At the very least, the error message is very unclear. It works as a normal slider, but not as a range slider. This error seems like it contradicts itself… it’s coming from my attempt to create a uislider in the following line:
sld = uislider(g,"range","Limits",[0 1], "Value",4/9, "ValueChangedFcn",@(src,event)updateFractal(src,event,iterations, pattern,resolution,ax,ax2));
(of course, ignore the function call, unless it could somehow be relevant)
I don’t see why I would be unable to create a range slider in this instance. At the very least, the error message is very unclear. It works as a normal slider, but not as a range slider. uislider, gui MATLAB Answers — New Questions
How can I make x and y axis dates with contour?
Hello,
I’m trying to make a Porkchop Plot for a journey from Earth to Mars. The purpose here is to plot the C3 values that will coincide with the date of departure from Earth and the date of arrival on Mars as a contour function. In other words, there should be departure dates from Earth on the x-axis and arrival dates on Mars on the y-axis, but the contour function does not accept input as datetime. How can I solve this problem?Hello,
I’m trying to make a Porkchop Plot for a journey from Earth to Mars. The purpose here is to plot the C3 values that will coincide with the date of departure from Earth and the date of arrival on Mars as a contour function. In other words, there should be departure dates from Earth on the x-axis and arrival dates on Mars on the y-axis, but the contour function does not accept input as datetime. How can I solve this problem? Hello,
I’m trying to make a Porkchop Plot for a journey from Earth to Mars. The purpose here is to plot the C3 values that will coincide with the date of departure from Earth and the date of arrival on Mars as a contour function. In other words, there should be departure dates from Earth on the x-axis and arrival dates on Mars on the y-axis, but the contour function does not accept input as datetime. How can I solve this problem? contour, plot, datetime, xaxis, yaxix, porkchopplot MATLAB Answers — New Questions
How to get details on “Unable to load interface library” in using clib.
We’re building interface modules for a C++ library with clibgen – in different linux setups (server, local PC, cloud). We have a build that works, but not on all setups – and we use Matlab2021a. The error message on computer is:
Unable to load interface library:
‘/home/de11/sim_interfaces/matlab_to_cpp/cpp_lib/WI_MatlabInterface/WI_MatlabInterfaceInterface.so’.
Reason: The specified module could not be found.
Ensure the C++ dependent libraries for the interface library are added to
run-time path.
is telling me something, but not the full picture. As far as i can see all C++ dependent libraries are on the LD_LIBRARY_PATH and calling ldd on the command line tells me all dependencies can be resolved.
But how can i check the same within the matlab session?We’re building interface modules for a C++ library with clibgen – in different linux setups (server, local PC, cloud). We have a build that works, but not on all setups – and we use Matlab2021a. The error message on computer is:
Unable to load interface library:
‘/home/de11/sim_interfaces/matlab_to_cpp/cpp_lib/WI_MatlabInterface/WI_MatlabInterfaceInterface.so’.
Reason: The specified module could not be found.
Ensure the C++ dependent libraries for the interface library are added to
run-time path.
is telling me something, but not the full picture. As far as i can see all C++ dependent libraries are on the LD_LIBRARY_PATH and calling ldd on the command line tells me all dependencies can be resolved.
But how can i check the same within the matlab session? We’re building interface modules for a C++ library with clibgen – in different linux setups (server, local PC, cloud). We have a build that works, but not on all setups – and we use Matlab2021a. The error message on computer is:
Unable to load interface library:
‘/home/de11/sim_interfaces/matlab_to_cpp/cpp_lib/WI_MatlabInterface/WI_MatlabInterfaceInterface.so’.
Reason: The specified module could not be found.
Ensure the C++ dependent libraries for the interface library are added to
run-time path.
is telling me something, but not the full picture. As far as i can see all C++ dependent libraries are on the LD_LIBRARY_PATH and calling ldd on the command line tells me all dependencies can be resolved.
But how can i check the same within the matlab session? clibgen MATLAB Answers — New Questions
