MATLAB

() * !" # $% #& ' (matrix laboratory) MATLAB -

. + , - !" ) . /( ( $
. MathWorks 1 + % , -
= @ A 1x1: $ ! ; ) <; 9=> 9 8$ 1 2 345 1 6 % 78 -
. ,9
: + , % )! = 1 ) (8B % -
: $ & ) .& o

9- /8= D# % + E & 9> F /B 9 ) () * o

G) & G) & H) $ 6 $ o

+ . & Java 2Fortran 2C++ A 9 ) / )% ) I $ *) o

(Graphical User Interface) 9J @ : =; 9 ) #$ o

K(=4 /= # G /B G= ) GB L , o

MATLAB

. ;J 6 !@ @ ; % ) ) 9<(=4 J G= (Toolbox) ) GB MATLAB

: ;= ;) = L) * % N ) MATLAB ) GB % 9M )

O () $ 9J @ P) !" Q(= J % 9M ) Q(= 9 B 97 SIMULINK : B$

. + <= & 9 = L) * 2+ , 1 G$ R E %
 : + 5 % ;#% ! = MATLAB ) GB % 9( #

[Link]

MATLAB

! % )& , M ) M = ! = <N 1 . 5 B <N 1 MATLAB = :Command Window -

. S B / " Enter 8
. , 5 + MT m file N ) Function & Script /( ) <N 1 :Editor -
. !5 B U> !" ->U> = <N 1 :Current Folder -
 # P) <N 1 + , B & + , =, Command Window <N 9$ = :Command History -
. 5 + !5
. !5 ;= B A # MATLAB 9 V= % 9= # :Workspace -
. R &! = Workspace B V= <N 1 :Variable Editor -

MATLAB Help

. Help MATLAB

= 1=, ) !" ) 9 = + & + ,6 P=4 8; 2MATLAB Function ) -1

<N ! OE ) clc = ) XY . , Command Window Function 6 1=, !" U ) & help
: % P) 2 5 + <= Command Window
& ,)9 [ % X Help MATLAB 2Function ) P=4 8; ) +&X -2
Function & MATLAB Toolbox 1 ;\8 2Help 1 . B !" Function ) =( -$
. + ,+ ]-$ K(=4 /# Y + 8 ) /" ) I )

Demo ! ; $ + , /$ _ ( N ) 1 J(8 % 9`G) ) I ) XF -3

. ;= 9 = L) * MATLAB
 MATLAB ! # ! "

MATLAB $%# &

. , M !" % =Ja b MATLAB Y 1 =Ja :eps

(3.141593) 9E b :pi

+ 5 345 ) $ ; D ) b H) $ 9B& M ! @ " :ans

/ 9) b :inf

. 5 + @) b 1 , ) 6/ f = J %& Not a Number K<4 :NaN

 (Script) ' (% )% *%

Command Window 6 5E * 9$ = 1=, % 9;= /( ) 2 , ! ) b) 78

. D * + <= & 5 + MT m file Q# * 9;= /( .+ , L + J S B & 2L ,S B

:U Y

:Workspace ` V= & Tamrin01 L B % :E Workspace + ;E

 (Function) - *%

29;= /( ; 8 . 5 (generic) 9 8 P) a D * + <= ) MATLAB Function

9 8 2 L Script ) Function 9 & <$ . D * + <= & 5 m file Function

) % 2 F ) m file U& 7 2H) $ 1=, ). 9B& M & && () * J $ P) /" ! )

: ,+ *

function [argout1, argout2, ...] = function_name (argin1, argin2, ...)

. 5 + R 8 9)" g ) M P) & MATLAB ( 8( % function

:U Y

) B % G) & ;= 9( 2H) $ LM V= . , B $ ( ` ) Workspace LM V= )2 ) B % G)

. B& ` 2 + , K G$ Command Window + ;E 9 V= b & B& ` /" % h
. Tamrin02 L LM F) M /= (;) U @ & 7b D(8 2d & c V= R 8 (
. ,) &&! @"& G) 8 ) ) ) U& 9B& M . , ) =, 9B& M & &&! @"& ;) H) $ :.%
. && & V= ! $ ) U& V= ) , ) ) ) 9B& M & , ) & & & iB 8 ) ) ) & 9B& M

R = 82.04 ) ; D= B 25 & < 8$ 5 98 > B U" % @ ! * % + <= ) :.%

. ; 98 > B 2 4# V$ F ) .(atm.cm3/mol/K

V$ F ). ; D= B 78.32 U $ 4) 5 2U $ ) ! = " 7) % + <= ) :1 8$

. ; 4) 5 2 4#

1642.89
= 8.20417 −
230.3 +
 012 # 3 4

. , I,l B & I , (8B L , 9F , = . 5 345 if ) MATLAB I, )I) =

. 5 B2 * if % G) M (N X) I ,l B2 ,) * ) I , (8B J $ N

|| && == >= <= > < ~= :% ;$ 2I , (8B ) I ) D(8

 )& ; y b 2 % + ,+ +& ) B $ )& ; &&% x V= b ;) 9G) $ :1 8$
. ; U 7M V E 2 , B& + , + +& x b J$ N . ; U 9B& M

−10 < ≤ −1 =

−1 < <0 =√

=0 =0

>0 = exp

Otherwise, display this error message: Please enter a real number greater than -10.

.while b(> & for b(> % ;$ MATLAB b(>

. ;J E + 8, [ 8, /= $ b(> LM = S B & 5 K G$ + 8, for b(>

:U Y

. 5 = E b(> LM = 2 * )I ,1 J % $& 5 K G$ b(> )I , while b(>

. ; U 9B& M ) & l - 6 ) !" $ <N 1 ) 9G F & D) &&! ; ) ;) 9G) $ :U Y

. ;) For b(> ) m ) :1 8$
 MATLAB 5 %

! ; ) <; 9 8$ 1 ) ;) . 5 = @ A : $ ! ; ) MATLAB V= 9 8$ 2 , ! ) b) 78
. 5 = @ A 1x1: $

:U Y

(N M "," U @ & ) 2 7 B ". 5 + * ) 1) : $ " 2: $ )

. 5 + <= Enter ! % ";" U @ & 7b % 2 B 7 )& 5 B6 %

:U Y

. + <= linspace = % & ":" 1# % ! = 2 ) ) (N ) " )

:U Y

. ;, n 4= & 9 = L) * = E % + <= ): $ "

 :U Y

. ;J ! " / $ 1 + , K G$ MATLAB 9G) $ 2o M / $ % )

:U Y

. D 9 $ P) MATLAB 9- U8 ( 2 , ! ) b) 78

:U Y
 . D " ) " N ) 9- U 8 2(^ / \ * + -)9- D(8 % L * "." 7b ! * )

:U Y

. D T !" % 9 /# Y 8 ! ; ) + , K G$ MATLAB / $ ) $ G= function

:U Y
 . , ) 10 $ 1 % L , !" 7 ; 10 10 : $ :1 8$

) & ; + MT : $ " % =Ja 9G F & ; &&! ; ) Y ;) 9G) $ :1 8$

)&l > !" " l `(N > & 8 & ; 9B& M : $ 1 ;) D 9G) $ . ; U H) $ % n M
. ; U 9B& M
MATLAB 6% - ! 7

:9$ Y(Y H) $

:9 8 H) $
 :! @ H) $

: = @ - H) $

: (=4 H) $
 MATLAB 84 4 8 9 Debug

Q(. MATLAB .runtime error & syntax error % 6 9 ) B r J, ! H 9;G debug qX7N
. 5 ! 8 command window + ;E error V E & 5 345 ! , L s ! % syntax )I ) r J,

:U Y

: + 5 ) B % L * error 1 ! = 6 Editor <N B X & ) cursor ! ) ) = #

L) * 9@ ) r J, 1 . 5 8 N 6 9 7M V E & 5 8 345 ! , L s syntax )I) r J, % 9`G)

+ , =, 2 A; 1 ). % + , =, ) M ) M&S B )S B9 ) 9 @ & ;= 9 ;,
break % J; ) .6 ; 9 ) M ) M G) ) "% &6; B 6 !" ! ) ] N ) ! ; 8F ) % 97M $ )
. + <= ! = keyboard = % point
 :U Y

B ) (runtime error) ) B! % 9$r J, 2 5 8 345 ! , L s ! % syntax )I ) r J, ; 8

6 !" N ) ! ; 8F )% 9 B $ )& L8 )L* #> ; & ;= 9 ;, L) * 9@ ) ;"9
. M ) M9 ) G) & B )

:U Y

( s (Warning) 9 7M ) I ) X 1 . 5 + 5 @ 9D 9 X 2Editor H) $ 1=, H*

9$r J, & ; 6 ) = /" ) B $ ;= 9 /=; #. 7M 1 . ;J MATLAB

. + 5 L) * 7M !" = 2 X 1 & ) cursor ! * ). ; ) B1>

 MATLAB ) -

1 . . V= ) 9- /$ (8 () * MATLAB 2 , ! ) MATLAB 9 G 78
. 5 =M ;, MTALAB # 8 H) $ ! ; $ J

. 5 sym = $ # 8 tX V= K G$

:U Y

. B = $ & )5 = syms = % ! = 2 5 A>X 78

. . 9( ( $ P) MATLAB 9- U8 % )! = 2 # 8 tX V= $

:U Y

. 9 b subs = )! = # 8 tX V=
 :U Y

V= ) I ) = $ H) $ 1 . D) A % H) $ :.%

. ; 4# 7b ;a H) $ b . ; # 8

, ," = + . − exp √"

3
 11 @ 82 # : ; < => ?

. 9( ( $ P) 9 @ @ 9- /$ (8 ! = # 8 tX V= % + <= )

u=5 . D * + <= # 8 tX V= ) 9- % @ u=5 ) diff = : =>?

. , K(=4 Q$ = @

:U Y

7b b=5 b & ; H) $ % ) 9( ( $ P) + ,! ; b=5 . D) A % H) $ :.%

. ! 5 y=4 & x=2

, = tanh + coth

? ? ? ?
)* )- * )- * )- *
)+ )+- ). - ).)+
U D= . D * + <= # 8 tX V= ) 9- % @ U D= ) int = : =: ;
. + <= @ ;a U D= ) = 1 % ! = 1 ;\8 . , 1 G &1 G P) = @

:U Y
y ) =) i H) $ 1 G @& U D= . &" )x&y ) L* 8 * + , K G$ H) $ 1 G U D= :.%
. ; 4 $1% x )&x $1%

8 & 8 > ) * = 1 . , H) $ % 7b 9D 8 > ) I ) limit = :82

. 7b _a

:U Y
 8 # A

M ) 9r) < % 1 )+ ,6 8 . G) & G) & 8 6 ) * MATLAB

. ) ) 8 6 9<(=4 J 2 =M 1 8 D5 & 9 ) & ;=

L8 l (7 I b b &> / $ & ) = 1 . 5 plot = $ MATLAB G) & 8 6

. ;J

:U Y
 . 6 LJ, = 1 ) 8 % =5 ) ! =

:U Y

. ; GB 8 6 = 1 J ) 9) = ) = 1 =5 ) v B ) plot = % help 1= @ )

! = l 4= plot = ! = 9 / & /D . ; 6 + 4# I b G$ ) x=[0,5] +% ) % H) $ :1 8$

. ;

f(x)=exp(x) g(x)=x+1

; @ Tools ; %2 8 D5 & ) 9 = ). + ,6 8 V$ ! = 2 8 D5 & $

% LJ, )) D5 & ) I ) 8 & ) ( 2 D [& . E9 = D5 & ) & l 4= ! = Edit Plot

.( ; B$

. 5 ! 8 8 % D5 & ) I ) ; 2 D5 & ) I ) 8 & )! ( )

 8 8; 2y & x & ) ! ; 2LJ, )! ; ! *2 + ,6 ) ) 8 6 D J
. 9 = L) * insert ; % J 1 9 8$ . , (Colorbar) g GB & (Legend)

+ , + &" % U& B = 1 % 9`G) . + , K G$ 9$ = w ;= 2 , ! ) 8 R & ) 9# 8 9 8$ )

.

title(‘expr’); 8 ! ;
xlabel(‘expr’) x ! ;
ylabel(‘expr’) y ! ;
axis([xmin,xmax,ymin,ymax]) y&x +&
figure; B LJ,
figure(n); n+ 8, LJ, ) 9 =
legend(‘expr1’,’expr2’,...) 5b 8;

+ ,6 8 ) 2MATLB + , K G$ = & 8 6 )I) + ;E LJ, D5 & % + <= ) :1 8$

+ 4# N ) 8 R 8 +& 1 ;\8 . ) * 5b 8; &y ! ; 2x ! ; 2LJ, ! ; 2L * 1 8$
. V$
6 !" G) LJ, ! = & 5 LJ, 6 ) B + ;E 2figure = % + <= ) 2LJ, % =5 ) 6 )
.

. 5 + <= subplot = % 2(B4 )% P) 8 & 6 )

:U Y

9 8 2 6 LJ, ! = + ;E- % 1 % + ;E- % ;a ) LJ, R 8 + ;E 6 b$ )

% ) 2 + ;E R 8 <N ! = 2 / " & ) ! ( ) & ( % LJ, g 9)" L 7= ) + ,6 LJ, R 8 <N
. 6 b$ 4# + ;E

:U Y
 inline -

. + <= H) $ K G$ ) inline H) $ % ! = & m file 1=, ) % 2+ $ H) $ )

:U Y

. 5 + <= symbolic & inline H) $ G) & 6 ) ezplot =

:U Y
!" B LJ, & 5 OE 8 + ;E = 2 5 B ezplot & plot = +@ MATLAB 8 6
hold = % 2+ ;E ! 8 6; 6 B LJ, & 6 ; x<> LJ, < N 9( * = 6 4) J $ N . 5 6
. 5 + <=

:U Y

)I ) v B. + ,! ; % U& B 8 1 )I ) = 6 $ L) * MATLAB D 8
. A>X MATLAB 8; % + <= )! = = %

polar 9 7* P=4
loglog 98= D# P=4
stairs (E 8
bar ( 8
hist 9& 8
 8 A

& 5 6 G) ;= 9 ` 6M U& = . 5 6 b$ & ) 9( ;) 6 b$ 2 G) 8

) z & y 2x ) 2 = 1 . 5 plot3 = $2 G) /8M 8 6 . ;= & &q7 & =
. ;J 6 ) 1 % LN > 6M & D &&! ;

:U Y

. , 2 ,! ) G) & 8 6 8 * \ " ; 8 2LJ, R & ) I ) J )9 =

! ; ) z & y 2x : $ = & 1 . L) * surf & mesh = & $ l (7 I b & ) ]7 & & 6

)& ; = y&x & & = # . ;J 6 2 + + MT " P=4 : $ 1 & & D &&
. 5 L , y&x l (7 I b ;, )

= % 2 @ * + <= & 6 ) :s & ;, ) z ) l (7 I b L , J $ P) y & x / $ )

. *y&x J, Ib P=4 2 ;J l 4= ) 9F b = 1 . 5 + <= meshgrid

:U Y
. ;J U =; & & ) g R 8 + shading = . ;J U =; & " g colormap =

:U Y
 . 5 + <= ezmesh & ezsurf = % 2H) $ G) 8 6 )2 G) & 8 6 ) )5

:U Y

! ; 2Q ; ! ; 2 /(J, % ) .( = E P) 6 & P) 6 ) ; 6 + 4# + & % H) $ :1 8$

. ; ! = % H) $ & ) M &g . * 5b 8; &

, = cos 0
+ exp sin y

= 3
−1 = 0
−1

−1 +1
+ = 10
2 4
 MATLAB 0C' #

b & 1 U G . , ) (false or true) /" b [% 5 mXF 9$ ) 9b7; )&

H) $ & D(8 % G$ $ MATLAB 9b7; )& . ) <N & )1 ) ) ) MATLAB
. 5 ;) 6 b$ (Logical operators) 9b7; D(8 & (Relational operators) 7) D(8 = & ) 5

MATLAB V= 9 8$ a ;;J L8 9 $ & ) MATLAB 9b7; & 7) D(8 & H) $ 9 8$ : B

. 5 y : $

:% ;$ ! , w ;= H) $ & MATLAB 7) D(8 : C ;1 D

eq() H) $ & == : & $

ne() H) $ & ~= :{ b K# 4

gt() & lt() H) $ & > < : =@ ) & =Ja

ge() & le() H) $ & >= <= : & =@ ) & & =Ja

) (1 ) b 2 ,) *) A 7) J$ N N1 ) ! , w ;= H) $ & D(8 1 ) +
. D ) ( <N ) b 2 , *) A 7) J$ N & D
<N ) ) 9 = & 7 + 8, & ; ) % &&! ; ) : $ ;) 9G) $ :1 8$
. ; U 9B& M ) & ; + MT D : $

:% ;$ MATLAB 9b7; H) $ & D(8 : 0C' ;1 D

 and() H) $ & & :9b7; And

or() H) $ & | :9b7; Or

not() H) $ & ~ :9b7; Not

xor() H) $ :(exclusive or) P Xor

L* > J $ N 2 ;J L8 ) & ) any H) $ . ;;J L8 & ) 9b7; H) $ S B all & any H) $ &
J $ N all H) $ . D ) ( <N ) b r & D ) (1 ) b 2 , ) <N . ) " % 9J
. D ) ( <N ) b r& D ) (1 )] N b , ) <N . ) " 9 8$

<N b : $ / = G$ ) & ;;J 9 ) P) / = b ;;J L8 / $ & ) 2H) $ & 1

. ; D )
! -1 ! $ ) ) J & inv = $ J : $ 1 t JG . D) A + 4# 9G) : $ :1 8$
. 5 ! J & 1 = ; 9 ) 9b7; H) $ & D(8 % + <= ). ;
 MATLAB 3 - # 4

MATLAB & 5 K G$ char tX = =, . 5 mXF =, ) MATLAB = =,

. 9 $ P) = /$ (8

:U Y

) )2 ,L , 4) % =5 ) V= J$ N 2= tX V= 9 $ 1=, ) B$ )
+ MT /" 6 4) V= 7 & & 6 , ) =, & L , V= XY @ D ) . , B $ V= "
. ,) ) ) & = G$ ) 26 ;

:U Y

+% V$ A; ) =, /= 9# M B! - 2 ! = /;$ 2 = =, ! +% 6 )
9# M B =, /= 2 M N ) = 1 . - A; 1 ) char = % ! = 1 ;\8 . = =,
. ,: $ B =, +% 6 $ ;J -

:U Y
 : + , + &" % U& B = =, ) I ) H) $ % 9M )

strcat(‘str1’,’str2’,...) \ D J ) = =, ;a &
upper(‘str’) ;J L $ | ) '& > ) = =, '& > 9 8$
num2str(n) = =, ) L $
eval(‘str’) = =, B

. ; U 9B& M ) & ; %: $ & ; &&! ; ) (n XY ) 9G F ;) 9G) $ :1 8$

611 … 618
4=5 9
681 … 688

% & ! @"! ; ) ) + 8 x % + V= $ 4# 9- H) $ &> = =, ;) 9G) $ :1 8$

. ; U 9B& M ) & ; 2 & ! @" + , & b H) $ !" & ; )
 1 8'E

.6 ;J & /" 8 * 1 9- (8B ;a ) u(G= 9N P=M = MATLAB

. 5 q 7 H) $ ) &&! ; ) ) Q# * !" Q - $ (8B ;a K G$ r 8G

:U Y

;a Q - m U Y . ;J 2 5 & & ! @"! ; ) b + , & (8B ;a b 2m =

. + ,+ H) $ ) U& ! @ " ! ; ): = +2 +1 (8B

. ;"9 ) roots = $ (8B ;a 5

% ;$ : =4 0
+ 2x − x + 3 # G 5 :U Y

;a Q - & ;J ) Q# * &&! ; ) (8B ;a 5 poly = roots = :J

. ;J 9B& M ! ; ) 5 !" ) I ) (8B

:% $ 5 & -2 23 !" 5 (8B ;a Q - :U Y

. D deconv & conv = $ (8B ;a ) 6 b$ & l -

 :% $ ; = −3 0
+x& : =2 0
− x + 6x (8B ;a & l - :U Y

:% $ #& (8B ;a & % ) (N > 6 b$ &

) H8B ) 9@ ) MATLAB (8B ;a u <$ & H8B 2 5 K G$ ) $ (8B ;a J; ( )

2 ) u <$ & H8B a ) ! ) +% 6 2 B$ ) =J /;$ . 5 (8B ;a ' G
. ,)9 N ) %r I , )! ) 6

:% $ ; = −3 0
+x&: =2 0
− x + 6x (8B ;a & u <$ & H8B :U Y

= $ D# D= . + , K G$ polyint & polyder = (8B ;a D# D= & u=5 ) MATLAB

. 5 9( ( $ P) polyint

:% $ : =2 0
− x + 6x (8B ;a U D= & u=5 :U Y
 ( ;J &&! ; ) ) 1 G Ib 7b (8B ;a b 2 , + , !" ) r ) polyval =
. ;J

:K(=4 I b : =2 0
− x + 6x (8B ;a b :U Y

2H8B LN > . ; U D= & u=5 & 5 & 2 ; l 4= 3 % =5 ) & & <= B % (8B ;a & :1 8$
. ; K(=4 + 4# I b (8B ;a & b 1 ;\8 . &" ) / " 6 b$ & l - 2u <$

. , ) =, polyder A J(8 ;) 9G) $ :1 8$

(8B ;a Q - 2 = 1 . (8B ;a 3=4 2 , ) 9 9; ; [% ) ) I ) = % polyfit =

y &xIb ' G ) & 2H) $ 1 . ;J G) & <N 7b G$ ! % + ;; M 1 =/) ' G
Q - & ;J &&! ; ) ; Ib 1 % 6 4 9 4# (8B ;a B & ;J ;=
. ;J G) L* > [& % l (7 (8B ;a

:U Y

40

35

30

25

20

15

10

0
0 1 2 3 4 5 6 7 8 9 10
9(b= V= & H) $ ) : handle . 9 t " handle $! = H) $ MATLAB :handle F%
!" 6 u F % + <= H) $ ) 2 &&! ; ) D H) $ )! )% MATLAB H) $ % ).
handle .H) $ !" handle u F % & ;;J E9 = ( , } ) eval = 9 /# Y A) = =, P)
. 5 345 @ D(8 $

:U Y
 6% G % ?%

. 5 fzero & roots = & $ MATLAB r G 9) 5

. @ *9 ) L* (B 2 / (8B ;a ) I ) roots =

. &" ) : =3 <
− 2 x 3 − x + 6x + 1 (8B ;a 5 :U Y

#& I , )% = 1 . + <= L) * (8B ;a . & (8B ;a % 6 B # G @ ) fzero =

. ;J E+ , & #& I , ) 5 1 =J & l (7 # G 5 ! E )

. ,+ = fzero H) $ ) & &! @ "! ; ) ) H) $ 1 handle & ;J L8 H) $ & ) = 1 : B

:U Y
. = fzero )6 MATLAB ) + , K G$ H) $ handle ! = 1 ;\8

:U Y
7M V E r & , 9 J X ) 2+% ) 1 /= & H) $ 1 !" I , = # . , ) +% ) = fzero H) $
. 5 + =

) l (7 8 A;$ &> \ = 2optimset = $ . ,9 optimset = ) fzero = 8 A;$

.(optimset & fzero = 9 8; ) , B) ,9 + fzero ) ! @"! ; )& ,9 =M fzero =
 : &" ) K(=4 9> % r G 5 :1 8$

= 2 − => .

= cos − 2sin

= 1 + logA√1 + xB + logA√1 − xB

G ;

. 6 b$ L) * 97M . r G + D= & 97M r G + D= 8 *& ) MATLAB B r G + D=

:% $ 97M r G + D= 9( LJ, . L> MATLAB 9 $ = $ 9@ )! = 97M r G + D=

4 =C

. , r G + D= 8 &> )B& r / ) x 2Q - : $ A 7) 1

; L> % r G + D= :U Y

2 + − 3" = 2
:D− + 3 + 6" = 0F
− 4 + 5" = −1
H) $ F ) 9- r G ) J; ). L> P) ! = fsolve = ) 97M . B r G + D=
t > ) 2 r G + D= Lb= V= G$ ) 1 ;\8 . U fsolve H) $ ) &&! ; ) !" & , @ B P)
. ,+ fsolve H) $ ) & &! ; ) #&

:U Y
t > V$ ) 1J8 r G L> l B 2 ;, ) =, l B % =5 ) 1J8 97M . r G + D= 9 "%
. ; & <$ #&

. "9 ) & <= /) B 2 & <= #& t > ) m r G + D= ) :U Y

. , M fzero H) $ ) 5 J(8 N1 . ; L8 #G & ) = fzero A 6 fsolve = : B

:U Y
 8 A;$ &> \ = 2 optimoptions = $ . ,9 optimoptions = ) fsolve = 8 A;$
& fzero = 9 8; ) , B) ,9 + fsolve ) ! @"! ; )& ,9 =M fsolve = ) l (7
.(optimoptions

; L> % B r G + D= :1 8$

=G
=−

− =2

sin " = 1

# 8 P) 9- r G & ) + D= 1 . D * + <= 97M . r G + D= L> ) solve H) $

. ;J L8 + , =,

. ; b fsolve H) $ /4 E ) /4 E & ; L> solve = ) u) 1 8$ B r G :1 8$

 8D =: ;

. + , K G$ MATLAB quad & trapz = & 2 D# D= )

[& ) D# D= = 1 . 5 ;= y &x b D ) % 9 " & ) trapz $ D# D=

. b %&T

! = 6 A H) $ 2fzero & fsolve A u) = D ; 8 . 5 H) $ & ) quad $ D# D=

. U quad H) $ ) " handle & , @ BL "! = 6 & U quad H) $ ) 6 b=

:U Y
 # % 4 D( H% G ; (2

2L < r G 9( ( $ L> = # . 9( ( $ P) 6 & [& )6 L < r G L> 9 $ MATLAB

L < r G + D= L> ) G= = . ,) B # G ) 9( ( $ L> oM ) &
. + , + &" % U& B NXM P) + , K G$ MATLAB #& , )

6 -
, 1 E $ 9 @ $ $ g [& ode23

, 9 $ 9 @ $ $ g [& ode45

, stiff Lv L> )9 @ $ 9;8- L> [& ode23s

, b %&T + *9 @ $ 9;8- L> [& ode23t

, b %&T + *9 @ $ 9;8- L> [& ode23tb

, stiff Lv L> ) 9 @ ;a 9;8- L> [& ode15s

, 1=# <5) " 9 @ ;a 9;8- L> [& ode113

@ B P) )L < r G + D= )I) r G 2r ) = $L < r G + D= L> )

% , & #& , 1 ;\8 . , U &&! ; ) r ) H) $ ) 2!" handle & H) $ 6 & , =, D L
. ,U ) r ) H) $ ) 2 , L> !" L < r G 6 4 +% ) ) + 8 )

H
= − 1 =1
H

+ 2x + 5y = 0 2 =0 =1
)-+ IJ )+
). - IK ).
 6 ;J L $ U& B # G & ) & B # G =) :L> [&

dy
="
L dx F
dz
= −2xz − 5y
dx

. ; L> 9( ( $ P) %L < r G :1 8$

− 2x + 5y = 0 1 =2 =1 =0
)O+ IJ )+ P )-+ P
). O IK ). ). -

+5 = −1
)+ )G

Q )-G
). ). F 0 =2 " 0 =1
)G R
= −1
− =1
).
). -

6 & ! % ) K(=4 /=A(. V$ + . 5 % /5; & 2 = E L IX=M =" :1 8$

. ;
4+C →T+U U ⇆W+X

Z [ Z [
TYR = 1 TYR = 0.5
[ [

. 5 + <= dsolve = % 2L < r G 9( ( $ L> )

:U Y

. ; b 8 * ) /) B & ; L> 9( ( $ [& ) 8 * + ,q 7 L < r G :1 8$

 ' ' I!

) GB % 9J 2o M 7) & 5 9<(=4 ) GB $ MATLAB (curve fitting) 9; ; [% )

. , 9; ; [% ) 3=4 MATLAB

. + , K G$ polyfit = 2 (8B ;a $ 9; ; [% ) )

:U Y

) ; + 1 % Q ; (8B ;a & r G . + ,+ % U& B 1v # $ & ! ;) 9# G$ + :1 8$

. &"
97M . 9; ; [% ) ) nlinfit = 2 + , K G$ 9; ; [% ) ) MATLAB 9$ = D %
) 26 !" & ) + 9; ; [% ) 6 4 # G 2 = 1 % + <= ). D * + <=
. ,+ = nlinfit H) $ ) H) $ !" handle 6 & , =, H) $

D J ) f = & ; ! = nlinfit = $ 1v # $ & ! ;) 9# G$ + [% ) ) 9<(=4 9- r G :1 8$

. ; b
 #

. 5 G) & 6 & G) 6 MATLAB 9) &

. 5 + <= interp1 = %2 G) 9) & )

:U Y
!5 ) +&X ;J + <= 9) & ) 9<(=4 /,& % interp1 = 2 + ,+ !5 UY 78
. @ * + <= = v5cubic [& 2+ , +

= 1 K(=4 /,& % . ; 9) & 1v # $ ! ;) 9# G$ + % 9 4# b 2interp1 = % + <= ) :1 8$

; + <=

% UY ;J + <= 9) & ) 9<(=4 /,& = 1 . 5 interp2 = $ G) & 9) &

. + ,+ R 8

:U Y
. ; 1 8$ !" & ) G) & 9) & )I ) = & 2 ; l 4= + 4# + V= & 9G) $ :1 8$
 ! % 4 D( H% G ; (2

= & 1 % + <= + . + <= ! = bvp5c & bvp4c = & % 2 % , )L < r G + D= L> )

Command Window = & 1 % + <= 9( N. , & <= / " + , + <= 6= D# D J ; 8
: , % P) MATLAB

>> sol = bvp4c(odefun,bcfun,solinit)

: ;;J &&! @" 2 = & 1

. + , =, !" L < r G 9G) $ U ; 2U& ! @ " -

. v = 1 ) 3=4 =M ) % , 9G) $ U ; 2 & ! @ " -
. 5 bvpinit = $ \ = % V= 2 ! @" -

:# 9 = J

L < r G + D= L> ) 9( < r G &> L ; 8 b * !" =, =M &L 1 +

. 5 + + . & ode23 = )& 5 /$ #& , )

# 9 =J

: D K G$ % 9( P) 2 # % , & > H) $

function res = bvp4bc(ya,yb)

res=[];

# G = )& V= b & > 2(ya) U& &&. 9B& M & &&& H) $ 1 2 D A>X 78
(b % ) 8 % L < # G = )& V= b & > 2(yb) & & & & (a % ) _a 8 % L <
, ) & $ Q# * % , 2 #> 1 . , % , + 8*) & > (res) H) $ 1 9B& M . ;,
. ,+ res ) & $ D 8 & , ) <N ) ) ) & $ !" 8 J 7)

9 =J

l B Lb= V= % b U& ! 8@ " . ! 8@ " & = 1 . 5 bvpinit = $ ! 8@ " 1

Ib = )& V= ) #& t > 2 & ! 8@ " . ;, l (7 2Lb= V= b !" = )& V= ) #
( # G 2 & >) 9( < r G + D= = )& V= @ XY . , H) $ 1 8 U& ! 8@ " Lb= V= % + , 345
: 5 D; solinit V= 2 , ) l (7 10 9# 3 22 21 I b 10 $ 1 1 ) % = ) & V= & % LJ5=

>> solinit = bvpinit ( [1:10] , [1 -1] );

= @ A -1 ) ) ) & = ) & V= ) #& t > & 1 ) ) ) 10 $ 1 I b 9 8$ U& = ) & V= ) #& t > b
. + ,

:U Y

:L>

: 5 =, N1 ) % ,& 5 L $ % P) U& B # G & )2 & B # G

 K% L M ( H% G (2

. 9) = L) * V= % ] N H) $ P) u=5 9$r G 2L < ] N . r G

: , % H) $ P) ] N . r G ! $9 29( P)

\, , : = 0

=, ] N . L < r G + D= J$ N . = ) & V= u=5 yp & = ) & V= y 2Lb= V= t 2 7) 1

. ) ; M !" u=5 & = ) & V= ) P) yp & y 26 , )

: ,9 + <= % P) ode15i = % 2 r G + D= 1 L> )

[t,y] = ode15i(@odefun, tspan, y0, yp0); yp( t=t0 ) = yp0

decic = $ yp0 b . = )& V= u=5 ) #& b yp0 & = ) & V= #& I , y0 2 7) 1

: ;; m N % # G ) yp0 b . L) *

\0, 0, :0 = 0

: % P) yp0 ) decic = % + <= +

[y0mod,yp0mod] = decic(@odefun, t0, y0, fixed_y0, yp0, fixed_yp0)

& fixed_y0 % t > &I, J; & # L> I , D & #& t > yp0 & y0 % 9J 2 7) 1
. ,9 1 G$ fixed_yp0

:1 8$

y*sin(y’)+z*y’=-2x

z*-z’=-y^2+1

y(1)=5, z(1)=-2
 # % 4 1 H% G ; (2

) = )& V= 2 r G 1 . ,9 L> B r G G$ ; L < r G + D= 9 @ 29- Lv

;)6 b$ Lb= V= ) u=5 * = )& V= & Lb= V= ) u=5 = )& V= = &
. ,9

L> (> ! $9 2 ;, ) L) * u=5 = )& V= % ] N P) 2u=5 * = )& V= J$ N

9# 8G L < r G + D= L> [& % 9J ) L < r G L> 2:s & L> B r G =) 2
.

u=5 V= % ] N P) L) * u=5 * V= J 7) 2 ;, ) 9;8- P) B r G 9 %

. , L> ! 8 ) B&L < r G 2 ;,

]P = P , P , … , ^ , "P , … , "_
⋮
]
^ = ^ , P, … , ^ , "P , … , "_

0= P , P , … , ^ , "P , … , "_
⋮
0= _ , P, … , ^ , "P , … , "_

) 2!" ; . + <= ) #& , ) 9# 8G L < r G L> = ! 8 % 2Q(= r G 1 L> )

. , K G$ ; 9 345 !" * & u=5 = )& V= Mass : $

. ,9 odeset = $ Mass : $ K G$

] ^∗P ^∗P
d^∗^ 0^∗_
a. 5 ⋮ 9 = 5 ⋮ 9 a=c e
0_∗^ 0_∗_
"]_∗P _∗P

:1 8$

y'=x*z+y^2, y(1)=-2

z*exp(-z)=y