# Convert a working Matlab to Fortran

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!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

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}(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}(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.

Навички: Техніка, Архітектура ПЗ, Тестування ПЗ

Про роботодавця:
( 105 відгуки(-ів) ) Orlando, Canada

ID Проекту: #3603634

## 4 фрілансерів(-и) у середньому готові виконати цю роботу за \$69

koles23

See private message.

\$50.15 USD за 10 дні(-в)
(17 відгуків(и))
6.0
matfizvw

See private message.

\$90 USD за 10 дні(-в)
(46 відгуків(и))
5.9
spiderdevelopers

See private message.

\$100.3 USD за 10 дні(-в)
(16 відгуків(и))
4.3
hervebags

See private message.

\$34 USD за 10 дні(-в)
(4 відгуків(и))
3.4