数学月間の会SGKのURLは,https://sgk2005.org/
数学月間の会SGKのURLは,https://sgk2005.org/
% 2021-8-13 Ver 1.00
% 変数のリストを表示する
who
% 名前付き変数のクリア
% 3 つの変数、a、b、c を定義し、その後、a と c をクリアする。
a = 1;
b = 2;
c = 3;
who
clearvars a c
who
% whos 変数リストの詳細表示
a = 1;
b = 2;
c = 3;
who
whos
% clear すべての変数を削除する
a = 1;
b = 2;
c = 3;
who
clear
who
% 変数 a b c のうち、b を削除し、a と b は残す
a = 1;
b = 2;
c = 3;
who
clear b
who
% dir フォルダー内容の一覧表示
dir
% ls 現在のフォルダー内容の一覧表示
ls
% pwd 現在のフォルダーの識別
pwd
% 2021-8-13 Ver 1.00
% 加減乗除
a = 1;
b = 3;
c1 = a + b
c2 = a - b
c3 = a * b
c4 = a / b
c5 = a \ b
% 可変(任意)精度の演算
c4
disp('digits(10)')
digits(10)
vpa(c4)
% 可変(任意)精度の演算
disp('digits(30)')
digits(30)
vpa(c4)
% 可変(任意)精度の演算
disp('digits(3)')
digits(3)
vpa(c4)
% 階乗
a = 11;
b = a^2
c = a^3
% 平方根
a = 64;
b = sqrt(a)
c = -1;
d = sqrt(c)
% 実数のn乗根
a = 8;
b = nthroot(a,3) % 3乗根
c = -27;
d = nthroot(c,3) % 3乗根
% vpasolve 方程式の数値的な求解 1
% 例:方程式 x + 2 = 23 を数値的に解き、x を求める
% syms はシンボリック変数の作成
syms x
vpasolve(x + 2 == 23,x)
% vpasolve 方程式の数値的な求解 2
% 例:方程式 x^2 + 4 = 3 を数値的に解き、x を求める
syms x
vpasolve(x^2 + 4 == 3,x)
% vpasolve 方程式の数値的な求解 3
% 例:方程式 x^2 + 3 = 13 を数値的に解き、x を求める
syms x
vpasolve(x^2 + 3 == 13,x)
% 2021-8-13 Ver 1.00
% 行の入力
clear
A = [1, 2, 3, 4, 5]
% 行の結合
A = [1, 2, 3, 4, 5]
B = [6, 7, 8, 9, 10]
C = [A, B]
% 列の入力
A = [1; 2; 3; 4; 5]
% 行から列へ
clear
clc
A = [1, 2, 3, 4, 5]
% B = A'
B = A'
% 固有値 1
% 物理と工学で使う行列と固有値 P38
%
% | 0 -1 1 |
% A = | -1 2 1 |
% | 1 1 2 |
%
A = [0 -1 1 ; -1 2 1 ; 1 1 2]
e = eig(A)
% 固有値 2
% 物理と工学で使う行列と固有値 P38
%
% | 0 -1 1 |
% A = | -1 2 1 |
% | 1 1 2 |
%
A = [0 -1 1 ; -1 2 1 ; 1 1 2]
e = eig(A, 'matrix')
% 2021-8-13 Ver 1.00
% sin(x) の微分 高校数Ⅲ
syms x
f = sin(x)
diff(f)
% cos(x) の微分 高校数Ⅲ
syms x
f = cos(x)
diff(f)
% tan(x) の微分 高校数Ⅲ
syms x
f = tan(x)
diff(f)
% y = 5*e^x * sin(x) のx= 4.7におけるyの微分
syms x
y = 5*exp(x)*sin(x)
g = diff(y)
vpa(subs(g,x,4.7))
ezplot(g)
% 台形則を用いた数値積分 1
d = pi/10000;
x = 0:d:pi;
y = sin(x);
a = d*trapz(y)
plot(x,y)
% 台形則を用いた数値積分 2
d = pi/5;
x = 0:d:pi;
y = sin(x);
a = d*trapz(y)
plot(x,y)
% 2021-8-13 Ver 1.00
% 円周率
fplot(@(x) sin(x))
% 円周率
pi
% 円周率
disp('digits(30)')
digits(30)
vpa(pi)
% 虚数
i
j
a = i
b = a^2
% ネイピア数
exp(1)
% オイラーの恒等式
% e^(π×i) = -1
exp(pi * i)