N = getenv('環境変数名')

subject = strcat(['[', getenv('COMPUTERNAME'), ']', 'テストメール']);

diary2　…　diary('on')と同じ
diary new　…　以後の出力をdiaryファイルに書き込みます。(ファイルが存在しても新規書き込み)
diary on　…　以後の出力をdiaryファイルに書き込みます。(ファイルが存在すれば追加書き込み)
diary off　…　書き込みを中止します。

diary stat　…　diaryモードのON/OFFを表示します。
diary fname　…　現在指定されているdiaryファイル名を表示します。

diary('ファイル名')　…　diaryファイル名を指定して書き込みを開始します。
diary('ファイル名','on')　…　diaryファイル名を指定し、追加書き込みを行うことができます。
diary('ファイル名','new')　…　diaryファイル名を指定し、新規書き込みを行うことができます。

※diary2('file')は、diary2('file','on')と同じです。

>> diary2('diary.txt','new')
>> diary2 fname
現在のdiaryのファイル名 ... [ diary.txt ]
>> diary2 stat
現在のdiaryの状態 ... [ on ]
>> diary2 off
>> diary2 stat
現在のdiaryの状態 ... [ off ]

function diary2(varargin) 

%% 引数のチェックと処理
error(nargchk(0,2,nargin));

opt = 0;
diary_path = get(0, 'DiaryFile');
mode = {'on','new','off','stat','fname'};

for ii = 1:nargin 
    in = varargin{ii};

    if ischar(in) 
        hit = find(strcmpi(lower(in), mode));
        if isempty(hit) 
            if ~strcmpi(diary_path,in) 
                diary2 off;
                diary_path = in;
            end;
        else 
            opt = hit(1);
        end;
    else 
        error('引数は文字列である必要があります');
    end;
end;

if opt==0 
    if strcmpi(get(0,'Diary'),'on') 
        opt=3;
    else 
        opt=1;
    end;
end;
%% 


%% 実際の処理
switch opt 
    % 書き込み開始 
    case {1,2} 

        % diaryファイルへの書き込み開始メッセージ書き込み 
        if strcmpi(get(0,'Diary'),'off') 
            if opt==1 
                perm = 'a+';
            else 
                perm = 'w+';
            end;
            fid = fopen(diary_path,perm);
            fprintf(fid, '%s¥n', '%%% =================  diary 書き込み開始 ==================== ');
            fprintf(fid, '%s%s¥n¥n', '+ 現在時間 : ', datestr(now,31));
            fclose(fid);
        end;

        % diary on
        diary(diary_path);
    % 書き込み停止 
    case 3
        % diaryファイルへの書き込み終了メッセージ書き込み 
        if strcmpi(get(0,'Diary'),'on') 
            diary off;

            fid = fopen(diary_path,'r+');

            body = {};
            frewind(fid);
            while 1 
                tline = fgetl(fid);
                if ~ischar(tline), break, end;
                body{length(body)+1} = tline;
            end;

            body{end} = [];

            frewind(fid);
            for ii=1:length(body)
                fprintf(fid, '%s¥n', body{ii});
            end;
            fprintf(fid, '%s%s¥n', '+ 現在時間 : ', datestr(now,31));
            fprintf(fid, '%s¥n', '%% ========== / diary 書き込み停止 ==========');

            fclose(fid);
        end;
    % ステータス 
    case 4 
        disp(sprintf('現在のdiaryの状態 ... [ %s ]', get(0,'Diary')));
    % ファイル名表示 
    case 5 
        disp(sprintf('現在のdiaryのファイル名 ... [ %s ]', get(0,'DiaryFile')));
    otherwise
        help diary2;
end

2月最後の授業で言ったように､Matlabのoptimization toolboxを持っている場合､この中の関数fsolveを用いることによって非線形の連立方程式の解を求めることができます。しかし、その後調べたところ､数値計算に関する教科書を著しているKenn Juddは彼のホームページでこの関数はイマイチ（”poor”）であると述べ(リンク先参照)、それよりもChristopher Simsの書いた関数コードであるcsolve.mの利用を推奨しています。…
…
参考：コメント：Matlab関数fsolveについて

A Matlab Nonlinear Equation Solver The nonlinear equation solver in Matlab, fsolve, is a poor one based on minimizing the sum of squares of the functions. Chris Sims offers an improvement; see his ftp page.…
…
参考：Chapter 5 public code

