Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- /*
- * To change this license header, choose License Headers in Project Properties.
- * To change this template file, choose Tools | Templates
- * and open the template in the editor.
- */
- package kwadratura.gaussa.legendre.a;
- import java.util.Scanner;
- /**
- *
- * @author Student
- */
- public class KwadraturaGaussaLegendreA {
- public static double caleczka (double x)
- {
- return Math.sin(1.1*x -0.3)/(2.5+Math.pow(x, 2));
- }
- /**
- * @param args the command line arguments
- */
- public static void main(String[] args) {
- // TODO code application logic here
- double a=0.4;
- double b=2.2;
- double wynik,s=0;
- int n;
- Scanner el_wejsciowy =new Scanner(System.in);
- System.out.println("Podaj n");
- n=el_wejsciowy.nextInt();
- switch (n)
- {
- case 2:
- {
- double tab2[][]={{1,-0.57735027},
- {1,0.57735027}};
- for (int i=0;i<tab2.length;i++)
- {
- s += (tab2[i][0]* caleczka((((b-a)/2)*tab2[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 2 wynosi: "+ wynik);
- break;
- }
- case 3:
- {
- double tab3[][]={
- {0.88888889,0.0},
- {0.55555555,-0.77459667},
- {0.55555555,0.77459667}};
- for (int i=0;i<tab3.length;i++)
- {
- s += (tab3[i][0]* caleczka((((b-a)/2)*tab3[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 3 wynosi: "+ wynik);
- break;
- }
- case 4:
- {
- double tab4[][]={{0.65214515,0.33998104},
- {0.65214515,-0.33998104},
- {0.34785485,0.86113631},
- {0.34785485,-0.86113631}};
- for (int i=0;i<tab4.length;i++)
- {
- s += (tab4[i][0]* caleczka((((b-a)/2)*tab4[i][1])+((b+a)/2)));
- }
- wynik=((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 4 wynosi: "+ wynik);
- break;
- }
- case 5:
- {
- double tab5[][]={{0.56888889,0.0},
- {0.47862867,-0.53846931},
- {0.47862867,0.53846931},
- {0.23692689,-0.90617985},
- {0.23692689,0.90617985}};
- /*double tab5[][]={{0.0,0.0,0.53846931,-0.53846931,0.90617985,-0.90617985},
- {0.56888889,-0.56888889,0.47862867,-0.47862867,0.23692689,-0.23692689 }};*/
- for (int i=0;i<tab5.length;i++)
- {
- s += (tab5[i][0]* caleczka((((b-a)/2)*tab5[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 5 wynosi: "+ wynik);
- break;
- }
- case 6:
- {
- double tab6[][]={{0.46791393,-0.23861918},
- {0.46791393,0.23861918},
- {0.36076157,-0.66120939},
- {0.36076157,0.66120939},
- {0.17132449,-0.93246951},
- {0.17132449,0.93246951 }};
- /* double tab6[][]={{0.23861918,-0.23861918,0.66120939,-0.66120939,0.93246951,-0.93246951},
- {0.46791393,-0.46791393,0.36076157,-0.36076157,0.17132449,-0.17132449 }};*/
- for (int i=0;i<tab6.length;i++)
- {
- s += (tab6[i][0]* caleczka((((b-a)/2)*tab6[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 6 wynosi: "+ wynik);
- break;
- }
- case 7:
- {
- double tab7[][]={{0.41795918,0.0},
- {0.38183005,-0.40584515},
- {0.38183005,0.40584515},
- {0.27970539,-0.74153119},
- {0.27970539,0.74153119},
- {0.12948497,-0.94910791},
- {0.12948497,0.94910791}};
- /* double tab7[][]={{0.0,0.0,0.40584515,-0.40584515,0.74153119,-0.74153119,0.94910791,-0.94910791},
- {0.41795918,-0.41795918,0.38183005,-0.38183005,0.27970539,-0.27970539,0.12948497,-0.12948497}};*/
- for (int i=0;i<tab7.length;i++)
- {
- s += (tab7[i][0]* caleczka((((b-a)/2)*tab7[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 7 wynosi: "+ wynik);
- break;
- }
- case 8:
- {
- double tab8[][]={{0.36268378,-0.18343464},
- {0.36268378,0.18343464},
- {0.31370665,-0.52553241},
- {0.31370665,0.52553241},
- {0.22238103,-0.79666648},
- {0.22238103,0.79666648},
- {0.10122854,-0.96028986},
- {0.10122854,0.96028986}};
- /* double tab8[][]={{0.18343464,-0.18343464,0.52553241,-0.52553241,0.79666648,-0.79666648,0.96028986,-0.96028986},
- {0.36268378,-0.36268378,0.31370665,-0.31370665,0.22238103,-0.22238103,0.10122854,-0.10122854}};*/
- for (int i=0;i<tab8.length;i++)
- {
- s += (tab8[i][0]* caleczka((((b-a)/2)*tab8[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 8 wynosi: "+ wynik);
- break;
- }
- case 9:
- {
- double tab9[][]={{0.33023936,0.0},
- {0.31234708,-0.32425342},
- {0.31234708,0.32425342},
- {0.26061070,-0.61337143},
- {0.26061070,0.61337143},
- {0.18064816,-0.83603111},
- {0.18064816,0.83603111},
- {0.08127439,-0.96816024},
- {0.08127439,0.96816024}};
- /* double tab9[][]={{0.0,0.0,0.32425342,-0.32425342,0.61337143,-0.61337143,0.83603111,-0.83603111,0.96816024,-0.96816024},
- {0.33023936,-0.33023936,0.31234708,-0.31234708,0.26061070,-0.26061070,0.18064816,-0.18064816,0.08127439,-0.08127439}};*/
- for (int i=0;i<tab9.length;i++)
- {
- s += (tab9[i][0]* caleczka((((b-a)/2)*tab9[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 9 wynosi: "+ wynik);
- break;
- }
- case 10:
- {
- double tab10[][]={{0.29552422,-0.14887434},
- {0.29552422,0.14887434},
- {0.26926672,-0.43339539},
- {0.26926672,0.43339539},
- {0.21908636,-0.67940957},
- {0.21908636,0.67940957},
- {0.14945135,-0.86506337},
- {0.14945135,0.86506337},
- {0.06667134,-0.97390653},
- {0.06667134,0.97390653}};
- /*double tab10[][]={{0.14887434,-0.14887434,0.43339539,-0.43339539,0.67940957,-0.67940957,0.86506337,-0.86506337,0.97390653,-0.97390653},
- {0.29552422,-0.29552422,0.26926672,-0.26926672,0.21908636,-0.21908636,0.14945135,-0.14945135,0.06667134,-0.06667134}};*/
- for (int i=0;i<tab10.length;i++)
- {
- s += (tab10[i][0]* caleczka((((b-a)/2)*tab10[i][1])+((b+a)/2)));
- }
- wynik =((b-a)/2)*s;
- System.out.println("Wynik dla zadanego n wynoszacego 10 wynosi: "+ wynik);
- break;
- }
- }
- }
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement