Matlab - Matlab - ayuda con planteamiento de función

global var_1;

var_1 = 1;

var_2 = 2;

fun1();

fun2();

function fun1()

global var_1;

fprintf('fun1 %d \n', var_1);

end

function fun2()

global var_2;

fprintf('fun2 %d \n', var_2);

end

fun1 1

fun2

clear, clc

global b1 c1

xlswrite('datos traccion grupo b.xlsx',rand(120,40))

Carga=xlsread('datos traccion grupo b.xlsx', 'C10:C108');

Deformacion=xlsread('datos traccion grupo b.xlsx','B10:B108');

 plot(Deformacion, Carga);

 hold on

 X=xlsread('Datos traccion grupo b.xlsx','B11:B22');

 Y=xlsread('Datos traccion grupo b.xlsx','C11:C22');

 plot(X,Y,'g')

 A=xlsread('datos traccion grupo b.xlsx', 'Q26:X33');

B=xlsread('datos traccion grupo b.xlsx','AB26:AB33');

C=inv(A)*B;

b1=C(1);

c1=C(2);

a2=C(3);

b2=C(4);

c2=C(5);

a3=C(6);

b3=C(7);

c3=C(8);

fprintf('f1(X)=%5.3f   X %5.2f\n',C(1),C(2))

fprintf('f2(X)=%5.3f  X^2+%5.2f  %5.2f  \n',C(3),C(4), C(5))

fprintf('f3(X)=%5.3f  X^2 %5.2f +%5.2f  \n',C(6),C(7), C(8))

disp(' ')

disp('Simpson 3/8 compuesta')

disp('f1')

a=Deformacion(2,1);

b=Deformacion(5,1);

n=input('Ingrese el numero de composiciones:');

h=(b-a)/(3*n);

y=f1(a)

function y=f1(x)

global b1 c1

y=b1*x-c1;

end

>> finalborrador

f1(X)=-0.264   X -0.79

f2(X)=0.684  X^2+ 0.56   0.55

f3(X)=-0.469  X^2 -0.08 + 0.59

Simpson 3/8 compuesta

f1

Ingrese el numero de composiciones:12

y =

    0.6194