clc
clear all
H=char(input(' Mode 1 OR 2  '))
switch H
case '1'
disp(' This Mode Is For Lift Calculation ');
V=input(' enter your velocity: ');
P=input(' enter your denisity: ');
CL=input(' enter your cl: ');
S=input(' enter your wing area: ');
L=1/2*V^2*CL*P*S
L
F=char(input(' do you want to figur lift with velovity '));
switch F
    case 'yes'
        z=0;
        m=0
        h=input(' Enter Vleocity Range ')
        z=1/2*CL*P*S;
        m=h.^2;
        q=m*z
        plot(h,q)
        xlabel('Velocity')
        ylabel(' Lift')
    case 'no'
        disp(' OK ')
end
R=char(input('  do you want to know how max weight capasity ? '));
switch R
    case 'yes'
        W=0;
        W=L/9.85;
        W
    case 'no'
        disp('               Thank you for use me        ')
end
case '2'
    disp(' This Mode Is For Wing Area Calculation ')
V=input(' enter your velocity: ');
P=input(' enter your denisity: ');
CL=input(' enter your cl: ');
W=input(' enter your Weghit: ');
S=(W*19.72)/(P*CL*V.^2);
S
Lift=9.86*W
end
disp('     this software for geometrey calculation!    ');
Y=char(input('      Mode zozanaghei OR mostatili       '))
switch Y
    case 'zozanaghei'
        disp(' this methode is for zozanaghe wing ');
        S=input(' enter your area wing:    ');
        AR=input(' enter your aspect ratio:    ');
        T=input( '  enter your Taper ratio:    ');
        b=0;
        b=(AR*S)^0.5;
        A=0;
        A=T+1;
        C1=0;
        C1=S/(A*b);
        C2=0;
        C2=(T*C1);
        C1
        C2
        b
    case 'mostatili'
         disp(' this methode is for mostatil wing ');
        S=input(' enter your area wing:    ');
        AR=input(' enter your aspect ratio:    ');
        C=0;
        C=(S/AR)
        M=0;
        M=C^0.5;
        Chord=M
        b=0;
        b=AR*M;
        b
end
        disp('     Thank You For Use This Software       ')