دانلود پروژه کدنویسی: روش رانگ کوتای مرتبه چهار و بهره گیری از شوتینگ
قیمت:900000ریال
عنوان:
پروژه محاسبات عددی پیشرفته
روش رانگ کوتای مرتبه چهار و بهره گیری از شوتینگ
فرمت فایل اول ورد 9 صفحه
فرمت فایل دوم پی دی اف 9 صفحه
این پروژه که کدنویسی یک مقاله توسط برنامه فرترن میباشد(عکس مقاله نیز گذاشته شده است).در این گزارش با استفاده از روش رانج کوتای مرتبه چهار و شوتینگ غیر خطی این دستگاه معادلات حل شدهاند.
این یک کار منحصربفرد میباشد که در هیچ جای اینترنت قابل دسترس نیست،چراکه بنده آن را خود انجام داده و کپی برداری نیست.
به شما اطمینان خاطر میدهم که در ضمن منحصربفرد بودن ،یک پروژه ی کامل و بی عیب و نقض را دریافت میکنید.
قابل ذکر است که فایل مقاله و پروژه ی انجام شده کامل را میتوانید ،در یک فایل زیپ دریافت کنید.
پاره ای از حل این پروژه کامل را در زیر به نمایش گذاشته ام:
ابتدا معادله شماره 10 مقاله را به سه معادله و معادله شماره 11 را به دو معادله تبدیل میکنیم.با رانگ کوتای مرتبه چهار که کد آن را در بالا آوردیم،F , F’ , F’” را بدست می آوریم و F ها را در معادله 11 جاگذاری میکنیم.F’’’ را با شوتینگ بدست می آوریم و برای معادله 11 هم از شوتینگ بهره میگیریم.(چون’ θ را نداریم از شوتینگ میرویم)
با بهره گیری از نرم افزار گیت دیتا و با توجه به نمودارهای داده شده در مقاله،چند نقطه را از آن نمودارها در نظر گرفتم و نمودارهای مورد نیاز را توسط اکسل رسم کردم و با نمودار مقاله کنار در یک محور مختصات رسم کرده ام.
کد نویسی رانگ کوتای مرتبه چهار و شوتینگ در برنامه فرترن بشرح زیر انجام شده است :
Program rungkuta
implicit none
real::h,c,g,gg,n1,n2,a,xk,beta,q1,q2,u,gama,qq1,qq2
integer::b,d,m,i,j,x,bbb,bb
real,dimension(1000,1000)::k,k1
real,dimension(1000,1000)::w,w1,ww,ww1
real,dimension(0:1000)::yk
real,dimension(1000)::pp
open(23,FILE='23.txt')
!!!!!!!!!!!!!!!!!!!!!!**********************************************************
!!!!!!!!!!!!!!!!!!!!!!**********************************************************
!!!!!!!!!!!!!!!!!!!!!!**********************************************************
!!!!!!!!!!!!!!!!!!!!!!**********************************************************
a=0
c=3
beta=2
m=3
h=0.1
n1=(c-a)/h
!!!!!!
xk=1.75
11 w(1,1)=0
w(2,1)=1
w(3,1)=xk
k(:,:)=0
d=0
do j=1,n1-1
d=d+1
b=0
do i=1,m-1
b=b+1
k(1,b)=h*w(b+1,d)
end do
k(1,m)=h*g(w(1,d),w(2,d),w(3,d))
b=0
do i=1,m-1
b=b+1
k(2,b)=h*(w(b+1,d)+k(1,b+1)/2)
end do
k(2,m)=h*g(w(1,d)+k(1,1)/2,w(2,d)+k(1,2)/2,w(3,d)+k(1,3)/2)
b=0
do i=1,m-1
b=b+1
k(3,b)=h*(w(b+1,d)+k(2,b+1)/2)
end do
k(3,m)=h*g(w(1,d)+k(2,1)/2,w(2,d)+k(2,2)/2,w(3,d)+k(2,3)/2)
b=0
do i=1,m-1
b=b+1
k(4,b)=h*(w(b+1,d)+k(3,b+1))
end do
k(4,m)=h*g(w(1,d)+k(3,1),w(2,d)+k(3,2),w(3,d)+k(3,3))
b=0
do i=1,m
b=b+1
w(b,d+1)=w(b,d)+(k(1,b)+2*k(2,b)+2*k(3,b)+k(4,b))/6
end do
end do
w1(1,1)=0
w1(2,1)=0
w1(3,1)=1
k(:,:)=0
d=0
do j=1,n1-1
d=d+1
b=0
do i=1,m-1
.....
