Aptarimas:Gryno formulė

Iš Wikibooks.
Jump to navigation Jump to search

Ne tas polinomas padalintas (a ir x ne tas pats), bet polinomas gražiai pasidalijo į [keisti]

  • Taikydami Gryno formulę, apskaičiuokime kreivinį integralą

kai L - apskritimas (a>0), apeinamas teigiama kryptimi (prieš laikrodžio rodyklę).

Kadangi skritulyje funkcijos ir bei jų dalinės išvestinės ir yra tolydžios, tai duotajam kreiviniam integralui galima taikyti Gryno formulę.
Turime:
Randame apskritimo y išraišką:
Randame y išvestinę, o paskui ir dy:
Apskritimo spindulys nes pavyzdžiui, kai tai Žinome, kad šis apskritimas liečiasi koordinačių pradžios taške O(0; 0) ir kad ašis Ox dalina apskritimą pusiau. Vadinasi, kai ir tai gauname teisingą lygybę Vadinasi taškas (0; 0) priklauso apskritimui Kitas apskritimo taškas yra (3; 0), kuris yra ant Ox ašies. Įstačius taško (3; 0) koordinates į apskritimo lygtį gauname Žinome, kad taškas (3; 0) yra toliausias taškas ant Ox ašies. Todėl apskritimo spindulys yra
Kadangi apskritimo spindulys yra ir Ox ašis dalina apskritimą per pusę, tai didžiausia y reikšmė gali būti Vadinasi integravimas vyksta pirmame ir ketvirtame ketvirčiuose. Bet, kadangi, skritulio plotas yra vienodas ketvirtame ketviryje kaip ir pirmame, tai užtenka apskaičiuoti skritulio plotą tik pirmame ketvirtyje, o paskui gautą plotą padauginti iš dviejų. Kad apskaičiuoti skritulio plotą pirmame ketvirtyje, turime žinoti integravimo ribas. Nustatome, kad x kinta nuo 0 iki a, o y kinta nuo 0 iki
Taikydami Gryno formulę, integruojame (pasinaudodami internetiniu integratoriumi):
Pasitikriname (įstatydami y ir dy ir pasinaudodami internetiniu integratoriumi):
Pastaba, kad taip gauname dalyba iš nulio ir neįmanoma vietomis išintegruoti įstatant a arba 0. Bet gauname kažką panašesnio į teisingą atsakymą:
Toliau pasinaudojame internetiniu polinomų dalikliu parinkę gauname Toliau integruojame:
Kadangi tai seniau gautas atsakymas įstačius atitinka ir iš to darome išvada, kad Gryno formulė veikia teisingai (išskyrus minusiuką).

Kvaili sutapimai[keisti]

  • Apskaičiuosime darbą jėgos persikeliant materialiam taškui elipse teigiama kryptimi, jeigu jėga kiekviename taške (x; y) elipsės nukreipta į elipsės centrą ir pagal dydį lygi atstumui nuo taško (x; y) iki elipsės centro.
Pagal sąlyga, Jėgos F(x, y) koordinatės tokios: [ženklas "" paaiškinamas tuo, kad jėga nukreipta į tašką (0; 0)]. Pagal formulę turime kur L - elipsė Todėl

Jei t keistusi nuo 0 iki integralas butu lygus Tarkime, jei tai padarytas darbas pirmame ketvirtyje yra
Kad patikrinti ar apskaičiuota teisingai reikia sudėti visas x reikšmes ant elipsės linijos pirmame ketvirtyje. Taip pat reikia sudėti visas y reikšmes ant elipsės linijos pirmame ketvirtyje. Galiausiai reikia sudėti sumas x ir y reikšmių, kad gauti darbą A pirmame ketvirtyje.
Arba tiesiog darbas yra lygus čia yra visos x reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki a; analogiškai yra visos y reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki b.
Kuris iš šių variantų yra darbas A, paaiškės pasumavus ir paskaičiavus. Sutapimas ar ne, bet pagal pirmą variantą. Na, o pagal antrą variantą turime:
Na, pagal pirmą variantą tikrai ne sutapimas, nes, jei tai ir
Bet ne, vis dėl to tai sutapimas, nes jei tai ir
Integruojant gauname:

Pastebime, kad

Per daug kodų[keisti]

  • Apskaičiuosime darbą jėgos persikeliant materialiam taškui elipse teigiama kryptimi, jeigu jėga kiekviename taške (x; y) elipsės nukreipta į elipsės centrą ir pagal dydį lygi atstumui nuo taško (x; y) iki elipsės centro.
Pagal sąlyga, Jėgos F(x, y) koordinatės tokios: [ženklas "" paaiškinamas tuo, kad jėga nukreipta į tašką (0; 0)]. Pagal formulę turime kur L - elipsė Todėl

Jei t keistusi nuo 0 iki integralas butu lygus
Tarkime, jei tai padarytas darbas pirmame ketvirtyje yra
Kad patikrinti ar apskaičiuota teisingai reikia sudėti visas x reikšmes ant elipsės linijos pirmame ketvirtyje. Taip pat reikia sudėti visas y reikšmes ant elipsės linijos pirmame ketvirtyje. Galiausiai reikia sudėti sumas x ir y reikšmių, kad gauti darbą A pirmame ketvirtyje.
Arba tiesiog darbas yra lygus čia yra visos x reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki a; analogiškai yra visos y reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki b.
Kuris iš šių variantų yra darbas A, paaiškės pasumavus ir paskaičiavus. Sutapimas ar ne, bet pagal pirmą variantą. Na, o pagal antrą variantą turime:
Iš tikro, ko gero, mes apskaičiuojame tokiu budu ne darbą atlikta apeinant elipsės liniją pirmame ketvirtyje, o darbą atlikta apeinant tiesę pirmame ketvirtyje. Štai kodas programos "Free Pascal" (FreePascal IDE for Win32 for i386; Target CPU: i386; Version 1.0.12 2011/04/23; <Compiler Version 2.4.4>; <Debugger GDB 7.2>; Copyright <C> 1998-2009):
     var
     a:longint;
     b,d,c:real;
     begin
     c:=0;
       for a:=1 to 10
       do
       c:=c+0.1*sqrt(sqr(a*0.3)+sqr((10-a)*0.5)); 
       writeln(c);
     readln;
     end.
Gauname atsakymą
Arba darome taip (pradedame ne nuo 1, o nuo 0, bet tada reikėtų dalinti iš 11, o ne iš 10 kol a iki mažo skaičiaus skaičiuojamas):
     var
     a:longint;
     b,d,c:real;
     begin
     c:=0;
       for a:=0 to 10
       do
       c:=c+0.1*sqrt(sqr(a*0.3)+sqr((10-a)*0.5)); 
       writeln(c);
     readln;
     end.
Ir gauname atsakymą
Panaudojus šį kodą:
    var
    a:longint;
    c:real;
    begin
      for a:=1 to 10
      do
      c:=c+0.1*sqrt(sqr(a*0.3)+sqr((11-a)*0.5));
      writeln(c);
    readln;
    end.
gauname
Kad gauti tikslesnį atsakymą naudojame šį kodą:
     var
     a:longint;
     b,d,c:real;
     begin
     c:=0;
       for a:=1 to 100
       do
       c:=c+0.01*sqrt(sqr(a*0.03)+sqr((100-a)*0.05)); 
       writeln(c);
     readln;
     end.
Tokiu budu gauname gana tikslų atsakymą
Kad gauti labai tikslų atsakymą naudojame tokį kodą:
   var
   a:longint;
   c:real;
   begin
       for a:=1 to 1000000
       do
       c:=c+0.000001*sqrt(sqr(a*0.000003)+sqr((1000000-a)*0.000005)); 
       writeln(c);
   readln;
   end.
Gauname atsakymą
Kad gauti labai labai tikslų atsakymą naudojame šį kodą:
   var
   a:longint;
   c:real;
   begin
       for a:=1 to 1000000000  do
       c:=c+0.000000001*sqrt(sqr(a*0.000000003)+sqr((1000000000-a)*0.000000005)); 
       writeln(c);
   readln;
   end.
Gauname atsakymą po 22 sekundžių ant 2.6 GHz procesoriaus.
Jei naudojame šį kodą:
   var
   a:longint;
   c:real;
   begin
       for a:=0 to 1000000000  do
       c:=c+0.000000001*sqrt(sqr(a*0.000000003)+sqr((1000000000-a)*0.000000005)); 
       writeln(c);
   readln;
   end.
tai gauname atsakymą
Panaudojus šį kodą:
   var
   a:longint;
   c:real;
   begin
   for a:=1 to 1000000000  do
   c:=c+0.000000001*sqrt(sqr(a*0.000000003)+sqr((1000000001-a)*0.000000005));
   writeln(c);
   readln;
   end.
gauname atsakymą
Integruojant gauname:
Pastebime, kad Apskaičiuojame viską nepriekaištingai tiksliai su kompiuterio kalkuliatoriumi:
Randame, kad
Toliau bandome rasti atlikta darbą apeita tiese integruojant Randame Turime, kad x integravimo ribos yra 0 iki 5, o y integravimo ribos yra nuo 0 iki 3. Gauname:
Akivaizdu, kad atsakymas yra neteisingas, o atsakymas A=8 elipsei yra teisingas, nes elipsės lankas yra ilgesnis nei tiesės atkarpa pirmame ketvirtyje ir pagal sąlygą jėgos dydis proporcingas atstumui nuo centro, todėl elipsės lanko taškai yra toliau nuo centro nei tiesės atkarpos taškai. Ir todėl elipsės darbas A=8 yra didesnis už darbą A=3.2865795 tiesės atkarpa.
Pastebime, kad jeigu elipsė būtų apskritimas tai pirmame ketvirtyje atliktas darbas A=0. Taip yra todėl, kad kai x didėja, tada y reikšmės mažėja. Todėl skaičiuoti darbą, ko gero, galėtų būti teisingiau arba Tuomet akivaizdu, kad
Tą patį darba gauname ir apeinant tiese:
Lygiai tą patį darbą gausime ir integruojant taip:
Tą patį gausime ir taip integruojant nuo 0 iki 5:

O visgi vadovėliuose nėra kreivės masės skaičiavimo formulių (todėl ir atsakymas neteisingas)[keisti]

  • Nustatyti tiesės masę tik pirmame ketvirtyje. Tiesės tankis tolstant tiesės taškams nuo centro (koordinačių pradžios taško O) didėja proporcingai, t. y.
Sprendimas. Pasinaudosime masės skaičiavimo formule
Toliau integruodami gauname:

X failai (skyrius Jėgos darbas)[keisti]

Jėgos darbas padarytas judant kreive plokštumoje apskaičiuojamas pagal formulę Jėgos darbas padarytas judant erdvine kreive apskaičiuojamas taip:

  • Apskaičiuosime darbą jėgos persikeliant materialiam taškui elipse teigiama kryptimi, jeigu jėga kiekviename taške (x; y) elipsės nukreipta į elipsės centrą ir pagal dydį lygi atstumui nuo taško (x; y) iki elipsės centro.
Pagal sąlyga, Jėgos F(x, y) koordinatės tokios: [ženklas "" paaiškinamas tuo, kad jėga nukreipta į tašką (0; 0)]. Pagal formulę turime kur L - elipsė Todėl

