%% xU
% Frase en cuestión
f1 = 'If (x/1000*Cuerda) <= (Cuerda*(NACA2/10)) PartBody\Sketch.1\Offset.numerocota\Offset = (x/1000*Cuerda -NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*sin(atan(2*(NACA1/100)/((NACA2/10)**2)*((NACA2/10)-(x/1000)))))/1000 else PartBody\Sketch.1\Offset.numerocota\Offset = (x/1000*Cuerda -NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*sin(atan(2*(NACA1/100)/((1-(NACA2/10))**2)*((NACA2/10)-(x/1000)))))/1000'; 
%dejamos únicamente lo que será constante y ponemos nombre representativos
%a lo que va a cambiar, en este caso el numero de la cota (numerocota) y la
%coordenada x del perfil (x) esta coordenada va de 0 a 1000, ya que hemos
%decidido que así sea, pero podría estar en cualquier intervalo (de 0 a 1
%por ejemplo).

% Arreglar caracteres conflictivos con fprintf
f2 = strrep(f1, '\', '\\'); % función que sustituye en f1 los \ por \\ ya 
%que la función fprintf usa el comando \ para establecer propiedades como 
%el salto de linea (\n)

% Sustituir en la cadena los elementos que van a variar por %1.0f (ya que
% ninguno de los valores va a tener decimales) que es lo que la función 
% fprintf cambiará por lo que sea necesario.
t = '%1.0f';
f3 = strrep(f2, 'numerocota', t);
f4 = strrep(f3, 'x', t);

% valores por los que deben ser sustituidos los valores por defecto
% anteriores

% el valor numerocota queremos que vaya desde 61 hasta el 90 que es la 
%numeración de las cotas en longitud de los puntos del extradós que hemos 
%colocado en catia

% el valor de las coordenadas en x hemos decidido que tenga una division de
% 20 en los 5 primeros puntos y 36 en los 25 siguiente, para que quede 
% mejor definido el borde de ataque, variando x entre 0 y 1000 que era lo 
% que se pretendía

%el orden de la numeración de las cotas debe coincidir con el orden en el
%que se recorre el vector x
numerocota=61:90; x=[20:20:100,136:36:1000];

% Mostrar en pantalla
fprintf('\n')
for k = 1:length(x);
    C = [x(k),numerocota(k),x(k)*ones(1,7),numerocota(k),x(k)*ones(1,7)]; 
    % vector que contiene los carácteres ordenados que se van a sustituir 
    %en la frase. Tiene un orden especial como puede verse primero aparece
    %la x, despues numerocota, a continuación la x 7 veces etc.
    fprintf([ f4 , ' \n' ], C ) %sustituye los %1.0f por los valores de C
end

%% yU
% Frase en cuestión
f1 = 'If (x/1000*Cuerda) <= (Cuerda*(NACA2/10)) PartBody\Sketch.1\Offset.numerocota\Offset = (Cuerda/10+((NACA1/100)*(x/1000*Cuerda)*(2*(NACA2/10)-x/1000)/(NACA2/10)**2) +NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*cos(atan(2*(NACA1/100)/((NACA2/10)**2)*((NACA2/10)-(x/1000)))))/1000 else PartBody\Sketch.1\Offset.numerocota\Offset = (Cuerda/10+((NACA1/100)*(Cuerda-x/1000*Cuerda)*(-2*(NACA2/10)+x/1000+1)/(1-NACA2/10)**2) +NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*cos(atan(2*(NACA1/100)/((1-(NACA2/10))**2)*((NACA2/10)-(x/1000)))))/1000'; 
%dejamos únicamente lo que será constante y ponemos nombre representativos
%a lo que va a cambiar, en este caso el numero de la cota (numerocota) y la
%coordenada x del perfil (x) esta coordenada va de 0 a 1000, ya que hemos
%decidido que así sea, pero podría estar en cualquier intervalo (de 0 a 1
%por ejemplo).

% Arreglar caracteres conflictivos con fprintf
f2 = strrep(f1, '\', '\\'); % función que sustituye en f1 los \ por \\ ya 
%que la función fprintf usa el comando \ para establecer propiedades como 
%el salto de linea (\n)

% Sustituir en la cadena los elementos que van a variar por %1.0f (ya que
% ninguno de los valores va a tener decimales) que es lo que la función 
% fprintf cambiará por lo que sea necesario.
t = '%1.0f';
f3 = strrep(f2, 'numerocota', t);
f4 = strrep(f3, 'x', t);

% valores por los que deben ser sustituidos los valores por defecto
% anteriores

% el valor numerocota queremos que vaya desde 91 hasta el 120 que es la 
%numeración de las cotas en longitud de los puntos del extradós que hemos 
%colocado en catia

% el valor de las coordenadas en x hemos decidido que tenga una division de
% 20 en los 5 primeros puntos y 36 en los 25 siguiente, para que quede 
% mejor definido el borde de ataque, variando x entre 0 y 1000 que era lo 
% que se pretendía

%el orden de la numeración de las cotas debe coincidir con el orden en el
%que se recorre el vector x
numerocota=91:120; x=[20:20:100,136:36:1000];

% Mostrar en pantalla
fprintf('\n')
for k = 1:length(x);
    C = [x(k),numerocota(k),x(k)*ones(1,8),numerocota(k),x(k)*ones(1,8)]; 
    % vector que contiene los carácteres ordenados que se van a sustituir 
    %en la frase. Tiene un orden especial como puede verse primero aparece
    %la x, despues numerocota, a continuación la x 8 veces etc.
    fprintf([ f4 , ' \n' ], C ) %sustituye los %1.0f por los valores de C
end

%% xL
% Frase en cuestión
f1 = 'If (x/1000*Cuerda) <= (Cuerda*(NACA2/10)) PartBody\Sketch.1\Offset.numerocota\Offset = (x/1000*Cuerda +NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*sin(atan(2*(NACA1/100)/((NACA2/10)**2)*((NACA2/10)-(x/1000)))))/1000 else PartBody\Sketch.1\Offset.numerocota\Offset = (x/1000*Cuerda +NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*sin(atan(2*(NACA1/100)/((1-(NACA2/10))**2)*((NACA2/10)-(x/1000)))))/1000'; 
%dejamos únicamente lo que será constante y ponemos nombre representativos
%a lo que va a cambiar, en este caso el numero de la cota (numerocota) y la
%coordenada x del perfil (x) esta coordenada va de 0 a 1000, ya que hemos
%decidido que así sea, pero podría estar en cualquier intervalo (de 0 a 1
%por ejemplo).

% Arreglar caracteres conflictivos con fprintf
f2 = strrep(f1, '\', '\\'); % función que sustituye en f1 los \ por \\ ya 
%que la función fprintf usa el comando \ para establecer propiedades como 
%el salto de linea (\n)

% Sustituir en la cadena los elementos que van a variar por %1.0f (ya que
% ninguno de los valores va a tener decimales) que es lo que la función 
% fprintf cambiará por lo que sea necesario.
t = '%1.0f';
f3 = strrep(f2, 'numerocota', t);
f4 = strrep(f3, 'x', t);

% valores por los que deben ser sustituidos los valores por defecto
% anteriores

% el valor numerocota queremos que vaya desde 121 hasta el 150 que es la 
%numeración de las cotas en longitud de los puntos del extradós que hemos 
%colocado en catia

% el valor de las coordenadas en x hemos decidido que tenga una division de
% 20 en los 5 primeros puntos y 36 en los 25 siguiente, para que quede 
% mejor definido el borde de ataque, variando x entre 0 y 1000 que era lo 
% que se pretendía

%el orden de la numeración de las cotas debe coincidir con el orden en el
%que se recorre el vector x
numerocota=121:150; x=[20:20:100,136:36:1000];

% Mostrar en pantalla
fprintf('\n')
for k = 1:length(x);
    C = [x(k),numerocota(k),x(k)*ones(1,7),numerocota(k),x(k)*ones(1,7)]; 
    % vector que contiene los carácteres ordenados que se van a sustituir 
    %en la frase. Tiene un orden especial como puede verse primero aparece
    %la x, despues numerocota, a continuación la x 7 veces etc.
    fprintf([ f4 , ' \n' ], C ) %sustituye los %1.0f por los valores de C
end

%% yL
% Frase en cuestión
f1 = 'If (x/1000*Cuerda) <= (Cuerda*(NACA2/10)) PartBody\Sketch.1\Offset.numerocota\Offset = (Cuerda/10+((NACA1/100)*(x/1000*Cuerda)*(2*(NACA2/10)-x/1000)/(NACA2/10)**2) -NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*cos(atan(2*(NACA1/100)/((NACA2/10)**2)*((NACA2/10)-(x/1000)))))/1000 else PartBody\Sketch.1\Offset.numerocota\Offset = (Cuerda/10+((NACA1/100)*(Cuerda-x/1000*Cuerda)*(-2*(NACA2/10)+x/1000+1)/(1-NACA2/10)**2) -NACA34/(100*0.2)*Cuerda*(0.2969*sqrt(x/1000)-0.1260*(x/1000)-0.3516*(x/1000)**2+0.2843*(x/1000)**3-0.1015*(x/1000)**4)*cos(atan(2*(NACA1/100)/((1-(NACA2/10))**2)*((NACA2/10)-(x/1000)))))/1000'; 
%dejamos únicamente lo que será constante y ponemos nombre representativos
%a lo que va a cambiar, en este caso el numero de la cota (numerocota) y la
%coordenada x del perfil (x) esta coordenada va de 0 a 1000, ya que hemos
%decidido que así sea, pero podría estar en cualquier intervalo (de 0 a 1
%por ejemplo).

% Arreglar caracteres conflictivos con fprintf
f2 = strrep(f1, '\', '\\'); % función que sustituye en f1 los \ por \\ ya 
%que la función fprintf usa el comando \ para establecer propiedades como 
%el salto de linea (\n)

% Sustituir en la cadena los elementos que van a variar por %1.0f (ya que
% ninguno de los valores va a tener decimales) que es lo que la función 
% fprintf cambiará por lo que sea necesario.
t = '%1.0f';
f3 = strrep(f2, 'numerocota', t);
f4 = strrep(f3, 'x', t);

% valores por los que deben ser sustituidos los valores por defecto
% anteriores

% el valor numerocota queremos que vaya desde 151 hasta el 180 que es la 
%numeración de las cotas en longitud de los puntos del extradós que hemos 
%colocado en catia

% el valor de las coordenadas en x hemos decidido que tenga una division de
% 20 en los 5 primeros puntos y 36 en los 25 siguiente, para que quede 
% mejor definido el borde de ataque, variando x entre 0 y 1000 que era lo 
% que se pretendía

%el orden de la numeración de las cotas debe coincidir con el orden en el
%que se recorre el vector x
numerocota=151:180; x=[20:20:100,136:36:1000];

% Mostrar en pantalla
fprintf('\n')
for k = 1:length(x);
    C = [x(k),numerocota(k),x(k)*ones(1,8),numerocota(k),x(k)*ones(1,8)]; 
    % vector que contiene los carácteres ordenados que se van a sustituir 
    %en la frase. Tiene un orden especial como puede verse primero aparece
    %la x, despues numerocota, a continuación la x 8 veces etc.
    fprintf([ f4 , ' \n' ], C ) %sustituye los %1.0f por los valores de C
end