- ベストアンサー
プログラム言語について
このプログラムの言語が何かわかる方いらっしゃいますか? これはプログラムの一部です。 自分なりに調べたのですが よくわからずに困っています。 よろしくお願いします。 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 %
- みんなの回答 (2)
- 専門家の回答
質問者が選んだベストアンサー
Matlab ですね。 http://www.mathworks.co.jp/ ただで試したいなら、OctaveやFreeMat 等のクローンソフトを使えばいいでしょう。(完全互換ではない部分もありますが)
その他の回答 (1)
- takepon256
- ベストアンサー率41% (391/936)
よく知りませんが、同じような宣言文があったので貼っておきます。 http://www.obihiro.ac.jp/~suzukim/masuda/octave/man/octave-20060720.pdf