Jei t keistusi nuo 0 iki integralas butu lygus
Tarkime, jei tai padarytas darbas pirmame ketvirtyje yra
Kad patikrinti ar apskaičiuota teisingai reikia sudėti visas x reikšmes ant elipsės linijos pirmame ketvirtyje. Taip pat reikia sudėti visas y reikšmes ant elipsės linijos pirmame ketvirtyje. Galiausiai reikia sudėti sumas x ir y reikšmių, kad gauti darbą A pirmame ketvirtyje.
Arba tiesiog darbas yra lygus čia yra visos x reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki a; analogiškai yra visos y reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki b.
Kuris iš šių variantų yra darbas A, paaiškės pasumavus ir paskaičiavus. Sutapimas ar ne, bet pagal pirmą variantą. Na, o pagal antrą variantą turime:
Iš tikro, ko gero, mes apskaičiuojame tokiu budu ne darbą atlikta apeinant elipsės liniją pirmame ketvirtyje, o darbą atlikta apeinant tiesę pirmame ketvirtyje. Štai kodas programos "Free Pascal" (FreePascal IDE for Win32 for i386; Target CPU: i386; Version 1.0.12 2011/04/23; <Compiler Version 2.4.4>; <Debugger GDB 7.2>; Copyright <C> 1998-2009):
    var
    a:longint;
    c:real;
    begin
      for a:=1 to 10
      do
      c:=c+0.1*sqrt(sqr(a*0.3)+sqr((11-a)*0.5));
      writeln(c);
    readln;
    end.
gauname
Panaudojus šį kodą:
   var
   a:longint;
   c:real;
   begin
   for a:=1 to 1000000000  do
   c:=c+0.000000001*sqrt(sqr(a*0.000000003)+sqr((1000000001-a)*0.000000005));
   writeln(c);
   readln;
   end.
gauname atsakymą po 22 sekundžių ant 2.6 GHz procesoriaus.
Panaudojus šį (teisingesnį) kodą:
   var
   a:longint;
   c:real;
   begin
   for a:=0 to 1000000000  do
   c:=c+sqrt(sqr(a*0.000000003)+sqr((1000000000-a)*0.000000005))/1000000001;
   writeln(c);
   readln;
   end.
gauname atsakymą po 25 sekundžių ant 2.6 GHz procesoriaus.
Integruojant gauname:
(Update: Perskaičiavus gaunasi:
(9,1911764705882352941176470588235+3,6425041506487550296517491344981)-(-3,3088235294117647058823529411765-1,6139227719443049309711137193553)=
=12,833680621236990323769396193322-(-4,9227463013560696368534666605317)=17,756426922593059960622862853854.
Update 2: Ne -3,3088, o -1,9852941176470588235294117647059.)
Pastebime, kad Apskaičiuojame viską nepriekaištingai tiksliai su kompiuterio kalkuliatoriumi:
Randame, kad
Toliau bandome rasti atlikta darbą apeita tiese integruojant Randame Turime, kad x integravimo ribos yra 0 iki 5, o y integravimo ribos yra nuo 0 iki 3. Gauname:
Akivaizdu, kad atsakymas yra neteisingas, o atsakymas A=8 elipsei yra teisingas, nes elipsės lankas yra ilgesnis nei tiesės atkarpa pirmame ketvirtyje ir pagal sąlygą jėgos dydis proporcingas atstumui nuo centro, todėl elipsės lanko taškai yra toliau nuo centro nei tiesės atkarpos taškai. Ir todėl elipsės darbas A=8 yra didesnis už darbą A=3.2865795 tiesės atkarpa.
Pastebime, kad jeigu elipsė būtų apskritimas tai pirmame ketvirtyje atliktas darbas A=0. Taip yra todėl, kad kai x didėja, tada y reikšmės mažėja. Todėl skaičiuoti darbą, ko gero, galėtų būti teisingiau arba Tuomet akivaizdu, kad
Tą patį darba gauname ir apeinant tiese:
Lygiai tą patį darbą gausime ir integruojant taip:
Tą patį gausime ir taip integruojant nuo 0 iki 5:


  • Nustatyti tiesės masę tik pirmame ketvirtyje. Tiesės tankis tolstant tiesės taškams nuo centro (koordinačių pradžios taško O) didėja proporcingai, t. y.
Sprendimas. Pasinaudosime masės skaičiavimo formule
Toliau integruodami gauname:

Perskaičiavus gaunas m=41.4470084394957645512263519847-(-4.015301738572567547333986879663)=45.462310178068332098560338864363.
Wolframo integratoriui varžteliai atsisuko, nes integruojant tą patį taip ir taip nuo 0 iki 5 gauname skirtingus atsakymus atitinkamai:
ir
Beje, integruojant taip nuo 0 iki 5 gauname:
(13,865120473473883114631633205975+5,0555276680167165753752595311952)-(-0,66176470588235294117647058823529-0,64871162108625555582601097623286)=
=18,92064814149059969000689273717-(-1,3104763269686084970024815644681)=20,231124468459208187009374301638.

Pilna klaidų, bet vis tiek kažkas ne taip integruojant su ir be[keisti]

  • Nustatyti tiesės masę tik pirmame ketvirtyje. Tiesės tankis tolstant tiesės taškams nuo centro (koordinačių pradžios taško O) didėja proporcingai, t. y.
Sprendimas. Pasinaudosime masės skaičiavimo formule
Toliau integruodami gauname:

Wolframo integratoriui varžteliai atsisuko, nes integruojant tą patį taip ir taip nuo 0 iki 5 gauname skirtingus atsakymus atitinkamai:
ir
Beje, integruojant taip nuo 0 iki 5 gauname:
(13,865120473473883114631633205975+5,0555276680167165753752595311952)-(-0,66176470588235294117647058823529-0,64871162108625555582601097623286)=
=18,92064814149059969000689273717-(-1,3104763269686084970024815644681)=20,231124468459208187009374301638. Na, bent jau apytiksliai gauname kažką panašaus į :
Nors šita formulė ir skirta kreivės svorio centro koordinačių nustatymui (), bet sprendžiant iš kai kurių matematikos vadovelių, kreivės masę galima apskaičiuoti, bet tada, matyt, reikia pačiam sugalvoti formulę, pagal kurią skaičiuoti (nebent ten turima galvoje, kad kreivės ilgis lygus masei, o tankis toks pat visuose kreivės taškuose).
Tačiau, ko gero,

Wolframo integratorius integruoja skirtingai išskleistą ir neišskleistą funkciją ir skirtingai su skirtingomis konstantomis[keisti]

(Taip 1.16619*sqrt(9-18*x/5 +(34*x^2)/25) integruojant gauname 19,163865128930380119022361244004; čia 1,16619=(34/25)^(1/2).)
  • Nustatyti tiesės masę tik pirmame ketvirtyje. Tiesės tankis tolstant tiesės taškams nuo centro (koordinačių pradžios taško O) didėja proporcingai, t. y.
Sprendimas. Pasinaudosime masės skaičiavimo formule
Toliau integruodami gauname:

Wolframo integratoriui varžteliai atsisuko, nes integruojant tą patį taip ir taip nuo 0 iki 5 gauname skirtingus atsakymus atitinkamai:
ir
Beje, integruojant taip nuo 0 iki 5 gauname:
[34^0.5]*[(45,955882352941176470588235294118+18,21252075324377514825874567249)-(-9,9264705882352941176470588235294-8,0696138597215246548555685967763)]/25=
=[34^0.5]*[64,168403106184951618846980966608-(-17,996084447956818772502627420306)]/25=
=[34^0.5]*82,164487554141770391349608386914/25=19,163886975712642545064105297031. Na, bent jau apytiksliai gauname kažką panašaus į :
Nors šita formulė ir skirta kreivės svorio centro koordinačių nustatymui (), bet sprendžiant iš kai kurių matematikos vadovelių, kreivės masę galima apskaičiuoti, bet tada, matyt, reikia pačiam sugalvoti formulę, pagal kurią skaičiuoti (nebent ten turima galvoje, kad kreivės ilgis lygus masei, o tankis toks pat visuose kreivės taškuose).
Tačiau, ko gero,

