Simultaneous Non-Linear Equation
Simultaneous Non-Linear Equation
clc;
clear all;
close all;
A = [4 3 -4 -3;2 -1 -2 1];
b = [19; 7];
[x A] = gauss_jordan(A,b);
a1 = A(1,:);
h = 1;
for n = 2:2:length(a1)-1
a2(h) = A(2,n);
h = h + 1;
end
a2(1,length(a2)+1) = 0;
a3 = deconv(a1,a2);
a3(1,length(a3)) = A(1,length(A(1,:)))/A(2,length(A(2,:)));
a2(1,length(a2)) = A(2,length(A(2,:)));
% A_ = [[a2(1,1:length(a2)-1);a3(1,1:length(a3)-1)]];
[x A] = gauss_jordan([a2(1,1:length(a2)-1);a3(1,1:length(a3)-
1)],[a2(1,length(a2));a3(1,length(a3))]);
X_ = sqrt(x(1,:));
y_ = x(2);
for n = 1:length(X_)
x(n) = X_(n);
x(n+1) = -X_(n);
y(n,1) = y_(n,1);
y(n+1,1) = y_(n,1);
end
disp(['Solution set:']);
disp(['x = ' num2str(x')])
disp(['y = ' num2str(y')])
%%
Used Function: