How to get numerical output

7 Ansichten (letzte 30 Tage)
okoth ochola am 4 Jun. 2023
Beantwortet: KALYAN ACHARJYA am 4 Jun. 2023
Hello, I have been facing some problems with integration, matlab is not giving me numerical output but an output of the form given below. I have also attached part of the code responsible for the output (the loop is just a part of a larger script am using, so the Monthly_parameters(2,n) is caled from elsewhere in he code, which has no problem). How can I resolve this problem to get numerical output. Please can someone assist. What am I not doing well,
for n=1:c
syms x
expr=(Monthly_parameters(1,n)/Monthly_parameters(2,n))*((x/Monthly_parameters(2,n))^(Monthly_parameters(1,n)-1))*exp(-(x/Monthly_parameters(2,n))^Monthly_parameters(1,n))*density*x^3
Power_density_monthly(n,:) = int(expr,[0 25]);
if n==c
display(Power_density_monthly)
Max_power_monthly=0.5*density.*(Wind_speed_MaxE.^3)
Max_Power_annual=0.5*density.*(Wind_speed_MaxE_annual.^3)
Mp_power_monthly=0.5*density.*(Most_probable_wind_speed.^3)
syms x
expr=(Overal_parameters(1)/Overal_parameters(2))*((x/Overal_parameters(2))^(Overal_parameters(1)-1))*exp(-(x/Overal_parameters(2))^Overal_parameters(1))*density*x^3
Mp_power_Overall=int(expr,[0 25])
end
end
Power_density_monthly =
int((780035855555123397*x^3*exp(-((4503599627370496*x)/7214917075124989)^(4843579468129569/1125899906842624))*((4503599627370496*x)/7214917075124989)^(3717679561286945/1125899906842624))/225179981368524800, x, 0, 25)
int((876700907771204457*x^3*exp(-((2251799813685248*x)/3500040294331183)^(330107227435517/70368744177664))*((2251799813685248*x)/3500040294331183)^(259738483257853/70368744177664))/225179981368524800, x, 0, 25)
int((1011422054596262133*x^3*exp(-((4503599627370496*x)/6407789261932345)^(5577777155503949/1125899906842624))*((4503599627370496*x)/6407789261932345)^(4451877248661325/1125899906842624))/225179981368524800, x, 0, 25)
int((963101305032312717*x^3*exp(-((140737488355328*x)/203114773723167)^(1346864763691579/281474976710656))*((140737488355328*x)/203114773723167)^(1065389786980923/281474976710656))/225179981368524800, x, 0, 25)
int((899082092550608151*x^3*exp(-((281474976710656*x)/432556907928793)^(669411416785461/140737488355328))*((281474976710656*x)/432556907928793)^(528673928430133/140737488355328))/225179981368524800, x, 0, 25)
int((787803615475959759*x^3*exp(-((4503599627370496*x)/7016859725054415)^(4757527240969181/1125899906842624))*((4503599627370496*x)/7016859725054415)^(3631627334126557/1125899906842624))/225179981368524800, x, 0, 25)
int((388299016464805677*x^3*exp(-((4503599627370496*x)/7550763379329251)^(1261675515893479/281474976710656))*((4503599627370496*x)/7550763379329251)^(980200539182823/281474976710656))/112589990684262400, x, 0, 25)
int((109902315657116823*x^3*exp(-((4503599627370496*x)/7141945877977533)^(2702115500474055/562949953421312))*((4503599627370496*x)/7141945877977533)^(2139165547052743/562949953421312))/28147497671065600, x, 0, 25)
int((932214374257080879*x^3*exp(-((4503599627370496*x)/7077114042901685)^(5677962354422989/1125899906842624))*((4503599627370496*x)/7077114042901685)^(4552062447580365/1125899906842624))/225179981368524800, x, 0, 25)
int((30733691381488647*x^3*exp(-((4503599627370496*x)/7193181261676949)^(3044221384155719/562949953421312))*((4503599627370496*x)/7193181261676949)^(2481271430734407/562949953421312))/7036874417766400, x, 0, 25)
int((971959284150080301*x^3*exp(-((2251799813685248*x)/3251016897941623)^(2719492208413129/562949953421312))*((2251799813685248*x)/3251016897941623)^(2156542254991817/562949953421312))/225179981368524800, x, 0, 25)
int((900166504076381007*x^3*exp(-((4503599627370496*x)/6159683216945655)^(1193003600218407/281474976710656))*((4503599627370496*x)/6159683216945655)^(911528623507751/281474976710656))/225179981368524800, x, 0, 25)
>>
0 Kommentare-2 ältere Kommentare anzeigen-2 ältere Kommentare ausblenden

Melden Sie sich an, um zu kommentieren.

Akzeptierte Antwort

KALYAN ACHARJYA am 4 Jun. 2023
You can use vpa
https://in.mathworks.com/help/symbolic/vpa.html
0 Kommentare-2 ältere Kommentare anzeigen-2 ältere Kommentare ausblenden

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Number Theory finden Sie in Help Center und File Exchange

R2018a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by