Wolframo integratorius tikrai turi problemų tame kaip pateikti integravimo funkciją[keisti]

  • Nustatyti tiesės masę tik pirmame ketvirtyje. Tiesės tankis tolstant tiesės taškams nuo centro (koordinačių pradžios taško O) didėja proporcingai, t. y.
Sprendimas. Pasinaudosime masės skaičiavimo formule
Toliau integruodami gauname:

Wolframo integratoriui varžteliai atsisuko, nes integruojant tą patį taip ir taip nuo 0 iki 5 gauname skirtingus atsakymus atitinkamai:
ir
Beje, integruojant taip nuo 0 iki 5 gauname:
[34^0.5]*[(45,955882352941176470588235294118+18,21252075324377514825874567249)-(-9,9264705882352941176470588235294-8,0696138597215246548555685967763)]/25=
=[34^0.5]*[64,168403106184951618846980966608-(-17,996084447956818772502627420306)]/25=
=[34^0.5]*82,164487554141770391349608386914/25=19,163886975712642545064105297031.
Wolframo integratorius duoda keistus sutapimus:
Nors formulė () ir skirta kreivės svorio centro koordinatės nustatymui (), bet sprendžiant iš kai kurių matematikos vadovelių, kreivės masę galima apskaičiuoti, bet tada, matyt, reikia pačiam sugalvoti formulę, pagal kurią skaičiuoti (nebent ten turima galvoje, kad kreivės ilgis lygus masei, o tankis toks pat visuose kreivės taškuose).
Tačiau, ko gero,

Uždavinys baigtas spręsti[keisti]

  • Nustatyti tiesės masę tik pirmame ketvirtyje. Tiesės tankis tolstant tiesės taškams nuo centro (koordinačių pradžios taško O) didėja proporcingai, t. y.
Sprendimas. Pasinaudosime masės skaičiavimo formule
Toliau integruodami gauname:

Wolframo integratoriui varžteliai atsisuko, nes integruojant tą patį taip ir taip nuo 0 iki 5 gauname skirtingus atsakymus atitinkamai:
ir
Beje, integruojant taip nuo 0 iki 5 gauname:
[34^0.5]*[(45,955882352941176470588235294118+18,21252075324377514825874567249)-(-9,9264705882352941176470588235294-8,0696138597215246548555685967763)]/25=
=[34^0.5]*[64,168403106184951618846980966608-(-17,996084447956818772502627420306)]/25=
=[34^0.5]*82,164487554141770391349608386914/25=19,163886975712642545064105297031.
Wolframo integratorius duoda keistus sutapimus:
Nors formulė () ir skirta kreivės svorio centro koordinatės nustatymui (), bet sprendžiant iš kai kurių matematikos vadovelių, kreivės masę galima apskaičiuoti, bet tada, matyt, reikia pačiam sugalvoti formulę, pagal kurią skaičiuoti (nebent ten turima galvoje, kad kreivės ilgis lygus masei, o tankis toks pat visuose kreivės taškuose).
Tačiau, ko gero,
Kad išsiaiškinti teisybę pasinaudosime integralų lentele:
Tada integruojame:
Vadinasi, randame tiesės masę kaip prašyta sąlygoje:

Iš skyriaus Jėgos darbas neteisingi atsakymai (iš tikro masė yra p, o ne m)[keisti]

Jėgos darbas padarytas judant kreive plokštumoje apskaičiuojamas pagal formulę Jėgos darbas padarytas judant erdvine kreive apskaičiuojamas taip:

  • Apskaičiuosime darbą jėgos persikeliant materialiam taškui elipse teigiama kryptimi, jeigu jėga kiekviename taške (x; y) elipsės nukreipta į elipsės centrą ir pagal dydį lygi atstumui nuo taško (x; y) iki elipsės centro.
Pagal sąlyga, Jėgos F(x, y) koordinatės tokios: [ženklas "" paaiškinamas tuo, kad jėga nukreipta į tašką (0; 0)]. Pagal formulę turime kur L - elipsė Todėl

