├── LICENCE
└── README.md


/LICENCE:
--------------------------------------------------------------------------------
 1 | MIT License
 2 | 
 3 | Copyright (c) Hazrat Ali 
 4 | 
 5 | Permission is hereby granted, free of charge, to any person obtaining a copy
 6 | of this software and associated documentation files (the "Software"), to deal
 7 | in the Software without restriction, including without limitation the rights
 8 | to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 9 | copies of the Software, and to permit persons to whom the Software is
10 | furnished to do so, subject to the following conditions:
11 | 
12 | The above copyright notice and this permission notice shall be included in all
13 | copies or substantial portions of the Software.
14 | 
15 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
16 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
17 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
18 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
19 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
20 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
21 | SOFTWARE.
22 | 


--------------------------------------------------------------------------------
/README.md:
--------------------------------------------------------------------------------
  1 | # ✈ Hazrat Ali
  2 | 
  3 | # 🚢 Programmer || Software Engineering 
  4 | 
  5 | # Matlab
  6 | 
  7 | # Eular Method : Number 1
  8 | 
  9 | # a=0;
 10 | # b=0.1;
 11 | # n=5;
 12 | # h=(b-a)/n;
 13 | # f = @(y,x) x + y;
 14 | # xi = @(i) a+i*h;
 15 | 
 16 | # x_not=a;
 17 | # y_not=1;
 18 | # while xi(x_not) <=b
 19 |    # y_not+=h*f(xi(x_not),y_not);
 20 | # x_not++;
 21 | # end;
 22 | # printf("%f", y_not);
 23 | 
 24 | # Simpson 3/8 : Number 2
 25 | 
 26 | # clc;
 27 | # clear all;
 28 | # a=0;
 29 | # b=1;
 30 | # n=6; %must be an even
 31 | # delX=(b-a)/n;
 32 | # f=@(x) 1/(1+x^2);
 33 | # xi=@(i) a+i*delX;
 34 | # sum = f(xi(0))+f(xi(n));
 35 | # i=1;
 36 | # while i<n
 37 |     if rem(i,3)==0
 38 |         sum += 2*f(xi(i));
 39 |     else sum+=3*f(xi(i));
 40 |     end;
 41 |     i++;
 42 | # end;
 43 | # sum*=((3*delX)/8);
 44 | # printf("%f",sum);
 45 | 
 46 |  # Simpson 1/3 : Number 3
 47 | 
 48 | # clc;
 49 | # clear all;
 50 | # a=0;
 51 | # b=1;
 52 | # n=6; %must be an even
 53 | # delX=(b-a)/n;
 54 | # f=@(x) 1/(1+x^2);
 55 | # xi=@(i) a+i*delX;
 56 | # sum = f(xi(0))+f(xi(n));
 57 | # i=1;
 58 | # while i<n
 59 |     if rem(i,2)==0
 60 |         sum += 2*f(xi(i));
 61 |     else sum+=4*f(xi(i));
 62 |     end;
 63 |     i++;
 64 | end;
 65 | # sum*=(delX/3);
 66 | # printf("%f",sum);
 67 | 
 68 | # Trapezoidal : 4
 69 | 
 70 | # clc;
 71 | # clear all;
 72 | # f=@(x) 1/(1+x^2);
 73 | # a=0;
 74 | # b=6;
 75 | # h=(b-a)/6;
 76 | # xi=@(i) a+i*h;
 77 | # sum=f(xi(0))+f(xi(6));
 78 | # i=1;
 79 | # while i<6
 80 |     sum+=(2*f(xi(i)));
 81 |     i++;
 82 | # end;
 83 | # sum*=(h/2);
 84 | # printf("%f",sum);
 85 | 
 86 | # Runge Kutta : 5
 87 | 
 88 | # clc;
 89 | # clear all;
 90 | # f=@(x,y)x+y^2;
 91 | # h=0.1;
 92 | # for_x=0.2;
 93 | # init_x=0;
 94 | # init_y=1;
 95 | 
 96 | # while init_x<for_x
 97 |    # k1 = h*f(init_x,init_y);
 98 |    # k2= h* f((init_x)+h/2,(init_y)+k1/2);
 99 |    # k3 = h* f((init_x)+h/2,(init_y)+k2/2);
100 |    # k4 = h*f(init_x+h,init_y+k3);
101 |    # init_y+=(1/6*(k1+2*k2+2*k3+k4));
102 |    # init_x+=h;
103 | # end;
104 | # printf("%f",init_y);
105 | 
106 | # Hazrat Ali
107 | # Software Engineer
108 | # CSE 
109 | # ID : 221010050
110 | # 
111 | 
112 | 
113 | 


--------------------------------------------------------------------------------