Few lines of Matlab code need to be converted to Fortran. As you will see the lines are simple statements
## Deliverables
Matlab code
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C! Boundary Conditions - B.C.
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Er=0; % Radial Volt
sigmar=10; % Radial Stress
sigmart=0; % Shear Stress
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C! This is to read properly the knowns or inputs
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
b=cell(1,order);
A=cell(1,order);
TM=cell(1,order);
x=zeros(NoEq,NoEq);
% Applying concentrated-point (delta loads) @ theta0
theta0=pi/4;
%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%c! Fourier Expansion
%c! Eq1, Eq2 or u-r, sigma-r are multiplied by cos
%C! Eq3, Eq4 or u-theta, sigma-r-theta are multiplied by sin
%C!
%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
for n=1:order
b{n}(1,1)=(1/pi)*(-e33)*Er*cos(n*theta0);%Eq1 u-r r=Radial
b{n}(2,1)=(1/pi)*(-e33)*Er*cos(n*theta0);%Eq2 sigma-r
b{n}(3,1)=0; %Eq3= is always zero u-theta
b{n}(4,1)=(1/pi)*e33s*Er*sin(n*theta0); %Eq4 sigma-r-theta
b{n}(5,1)=(1/pi)*sigmar*cos(n*theta0); % Radial Stress Boundary
b{n}(6,1)=(1/pi)*sigmart*sin(n*theta0); % Shear Stress Boundary
end
%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%c! Fourier Expansion
%c! Eq1, Eq2 or u-r, sigma-r are multiplied by cos
%C! Eq3, Eq4 or u-theta, sigma-r-theta are multiplied by sin
%C!
%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
for n=1:order
%C! n can not be n=0 as t=pi/2*n
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C!
%C! Layer 0
%C!
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C Evaluating A unknown coeficients in layer 0
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
n0=0;n1=2;
t=0;% cosines coefficients in radial while sines=0
x(m+1,n0+1)=Eq1_0s(n,t,gm,dyna0,A0,prp0,prpp);
% Equation 2 in sigma-theta 0 layer
x(m+2,n0+1)=Eq2_0s(n,t,gm,dyna0,A0,prp0,prpp);
t=pi/(2*n);% cosines coefficients in radial while sines=0
% Equation 3 in sigma r-theta 0 layer
x(m+3,n0+1)=Eq3_0s(n,t,gm,dyna0,A0,prp0,prpp);
% Equation 4 in Dr 0 layer
x(m+4,n0+1)=Eq4_0s(n,t,gm,dyna0,A0,prp0,prpp);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C Evaluating F coeficients in layer 0
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% Equation 1 in sigma-r 0 layer
t=0;% cosines coefficients in radial while sines=0
x(m+1,n0+2)=Eq1_0s(n,t,gm,dyna0,F0,prp0,prpp);
% Equation 2 in sigma-theta 0 layer
x(m+2,n0+2)=Eq2_0s(n,t,gm,dyna0,F0,prp0,prpp);
% Equation 3 in sigma r-theta 0 layer
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+3,n0+2)=Eq3_0s(n,t,gm,dyna0,F0,prp0,prpp);
% Equation 4 in Dr 0 layer
x(m+4,n0+2)=Eq4_0s(n,t,gm,dyna0,F0,prp0,prpp);
m=0;n0=0;n1=2;
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C! Layer 1
%C!
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C! Physical unknown A
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% Equation 1 in sigma-r general layer
t=0;% cosines coefficients in radial while sines=0
x(m+1,n1+1)=Eq1_1s(n,t,gm,dyna1,A1,prp1,prpp);
% Equation 2 in sigma-theta general layer
x(m+2,n1+1)=Eq2_1s(n,t,gm,dyna1,A1,prp1,prpp);
% Equation 3 in sigma r-theta general layer
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+3,n1+1)=Eq3_1s(n,t,gm,dyna1,A1,prp1,prpp);
% Equation 4 in Dr general layer
x(m+4,n1+1)=Eq4_1s(n,t,gm,dyna1,A1,prp1,prpp);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C Boundary Conditions (known surface traction)
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
t=0;% cosines coefficients in radial while sines=0
x(m+5,n1+1)=sigmarn(n,t,r2,dyna1,A1,prp1);
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+6,n1+1)=sigmartn(n,t,r2,dyna1,A1,prp1);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C! Physical Equations Unknown C
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% Equation 1 in sigma-r general layer
t=0;% cosines coefficients in radial while sines=0
x(m+1,n1+2)=Eq1_1s(n,t,gm,dyna1,C1,prp1,prpp);
% Equation 2 in sigma-theta general layer
x(m+2,n1+2)=Eq2_1s(n,t,gm,dyna1,C1,prp1,prpp);
% Equation 3 in sigma r-theta general layer
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+3,n1+2)=Eq3_1s(n,t,gm,dyna1,C1,prp1,prpp);
% Equation 4 in Dr general layer
x(m+4,n1+2)=Eq4_1s(n,t,gm,dyna1,C1,prp1,prpp);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C Boundary Conditions (known surface traction)
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
t=0;% cosines coefficients in radial while sines=0
x(m+5,n1+2)=sigmarn(n,t,r2,dyna1,C1,prp1);
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+6,n1+2)=sigmartn(n,t,r2,dyna1,C1,prp1);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C! Physical Equations Unknowns F
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% Equation 1 in sigma-r general layer
t=0;% cosines coefficients in radial while sines=0
x(m+1,n1+3)=Eq1_1s(n,t,gm,dyna1,F1,prp1,prpp);
% Equation 2 in sigma-theta general layer
x(m+2,n1+3)=Eq2_1s(n,t,gm,dyna1,F1,prp1,prpp);
% Equation 3 in sigma r-theta general layer
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+3,n1+3)=Eq3_1s(n,t,gm,dyna1,F1,prp1,prpp);
% Equation 4 in Dr general layer
x(m+4,n1+3)=Eq4_1s(n,t,gm,dyna1,F1,prp1,prpp);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C Boundary Conditions (known surface traction)
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
t=0;% cosines coefficients in radial while sines=0
x(m+5,n1+3)=sigmarn(n,t,r2,dyna1,F1,prp1);
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+6,n1+3)=sigmartn(n,t,r2,dyna1,F1,prp1);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C! Physical Equations Unknown H
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% Equation 1 in sigma-r general layer
t=0;% cosines coefficients in radial while sines=0
x(m+1,n1+4)=Eq1_1s(n,t,gm,dyna1,H1,prp1,prpp);
% Equation 2 in sigma-theta general layer
x(m+2,n1+4)=Eq2_1s(n,t,gm,dyna1,H1,prp1,prpp);
% Equation 3 in sigma r-theta general layer
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+3,n1+4)=Eq3_1s(n,t,gm,dyna1,H1,prp1,prpp);
% Equation 4 in Dr general layer
x(m+4,n1+4)=Eq4_1s(n,t,gm,dyna1,H1,prp1,prpp);
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%C Boundary Conditions (known surface traction)
%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
t=0;% cosines coefficients in radial while sines=0
x(m+5,n1+4)=sigmarn(n,t,r2,dyna1,H1,prp1);
t=pi/(2*n);% cosines coefficients in radial while sines=0
x(m+6,n1+4)=sigmartn(n,t,r2,dyna1,H1,prp1);
TM{n}=x;
A{n}=TM{n}\b{n};
end
* * *This broadcast message was sent to all bidders on Wednesday Sep 28, 2011 5:10:04 PM:
As some of you have noticed the program calls functions and subrountines. I already have these functions or subroutines in Fortran. The Matlab code is a development on previous code that I had it translated previously. This is an on going work. I code the programs in Matlab then translate them.
* * *This broadcast message was sent to all bidders on Monday Oct 3, 2011 3:45:49 PM:
I used to code my programs then convert them to Fortran. Unfortunately, the same vworker is busy and unavailable. This is an on going project. I would like to pick someone who is a potential future worker. This will ease things for both of us.