[x,rc] = csolve(FUN,x0,gradfun,crit,itmax)
[x,rc] = csolve(FUN,x0,gradfun,crit,itmax,varargin)

% [引数]
FUN        :    解かれる非線形システム方程式。
            FUNはベクトルxを引数として計算した非線形方程式のベクトルFを返します。
            ベクトルFはベクトルxと常に同じサイズである必要があります。
            関数FUNは、関数ハンドルを使って指定することができます。
x0        :    解の探索のための初期値．
gradfun    :    ヤコビ行列の関数．
crit    :    FUNが返す絶対値の合計がこの値より小さいならば、解が収束したと判断されます。
itmax    :    許可する繰り返しの最大回数．もし繰り返しの回数がこの値に達した場合、第2出力rcに4が返ります。
varargin:    関数FUNとgradfunに渡す、ベクトルx以外のパラメータを指定します。

% 方程式の関数
function F = funcs(x) 
    F = zeros(numel(x),1);
    F(1) = x(1)^2+x(2)^2-4;
    F(2) = x(1)*x(2)-x(3);
    F(3) = x(1)*x(3)-1;
    F = F(:);

% 関数funcsのヤコビ行列
function J = funcsj(x) 
    J = zeros(numel(x));
    J(1,1) = 2*x(1);
    J(2,1) = x(2);
    J(3,1) = x(3);

    J(1,2) = 2*x(2);
    J(2,2) = x(1);
    J(3,2) = 0;

    J(1,3) = 0;
    J(2,3) = -1;
    J(3,3) = x(1);

>> [x,exitflag] = mcsolve(@funcs, ones(3,1), @funcsj, 1e-12, 100, 4, 1)
itct 1, af 1.0208, lambda 0.36, rc 0
   x        0.64       1.72       1.36 
   f      -0.632    -0.2592    -0.1296 
itct 2, af 0.0452359, lambda 1, rc 0
   x    0.735407    1.86822    1.35976 
   f   0.0310718  0.0141413 -2.28221e-005 
itct 3, af 8.75519e-005, lambda 1, rc 0
   x    0.733074    1.86082     1.3641 
   f  6.0167e-005 1.72542e-005 -1.01306e-005 
itct 4, af 3.64242e-010, lambda 1, rc 0
   x    0.733077    1.86081    1.36411 
   f  2.9954e-010 -4.20057e-011 2.26965e-011 
itct 5, af 8.88178e-016, lambda 1, rc 0
   x    0.733077    1.86081    1.36411 
   f  8.88178e-016          0          0 
x =
    0.7331
    1.8608
    1.3641


exitflag =
     0

% SMTPサーバの設定
setpref( 'Internet', 'SMTP_Server', 'SMTPサーバ');
% SMTPサーバでのユーザ名（必要なら）
setpref( 'Internet', 'SMTP_Username', 'ユーザ名');
% SMTPサーバでのパスワード（必要なら）
setpref( 'Internet', 'SMTP_Password', 'パスワード');
% 自分のアドレス（Fromになる）
setpref('Internet', 'E_mail', '自分のアドレス');

sendmail('相手のメールアドレス','題名','メッセージ')

sendmail('相手のメールアドレス','題名','メッセージ','添付ファイルのパス')

kk = curvature_numeric(ppp)　…　点列pppから各点における曲率kkを求めます。
kk = curvature_numeric(ppp,dim)　…　点ベクトルの形式をdimで指定して，点列pppから各点における曲率kkを求めます。

※点列pppは2次元，もしくは3次元である必要があります。
※第2引数dimは，どの次元でひとつの点を表現するかの整数値です。
※2次元の場合は曲率は正負の符号付で返ります。

rr = curvatureR_numeric(ppp)　…　点列pppから各点における曲率半径rrを求めます。
rr = curvatureR_numeric(ppp,dim)　…　点ベクトルの形式をdimで指定して，点列pppから各点における曲率半径rrを求めます。

※点列pppは2次元，もしくは3次元である必要があります。
※第2引数dimは，どの次元でひとつの点を表現するかの整数値です。

nvec = normal_numeric(ppp)　…　点列pppから各点における主法線nvecを求めます。
nvec = normal_numeric(ppp,dim)　…　点ベクトルの形式をdimで指定して，点列pppから各点における主法線nvecを求めます。