Jei t keistusi nuo 0 iki integralas butu lygus
Tarkime, jei tai padarytas darbas pirmame ketvirtyje yra
Kad patikrinti ar apskaičiuota teisingai reikia sudėti visas x reikšmes ant elipsės linijos pirmame ketvirtyje. Taip pat reikia sudėti visas y reikšmes ant elipsės linijos pirmame ketvirtyje. Galiausiai reikia sudėti sumas x ir y reikšmių, kad gauti darbą A pirmame ketvirtyje.
Arba tiesiog darbas yra lygus čia yra visos x reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki a; analogiškai yra visos y reikšmės ant elipsės linijos pirmame ketvirtyje nuo 0 iki b.
Kuris iš šių variantų yra darbas A, paaiškės pasumavus ir paskaičiavus. Sutapimas ar ne, bet pagal pirmą variantą. Na, o pagal antrą variantą turime:
Iš tikro, ko gero, mes apskaičiuojame tokiu budu ne darbą atlikta apeinant elipsės liniją pirmame ketvirtyje, o darbą atlikta apeinant tiesę pirmame ketvirtyje. Štai kodas programos "Free Pascal" (FreePascal IDE for Win32 for i386; Target CPU: i386; Version 1.0.12 2011/04/23; <Compiler Version 2.4.4>; <Debugger GDB 7.2>; Copyright <C> 1998-2009):
    var
    a:longint;
    c:real;
    begin
      for a:=1 to 10
      do
      c:=c+0.1*sqrt(sqr(a*0.3)+sqr((11-a)*0.5));
      writeln(c);
    readln;
    end.
gauname
Panaudojus šį kodą:
   var
   a:longint;
   c:real;
   begin
   for a:=1 to 1000000000  do
   c:=c+0.000000001*sqrt(sqr(a*0.000000003)+sqr((1000000001-a)*0.000000005));
   writeln(c);
   readln;
   end.
gauname atsakymą po 22 sekundžių ant 2.6 GHz procesoriaus.
Panaudojus šį (teisingesnį) kodą:
   var
   a:longint;
   c:real;
   begin
   for a:=0 to 1000000000  do
   c:=c+sqrt(sqr(a*0.000000003)+sqr((1000000000-a)*0.000000005))/1000000001;
   writeln(c);
   readln;
   end.
gauname atsakymą po 25 sekundžių ant 2.6 GHz procesoriaus.
Integruojant gauname:
Pastebime, kad Apskaičiuojame viską nepriekaištingai tiksliai su kompiuterio kalkuliatoriumi:
Randame, kad
Toliau bandome rasti atlikta darbą apeita tiese integruojant Randame Turime, kad x integravimo ribos yra 0 iki 5, o y integravimo ribos yra nuo 0 iki 3. Gauname:
Akivaizdu, kad atsakymas yra neteisingas, o atsakymas A=8 elipsei yra teisingas, nes elipsės lankas yra ilgesnis nei tiesės atkarpa pirmame ketvirtyje ir pagal sąlygą jėgos dydis proporcingas atstumui nuo centro, todėl elipsės lanko taškai yra toliau nuo centro nei tiesės atkarpos taškai. Ir todėl elipsės darbas A=8 yra didesnis už darbą A=3.2865795 tiesės atkarpa.
Pastebime, kad jeigu elipsė būtų apskritimas tai pirmame ketvirtyje atliktas darbas A=0. Taip yra todėl, kad kai x didėja, tada y reikšmės mažėja. Todėl skaičiuoti darbą, ko gero, galėtų būti teisingiau arba Tuomet akivaizdu, kad
Tą patį darba gauname ir apeinant tiese:
Lygiai tą patį darbą gausime ir integruojant taip:
Tą patį gausime ir taip integruojant nuo 0 iki 5:


  • Nustatyti tiesės masę tik pirmame ketvirtyje. Tiesės tankis tolstant tiesės taškams nuo centro (koordinačių pradžios taško O) didėja proporcingai, t. y.
Sprendimas. Pasinaudosime masės skaičiavimo formule