Final.rar_UWF_interpolation

% Guass Seidal Method syms x1 y1 z1 f1 f2 f3 f4 f5 f6 f7 f8 f9, syms a b c fc1 fc2 fc3 , disp( ' Values of 1st Eqs.... '); f1=input('Enter value of coefficient of a :'); f2=input('Enter value of coefficient of b :'); f3=input('Enter value of coefficient of c :'); fc1=input('Enter value of coefficient of B1 :'); disp( ' Enter values of 2nd Eqs : '); f4=input('Enter value of coefficient of a :'); f5=input('Enter value of coefficient of b :'); f6=input('Enter value of coefficient of c :'); fc2=input('Enter value of coefficient of B2 :'); disp( ' Enter values of 3rd Eqs : '); f7=input('Enter value of coefficient of a :'); f8=input('Enter value of coefficient of b :'); f9=input('Enter value of coefficient of c :'); fc3=input('Enter value of coefficient of B3 :'); A=[ f1*a f2*b f3*c ; f4*a f5*b f6*c ; f7*a f8*b f9*c;] x1=(fc1/f1); y1=(fc2-(f4*x1))/f5; % upper value of x1 is used z1=(fc3-(f7*x1)-(f8*y1))/f9; % upper values of x1 & y1 is used for i=1:1:5 a=((1/f1)*(fc1-(f2*y1)-(f3*z1))); % use above values b=((1/f5)*(fc2-(f4*a)-(f6*z1))); % use NEW X-value c=((1/f9)*(fc3-(f7*a)-(f8*b))); % use NEW X&Y-values x1=a; y1=b; z1=c; end disp(' Value of A is : ',a); disp(' Value of B is : ',b); disp(' Value of C is : ',c);