※点列pppは2次元，もしくは3次元である必要があります。
※第2引数dimは，どの次元でひとつの点を表現するかの整数値です。

% 3次元の例　：　螺旋曲線の接線（赤）、主法線（青）、従法線（緑）
divn = 200;
theta = linspace(0, 5*pi, divn);
theta2 = linspace(0, pi/3, divn);

ppp = 2*[cos(theta).*cos(theta2); sin(theta).*cos(theta2); sin(theta2)];
tvec = tangent_numeric(ppp);
nvec = normal_numeric(ppp);
bvec = binormal_numeric(ppp);

figure( 'InvertHardcopy', 'off', 'Color', [1 1 1], 'Position', [2 300 400 300]);
xlabel({'x'}, 'FontName','Times New Roman');
ylabel({'y'}, 'FontName','Times New Roman');
zlabel({'z'}, 'FontName','Times New Roman');
hold on;
nn=140;
plot3(ppp(1,:),ppp(2,:),ppp(3,:));
axis equal;
arrow3(ppp(:,nn)',[ppp(:,nn)+tvec(:,nn)]', 'r-', 0.1, 0.2, 0.1);
arrow3(ppp(:,nn)',[ppp(:,nn)+nvec(:,nn)]', 'b-', 0.1, 0.2, 0.1);
arrow3(ppp(:,nn)',[ppp(:,nn)+bvec(:,nn)]', 'g-', 0.1, 0.2, 0.1);
view(-18,26);

bvec = binormal_numeric(ppp)　…　点列pppから各点における従法線bvecを求めます。
bvec = binormal_numeric(ppp,dim)　…　点ベクトルの形式をdimで指定して，点列pppから各点における従法線bvecを求めます。

※点列pppは2次元，もしくは3次元である必要があります。
※第2引数dimは，どの次元でひとつの点を表現するかの整数値です。

>> ppp=[[0;0],[1;1],[3;2]]
ppp =
     0     1     3
     0     1     2

>> binormal_numeric(ppp)
ans =
     0     0     0
     0     0     0
     1     1     1

>> ppp=[[0,0];[1,1];[3,2]]
ppp =
     0     0
     1     1
     3     2

>> binormal_numeric(ppp,2)
ans =
     0     0     1
     0     0     1
     0     0     1

% 3次元の例　：　螺旋曲線の接線（赤）と従法線（緑）
divn = 200;
theta = linspace(0, 5*pi, divn);
theta2 = linspace(0, pi/3, divn);

ppp = 2*[cos(theta).*cos(theta2); sin(theta).*cos(theta2); sin(theta2)];
tvec = tangent_numeric(ppp);
bvec = binormal_numeric(ppp);

figure( 'InvertHardcopy', 'off', 'Color', [1 1 1], 'Position', [2 300 400 300]);
xlabel({'x'}, 'FontName','Times New Roman');
ylabel({'y'}, 'FontName','Times New Roman');
zlabel({'z'}, 'FontName','Times New Roman');
hold on;
nn=140;
plot3(ppp(1,:),ppp(2,:),ppp(3,:));
axis equal;
arrow3(ppp(:,nn)',[ppp(:,nn)+tvec(:,nn)]', 'r-', 0.1, 0.2, 0.1);
arrow3(ppp(:,nn)',[ppp(:,nn)+bvec(:,nn)]', 'g-', 0.1, 0.2, 0.1);
view(-18,26);

XYZ = bezierSurf(UU, VV, CPS)　…　CPSを制御メッシュとするN次ベジエ曲面のパラメータ(U,V)における座標を返す．
XYZ = bezierSurf(UU, VV, CPSX, CPSY, CPSZ) 　…　制御点のx座標をCPSX，y座標をCPSY，z座標をCPSZとした制御メッシュのN次ベジエ曲面のパラメータ(U,V)における座標を返す．
[X,Y,Z] = bezierSurf(UU, VV, CPS)　…　CPSを制御メッシュとするN次ベジエ曲面のパラメータ(U,V)における座標x,y,zごとに返す．

パラメータUU, VVは，0以上1以下の実数配列で与える必要がある．
制御メッシュ点列CPSは，1つの3次元の実数行列で与える必要がある．
制御メッシュ点列CPSX, CPSY, CPSZは，3つの同じ次元を持つ2次元の実数行列で与える必要がある．

>> [x,y] = meshgrid([0:5]);
>> z = zeros(size(x));
>> bezierSurf(1/2,1/2,x,y,z)

ans(:,:,1) =
    2.5000

ans(:,:,2) =
    2.5000

ans(:,:,3) =
     0

>> [x,y] = meshgrid([0:3]);
z = rand(size(x));
cps = zeros(4,4,3);
cps(:,:,1) = x;
cps(:,:,2) = y;
cps(:,:,3) = z;
bezierSurf([0, 0.3, 1], [0,0,0], cps)

ans(:,:,1) =
         0    0.9000    3.0000


ans(:,:,2) =
     0     0     0


ans(:,:,3) =
    0.2769    0.4894    0.1869

>> [bx,by,bz] = bezierSurf([0, 0.3, 1], [0,0,0], cps)

bx =
         0    0.9000    3.0000

by =
     0     0     0

bz =
    0.2769    0.4894    0.1869

% 制御点
nn=6;
pts = [[0,1,2,3,4,5,0,0,0,0,0,0,0.97358,1.6353,0.70145,0.41548,0.45184,0.86041];
[0,1,2,3,4,5,1,1,1,1,1,1,0.87172,1.5897,1.878,0.60249,0.34142,0.36963];
[0,1,2,3,4,5,2,2,2,2,2,2,0.89357,1.2886,1.7519,0.94185,0.45533,1.8098];
[0,1,2,3,4,5,3,3,3,3,3,3,0.6127,0.75722,1.1003,0.46098,0.8714,1.9595];
[0,1,2,3,4,5,4,4,4,4,4,4,1.017,1.6232,1.245,1.6886,0.6222,0.87774];
[0,1,2,3,4,5,5,5,5,5,5,5,1.0215,1.0657,1.1741,0.38953,1.8468,0.22224]];

cps = zeros(nn,nn,3);
cps(:,:,1) = pts(:,1:6);
cps(:,:,2) = pts(:,7:12);
cps(:,:,3) = pts(:,13:18);

% グラフィック表示
figure( 'Name','Bezier Surface', 'Color', [1,1,1],  'Position', [100,200,600,300]);

hold on;
[uu,vv] = meshgrid(linspace(0,1,30));
bez1 = bezierSurf(uu,vv,cps);
ax1 = subplot(1,2,1), surf(bez1(:,:,1),bez1(:,:,2),bez1(:,:,3), 'FaceColor', 'interp', 'EdgeColor', 'none', 'FaceLighting','phong');
set(ax1, 'Position', [0.05, 0.05, 0.45, 0.9]);
axis equal;
view(-17,40);
title( 'Bezier Surface', 'FontSize', 12, 'FontWeight', 'bold', 'Color', [0,0,0], 'FontName', 'Arial');

ax2 = subplot(1,2,2), mesh(cps(:,:,1), cps(:,:,2), cps(:,:,3), 1.6*ones(nn), 'Marker', 'o', 'FaceColor', 'none');
set(ax2, 'Position', [0.5, 0.05, 0.45, 0.9]);
axis equal;
view(-17,40);
title( 'Control Mesh', 'FontSize', 12, 'FontWeight', 'bold', 'Color', [0,0,0], 'FontName', 'Arial');

caxis([0,2]);

directory_path = uigetdir('表示されたときに選択されているディレクトリ','ダイアログのコメント')

※directory_path　…　選択したディレクトリのパス（失敗もしくは未選択、キャンセルなら0が返ります）

>> dpath = uigetdir
dpath =
C:¥sach1o¥MATLAB

>> dpath = uigetdir('C:¥sach1o¥MATLAB')

>> dpath = uigetdir('','プロジェクトフォルダの指定')

dpath = uigetdir('','プロジェクトフォルダの指定');

if dpath==0 
    error('プロジェクトフォルダが選択されていません。');
else 
    A = rand(3);
    B = rand(3);
    save(fullfile(dpath,'dataA.mat'), 'A');
    save(fullfile(dpath,'dataB.mat'), 'B');
    save(fullfile(dpath,'dataAll.mat'));
end;

>> dir(dpath)
.            ..           dataA.mat    dataAll.mat  dataB.mat