プログラム言語について

このプログラムの言語が何かわかる方いらっしゃいますか？
これはプログラムの一部です。
自分なりに調べたのですが
よくわからずに困っています。
よろしくお願いします。


function [U] = UA(XI,ET,Q,DISL1,DISL2,DISL3)
% DIMENSION U(12),DU% C***** INPUT OUTPUT COMMON /C0/ALP1,ALP2,ALP3,ALP4,ALP5,SD,CD,SDSD,CDCD,SDCD,S2D,C2D 06360000
% COMMON /C2/XI2,ET2,Q2,R,R2,R3,R5,Y,D,TT,ALX,ALE,X11,Y11,X32,Y32, 06370000
% * EY,EZ,FY,FZ,GY,GZ,HY,HZ 06380000

global ALP1 ALP2 ALP3 ALP4 ALP5 SD CD SDSD CDCD SDCD S2D C2D
global XI2 ET2 Q2 R R2 R3 R5 Y D TT ALX ALE X11 Y11 X32 Y32
global EY EZ FY FZ GY GZ HY HZ
global N_CELL

% DATA F0,F2,PI2/0.D0,2.D0,6.283185307179586D0/ 06390000
F0 = zeros(N_CELL,1,'double');
F2 = ones(N_CELL,1,'double').*2.0;
PI2 = ones(N_CELL,1,'double').*6.283185307179586;

DU = zeros(N_CELL,12,'double');
du1 = zeros(N_CELL,12,'double');
du2 = zeros(N_CELL,12,'double');
du3 = zeros(N_CELL,12,'double');

%C-----
%for I=1:1:12
U(1:N_CELL,1:12)=0.0;
%end
XY=XI.*Y11;
QX=Q .*X11;
QY=Q .*Y11;
% C====================================== 06460000
% C===== STRIKE-SLIP CONTRIBUTION ===== 06470000
% C====================================== 06480000
% if DISL1~=F0
c1 = DISL1 ~= F0;
du1(:,1)= TT./F2 +ALP2.*XI.*QY;
du1(:,2)= ALP2.*Q./R;
du1(:,3)= ALP1.*ALE -ALP2.*Q.*QY;
du1(:,4)=-ALP1.*QY -ALP2.*XI2.*Q.*Y32;
du1(:,5)= -ALP2.*XI.*Q./R3;
du1(:,6)= ALP1.*XY +ALP2.*XI.*Q2.*Y32;
du1(:,7)= ALP1.*XY.*SD +ALP2.*XI.*FY+D./F2.*X11;
du1(:,8)= ALP2.*EY;
du1(:,9)= ALP1.*(CD./R+QY.*SD) -ALP2.*Q.*FY;
du1(:,10)= ALP1.*XY.*CD +ALP2.*XI.*FZ+Y./F2.*X11;
du1(:,11)= ALP2.*EZ;
du1(:,12)=-ALP1.*(SD./R-QY.*CD) -ALP2.*Q.*FZ;
% for I=1:1:12
U(1:N_CELL,1:12)=U(1:N_CELL,1:12)...
+repmat(DISL1./PI2,1,12).*du1(1:N_CELL,1:12)...
.*repmat(c1,1,12);
% end
% end
%

No.2 ベストアンサー 回答者： rabbit_cat 回答日時： 2009/11/03 01:04

Matlab ですね。

http://www.mathworks.co.jp/

ただで試したいなら、OctaveやFreeMat 等のクローンソフトを使えばいいでしょう。（完全互換ではない部分もありますが）

No.1 回答者： takepon256 回答日時： 2009/11/02 20:29

よく知りませんが、同じような宣言文があったので貼っておきます。

http://www.obihiro.ac.jp/~suzukim/masuda/octave/ …