fortranでのプログラムについて

質問

初歩的な質問失礼します
func = x**3-3*(x**2)-2*x+14　に対して

real function ymax( minx, maxx)
real, intent(in) :: minx,maxx
real :: mx1,mx2
mx1 = (6.0+sqrt((-6.0)*(-6.0) - 4.0*3.0*(-2)))/(3.0*2.0)
mx2 = (6.0-sqrt((-6.0)*(-6.0) - 4.0*3.0*(-2)))/(3.0*2.0)
! mx2は極大 mx1は極小
if (minx <= mx2 .and. mx2<maxx ) then
　if (func(maxx)<func(mx2) .and. func(minx)<func(mx2)) then
　　ymax=func(mx2)
　elseif (func(mx2)<func(maxx) .and. func(maxx)>func(minx)) then
　　ymax=func(maxx)
　else
　ymaxx = func(minx)
endif
else
　if (func(minx)>func(maxx)) then
　　ymax=func(minx)
　else
　　yamx=func(maxx)
　endif
endif
end function ymax

こういう感じで関数の最小最大を決めてグラフを画面上に変換して表示しようとするのですがどうも最大値がうまく計算されません。
どこに問題があるのでしょうか？

回答No.3

doara_2011 · ベストアンサー

>　if (func(maxx)<func(mx2) .and. func(minx)<func(mx2)) then
>　　ymax=func(mx2)
>　elseif (func(mx2)<func(maxx) .and. func(maxx)>func(minx)) then
>　　ymax=func(maxx)
>　else
>　ymaxx = func(minx)
>endif

でymaxx= ... となっているのが気になる。

お礼

他の場所にあったタイプミスが原因でした…。
何度も見直したつもりでしたがこんな初歩的なミスで申し訳ないです。
本当にありがとうございました。

回答No.2

Tacosan

なんとなく
implicit none
してないのが敗因のような気がする.

お礼

何度も回答頂いて本当にありがとうございました。
もっとしっかり練習してきます(・∀・)

回答No.1

Tacosan

「どうも最大値がうまく計算されません。」とは, 具体的にはどのように「うまく計算されない」のですか?

あと, プログラムにコメントがほしい.... 「こういう感じ」というのをもっときちんと書いてほしいなぁ....

補足

すいませんあんまり長いと見づらいかと思い省いてしまってました…。

! 関数を決定する

real function func(x)
real, intent(in) ::x
func = x**3-3*(x**2)-2*x+14

return
end function func

! yの最大値を設定する

real function ymax( minx, maxx)
real, intent(in) :: minx,maxx
real :: y1,y2,y3
real :: mx1,mx2
!　func(x)の微分関数＝0 の解 =1±±√15/3=1±1.2　くらい
mx1 = (6.0+sqrt((-6.0)*(-6.0) - 4.0*3.0*(-2)))/(3.0*2.0)
mx2 = (6.0-sqrt((-6.0)*(-6.0) - 4.0*3.0*(-2)))/(3.0*2.0)

! mx2が極大　mx1が極小

!極大値が区間に入ってるか？
! 入っている
if (minx <= mx2 .and. mx2<maxx ) then
　if (func(maxx)<func(mx2) .and. func(minx)<func(mx2)) then
　 ymax=func(mx2)
elseif (func(mx2)<func(maxx) .and. func(maxx)>func(minx)) then
ymax=func(maxx)
else
　 ymaxx = func(minx)
! func(minx) func(maxx) func(mx2) の最大値が最大値
endif
! 入っていない
else
if (func(minx)>func(maxx)) then
　 ymax=func(minx)
else
　 yamx=func(maxx)

! func(minx) func(maxx)　の最大値が最大値
endif
endifend function ymax

! yの最小値を設定する

real function ymin( minx, maxx)
real, intent(in) :: minx,maxx
real :: y1,y2,y3
real :: mx1,mx2

mx1 = (6.0+sqrt((-6.0)*(-6.0) - 4.0*3.0*(-2)))/(3.0*2.0)
mx2 = (6.0-sqrt((-6.0)*(-6.0) - 4.0*3.0*(-2)))/(3.0*2.0)

! 極小値は区間に入ってるか？
! 入っている
if (minx <= mx1 .and. mx1<maxx) then
if (func(mx1)<func(minx) .and. func(mx1)<func(maxx)) then
ymin = func(mx1)
elseif (func(mx1)>func(minx) .and. func(minx)<func(maxx)) then
ymin = func(minx)
else
ymin = func(maxx)
! func(minx) func(maxx) func(mx1)の最小値が最小値
endif
! 入っていない
else
if (func(minx)<func(maxx)) then
ymin = func(minx)
else
ymin = func(maxx)
! func(minx) func(maxx)の最小値が最小値
endif
endif
end function ymin

! 画面上のy座標に変換
function transy(y,MINY,MAXY)
real,intent(in) :: MINY, MAXY
real,parameter :: CRTMINY=23.0, CRTMAXY=1
real,intent(in):: y
transy = (CRTMAXY-CRTMINY)/(MAXY-MINY)* (y-MINY) +CRTMINY
return
end function transy

このプログラムに対して 描画部内で
iy = transy(y,miny,maxy)

としたところ、正と負の数、負の数二つ、正の差の小さい二数を入力した時は正しく動くのですが、
正の3以上離れた数２つを入力したときiyが23を超えてしまいエラーが出てしまいます。
またyの値自体は正しく計算されているのに、ymaxを表示したときに　ymax=0.0000000E+00
となってしまいます。
至らない所が多くて恐縮ですが、もしよろしければ回答お願いします。