Matlab - Como obtener un valor de la expresión F en Matlab

clear all

d=67:0.1:95;

I=1./(d.^2);

n=0:0.1:28;

F=(I+n)/I(1),

clear all

d=67:0.1:95;

I=1./(d.^2);

n=0:0.1:28;

F=[];

for i=1:length(d)

  F=[F;I(i)+n/I(1) ];

end

F

function f

%Matriz de datos

 A = [1.013 0.983 0.953 0.872 0.795 0.727 0.666 0.608 0.557 0.509 0.467 0.427 0.392 0.359 0.329 0.302 0.277 0.254 0.233 0.214 0.196 0.180 0.165 0.152 0.139 0.128 0.118];

% d Distancia

 d = 67:.1:95;

% Inverso cuadrado

 I = 1./(d.^2);

% Factor de Correccion

 n = 0:.1:28;

 F =(I+n)/(I(1));

 a = 5;

 b = 2.5;

 p = 6;

 t = 0.02;

 w = 15;

 x = -14:.1:14;

 f =(abs(x)>w/2).*(t+(0.5-t)*exp((-b/p)*(abs(x)-(w/2))))+(abs(x)<=w/2).*(1-(0.5)*exp((-a/p)*((w/2)-abs(x))))-.5;

 plot(x,f,'g-'),axis([-15 15 -1.2 1.2]);

 hold on

 f1 =(-1)*(abs(x)>w/2).*(t+(0.5-t)*exp((-b/p)*(abs(x)-(w/2))))+(-1)*(abs(x)<=w/2).*(1-(0.5)*exp((-a/p)*((w/2)-abs(x))))+.5;

 plot(x,f1,'r-'),axis([-15 15 -1.2 1.2]);

 hold on

 find (f==f1);

 g2 =((f+n)*A(1,1))/F(2)-1.0;

 plot(x,g2,'ro');