1 function z = get_z(x,u) 2 3 4 5 if (x < 0 ) 6 disp('invalid X') 7 end 8 9 M = 0.02; 10 P = 0.4; 11 T = 0.12; 12 a0 = 0.2969; 13 a1 = -0.126; 14 a2 = -0.3516; 15 a3 = 0.2843; 16 a4 = -0.1015; 17 18 if x <P 19 z_camber = M/P^2*(2*P*x - x^2); 20 else 21 z_camber = (M/(1-P)^2)*(1 - 2*P +2*P*x - x^2); 22 end 23 24 %z_camber = M/P^2*(2*P*x - x^2); 25 z_thickness = (T/0.2)*(a0*x^.5+a1*x+a2*x^2+a3*x^3+a4*x^4); 26 27 if u==1 28 z = z_camber + z_thickness; 29 else 30 z = z_camber - z_thickness; 31 end 32 33 34 35